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THE  DAMPING  OF  FLEXURAL  AND  ACOUSTIC  WAVES 
BY  A  BIAS-FLOW  PERFORATED  ELASTIC  PLATE 
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SUMMARY 

An  analysis  is  made  of  the  damping  of  sound  and  structural  vibrations  by 
vorticity  production  in  the  apertures  of  a  bias  flow,  perforated  elastic 
plate.  Unsteady  motion  causes  vorticity  to  be  generated  at  the  aperture 
edges;  the  vorticity  and  its  energy  are  swept  away  by  the  bias  flow  and  result 
in  a  net  loss  of  acoustic  and  vibrational  energy.  In  this  paper  we 
investigate  the  interaction  of  an  arbitrary  fluid- structure  disturbance  with  a 
small  circular  aperture  in  the  presence  of  a  high  Reynolds  number,  low  Mach 
number  bias  flow.  By  considering  the  limit  in  which  the  aperture  is  small 
compared  to  the  length  scale  of  the  impinging  disturbance,  it  is  shown  that 
the  effect  of  the  interaction  can  be  represented  by  a  concentrated  source  in 
the  plate  bending  wave  equation  consisting  of  a  delta  function  and  two  of  its 
axisymmetric  derivatives.  A  generalized  bending  wave  equation  is  then 
formulated  for  a  plate  perforated  with  an  homogeneous  distribution  of  small, 
bias  flow  circular  apertures.  This  equation  is  used  to  predict  the 
attenuation  of  sound  and  resonant  bending  waves  by  vorticity  production. 
Acoustic  damping  is  found  to  be  significant  provided  the  fluid  loading  is 
sufficiently  small  for  the  plate  to  be  regarded  as  rigid  (e.g.,  for  an 
aluminium  plate  in  air  when  the  frequency  is  not  too  small) .  On  the  other 
hand,  a  bending  wave  is  effectively  damped  only  when  the  fluid  loading  is 

enough  for  the  wave  to  produce  a  substantial  pressure  drop  across  the 
plate;  when  this  occurs  the  predicted  attenuations  are  comparable  with  those 
usually  achieved  by  the  application  of  elastomeric  damping  materials. 

Numerical  predictions  are  presented  for  steel  and  aluminium  plates  in  air  and 
water. 
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1.  INTRODUCTION 

Aerodynamic  sound  and  vibrations  are  generated  by  turbulence  and  other 
unsteady  distributions  of  vorticity,  and  by  their  interactions  with  adjacent 
elastic  bodies  [1  -  11].  Acoustic  energy  can  also  be  dissipated  by  vorticity 
production,  however,  when  a  sound  wave  impinges  on  a  solid  surface  or  fluid 
region  of  nonuniform  mean  density.  For  example,  sound  incident  on  the  sharp 
trailing  edge  of  an  airfoil  in  a  mean  flow  will  cause  vorticity  to  be  shed 
into  a  wake  in  accordance  with  the  unsteady  Kutta  condition;  at  low  Mach 
numbers  this  leads  to  an  increase  in  the  kinetic  energy  of  the  mean  flow  and 
wake  at  the  expense  of  the  sound  [12,  13] .  Similarly,  one  or  more  perforated 
plates  aligned  with  the  mean  flow  through  the  tube  banks  of  a  heat  exchanger 
cavity  can  be  used  to  suppress  flow- induced,  cavity  acoustic  resonances  by  the 
production  of  vorticity  in  the  perforates;  the  kinetic  energy  of  the  vorticity 
is  extracted  from  the  sound,  converted  away  by  the  mean  flow,  and  ultimately 
dissipated  by  viscous  and  thermal  processes  [14,  15] .  The  same  mechanism  is 
responsible  for  the  greatly  improved  attenuation  of  "screech"  tones  in  the  jet 
pipe  of  a  jet  engine  when  a  bias  flow  is  maintained  through  the  perforated 
heat  shield  normally  used  to  protect  the  wall  liner  from  contact  with  hot 
combustion  products  [16,  17]. 

The  effectiveness  of  the  acoustic  attenuation  in  these  cases  is  critically 
dependent  on  the  presence  of  mean  flow  [18  -  21] .  Analytical  predictions  for 
rigid  perforated  screens  indicate  that  the  attenuation  is  linearly 
proportional  to  the  acoustic  amplitude  [16,  17,  20  -  24],  which  is  in  broad 
agreement  with  experiment  [17  -  19,  24,  25].  In  the  absence  of  flow  the 
attenuation  is  either  a  nonlinear  function  of  the  acoustic  amplitude  (and 
therefore  weak)  or  is  dominated  by  the  less  efficient  action  of  viscous 
dissipation  in  the  apertures  [26  -  30]. 

Theoretical  analyses  of  sound  absorption  by  vorticity  production  have 
hitherto  assumed  the  screen,  plate,  etc  (where  the  vorticity  is  generated)  to 
be  rigid.  However,  in  most  practical  configurations  high  acoustic  intensities 
are  accompanied  by  structural  vibrations ,  and  a  substantial  part  of  the  noise 
energy  is  actually  contained  in  vibratory  modes  of  the  structure.  These  modes 
would  also  be  expected  to  be  damped  by  vorticity  production,  and  the  primary 
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practical  objective  may  be  the  suppression  of  potentially  harmful  structural 
vibrations  rather  than  the  sound.  The  simplest  model  problem  that  exhibits 
structural  damping  by  this  means  (and  which  is  analogous  to  the 
acoustic  -  trailing  edge  problem  mentioned  above  [12,  13])  is  that  in  which  a 
bending  wa.ve  on  a  thin  elastic  plate  is  reflected  at  a  trailing  edge  in  a 
parallel  mean  flow  (see  [31]).  Incident  bending  wave  energy  is  dissipated  at 
the  edge  by  scattering  into  sound,  and  by  a  transfer  to  the  kinetic  energy  of 
the  mean  flow  and  wake.  Below  the  coincidence  frequency  of  the  plate 
(conventionally  defined  as  the  frequency  above  which  the  phase  velocity  of 
bending  waves  on  the  plate  in  vacuo  exceeds  the  speed  of  sound  in  the  fluid) 
the  efficiency  with  which  sound  is  generated  is  usually  very  small ,  but  much 
larger  net  losses  are  predicted  in  [31]  owing  to  the  generation  of  wake 
vorticity . 

In  this  paper  we  investigate  the  efficiency  with  which  both  sound  and  thin 
plate  bending  waves  are  attenuated  when  the  plate  is  drilled  with  an 
homogeneous  array  of  small  circular  apertures  through  which  is  maintained  a 
steady  bias  flow.  Damping  is  caused  by  vorticity  generation  at  the  aperture 
edges  triggered  by  the  fluctuating  pressure  drop  across  the  plate  produced  by 
the  incident  disturbance.  The  attenuation  of  sound  by  a  rigid  bias  flow 
screen  was  considered  in  [21] ;  it  will  be  concluded  from  the  present  analysis 
that  the  rigid  plate  approximation  is  applicable  only  for  small  fluid  loading, 
for  a  steel  plate  in  air,  say.  A  sound  wave  incident  on  the  same  plate  in 
water  generally  induces  a  negligible  pressure  differential  across  the  plate, 
and  the  acoustic  attenuation  is  also  negligible.  For  heavy  fluid  loading, 
bending  waves  are  predicted  to  be  significantly  damped  over  an  intermediate 
range  of  low  frequencies,  centred  around  an  aperture  Strouhal  number  wR/U 
0(1),  where  w  denotes  radian  frequency,  R  aperture  radius,  and  U  is  the  mean 
velocity  in  the  aperture.  The  strength  of  this  damping  will  be  compared  with 
that  normally  achieved  by  coating  the  plate  with  an  elastomeric  material  [32] , 
whose  use  frequently  entails  substantial  and  undesirable  increases  in  the 
structural  mass,  especially  in  applications  at  low  frequencies.  It  is 
predicted  that  comparably  large  attenuations  are  possible  by  perforating  the 
plate  (say,  with  5%  "open  area")  and  maintaining  a  bias  flow  through  the 
apertures.  In  principle  the  bias  flow  velocity  and  aperture  radius  can  be 
chosen  to  yield  optimal  attenuation  over  a  prescribed  frequency  range,  and  a 
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practical  damping  system  might  be  "tuned"  within  certain  limits  by  varying  the 
bias  flow  pressure  to  accommodate  varying  frequency  characteristics  of  a 
particular  noise  source. 

The  scattering  theory  for  an  arbitrary  flexural  disturbance  incident  on  a 
single,  small  aperture  in  the  presence  of  a  bias  flow  is  developed  in  §2.  The 
result  is  used  (§3)  to  formulate  a  generalized  bending  wave  equation  for  a 
plate  perforated  with  a  uniform  (but  sparse)  distribution  of  apertures. 
Numerical  results  are  given  in  §3  for  the  attenuation  of  sound  and  bending 
waves,  and  a  comparison  made  with  damping  produced  by  an  elastomeric  coating. 
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2.  SCATTERING  BY  AN  APERTURE  IN  THE  PRESENCE  OF  A  BIAS  FLOW 
2 . 1  The  governing  equations 

A  thin  infinite  plate  of  thickness  h,  mass  m  per  unit  area  and  bending 
stiffness  B  is  immersed  in  fluid  of  mean  density  and  sound  speed  .  In 
the  undisturbed  state  the  plate  occupies  the  plane  =  0  of  the  rectangular 
coordinate  system  (x^jX^jX^)  and  is  pierced  by  a  small  circular  aperture  of 
radius  R  whose  center  coincides  with  the  coordinate  origin  (see  Figure  1).  A 
mean  pressure  differential  is  maintained  across  the  plate  which  produces  a 
nominally  steady,  low  Mach  number,  high  Reynolds  number  "jet"  of  fluid  through 
the  aperture  from  the  "upper"  (x^  >  0)  to  the  "lower"  (x^  <  0)  side. 

Consider  small  amplitude,  time  harmonic  motions  of  frequency  w  (>  0) 
proportional  to  e"^^^  .  The  exponential  time  factor  is  henceforth  suppressed, 
and  the  flexural  displacement  of  the  plate  (in  the  positive  x^ -direction)  and 
perturbation  pressure  are  respectively  denoted  by  f^(x^,x^)  and  p(x) .  On 
=  0,  (x^+Xg)^^^  >  R,  P  satisfy  the  bending  wave  equation  [33] 

(BV^-  +  [p]  =  0,  >  R,  (2.1) 

where,  [p]  =  p(x^ ,+0 .Xj ) -p(x^ , -0 .Xg ) .  (2.2) 

and  V'*  s  (S^/Sx^+d^/ax^  )^  is  the  biharmonic  operator.  Except  in  the  immediate 

2  i  3 

vicinity  of  the  aperture,  where  the  influence  of  the  jet  becomes  important, 
the  small  amplitude  pressure  fluctuations  may  be  assumed  to  satisfy  the 
acoustic  equation 

{V^+  K^jp  =0,  x^  J  0,  (2.3) 

where  /c  =  w/c  is  the  acoustic  wavenumber.  This  equation  ignores  the 

o  '  o 

influence  of  convection  by  the  residual  effects  of  mean  flow  at  large 
distances  from  the  aperture,  but  this  is  negligible  when  the  Mach  number  is 
small . 

As  X  ->  ±0,  (x^+x^)^^^  >  R,  the  x  -component  of  the  fluid  velocity  becomes 

2  13  2 

equal  to  so  that  for  small  amplitude  flexural  motions  of  the  plate  the 

pressure  p  and  are  related  by  the  linearized,  x^ -component  of  the  momentum 
equation  of  the  fluid: 
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1  ap 


^  ±0.  (x2+x2)l/2  > 


(2.4) 


For  an  aperture  whose  radius  R  is  much  smaller  than  the  length  scales  of 
the  flexural  and  fluid  motions,  the  motion  in  the  immediate  vicinity  of  the 
aperture  may  be  regarded  as  axisymmetric .  The  boundary  conditions  that  no 
force  and  moment  are  applied  at  the  edge  can  therefore  be  taken  in  the  form 
[34] 


where 
plate . 


(1-a)  ^ 

r  a 


5-. 


=  0 


V  r  H  (x2+x2)i/2  ^  R+0, 


(2.5) 


a^/ax^+a^/aXg ,  and  o  is  Poisson’s  ratio  of  the  material  of  the 


2.2  The  composite  bending  wave  equation 

The  domain  of  validity  of  the  bending  wave  equation  can  formally  be 
extended  to  include  the  region  occupied  by  the  aperture  by  first  introducing 
the  indicator  function 


f(x^,X3)  =  r(x3,X3)  -  R,  (2.6) 

and  defining  a  composite  fluid-structure  displacement  f  on  x^  =  0  by 

^(x^.Xj)  =  r+H(f-£)  +  r.H(-f-0,  €  ^  +0,  (2.7) 

where  H(x)  is  the  Heaviside  step  function  (=1,  0  according  as  x  ^  1) ,  and 
denotes  the  perturbation  of  the  x^ -component  of  the  fluid  displacement  within 
the  aperture.  Equation  (2.1)  is  then  multiplied  by  H(f-£)  and  rearranged  to 
yield,  as  e  ->  +0, 

{BV^-  mu;M{H(f)r+}  +  H(f)  [p]  =  BV^  {V  { r+VH(f ) ) }  +  BV^  {VH(f ) } 

+  BV{(v2r^)VH(f))  +  BVH(f)»V(v2c^).  (2.8) 

where  VH(f)  =  d(f)Vf,  and  6(f)  is  the  Dirac  delta  function,  which  vanishes 
except  at  the  rim  r  =  R  of  the  aperture ;  and  its  derivatives  on  the  right 
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hand  side  are  evaluated  at  r  =  R+0. 

Continuity  of  pressure  within  the  aperture  implies,  for  a  plate  of 
infinitesimal  thickness  (more  precisely  for  h  ^  R) ,  that  [p]  —  0  within  the 
aperture,  and  therefore  (as  c  +0) 

nico2}{H(-f)f  }  +H(-f)[p]  =  {BV^  -  }  {H( -f )  } .  (2.9) 

When  this  is  combined  with  equation  (2.8)  the  composite  displacement  is  seen 
to  satisfy 

{BV^-  +  [p]  =  BV2{V{r+VH(f)}}  +  BV^  {VH(f) }  +  BV  { (V^  )VH(f ) } 

+  BVH(f)»V(v2$-^)  ^  }  ,  (2.10) 

where  the  "source"  terms  on  the  right  hand  side  depend  only  on  conditions  at 
the  edge  of  the  aperture  and  on  the  perturbation  flux  through  the  aperture. 
This  equation  is  applicable  everywhere  on  =0,  including  the  region 
occupied  by  the  aperture. 

Consider  the  limiting  form  of  this  equation  when  R  becomes  small  compared 
to  the  dominant  length  scales  of  the  fluid  and  structural  motions,  i.e.,  for 
K  R  K  R  «  1,  where  K  =  (mw^/B)^^''  is  the  in  vacuo  bending  wavenumber.  In 
this  limit  the  dominant  components  of  the  scattered  fields,  and  therefore  of 
the  source  terms  in  (2.10),  must  be  axisymmetric .  As  R  0  the  sources  vanish 
everywhere  except  at  r  =  0,  and  must  ultimately  reduce  to  combinations  of  the 
surface  delta  function  6(x)  =  and  its  axisymmetric  derivatives 

[35],  where,  however,  the  highest  order  derivatives  in  the  limiting  forms  of 
the  first  four  source  terms  in  (2.10)  must  be  restricted  to  ensure  that  the 
displacement  of  the  plate  remains  finite  as  r  R+0.  In  the  final  source 
{H(-f)f  }  ^  q6(x)  as  R  ->  0,  where  q  denotes  the  perturbation  volume 
displacement  through  the  aperture.  Thus,  in  the  leading  order  of 
approximation,  equation  (2.10)  becomes, 

{BV^-  ma)^}r  +  [p]  =  ajj6(x)  +  a^V\6{yi.)  +  qiBV^-  mw^}6(x),  R  -»  0,  (2.11) 

where  the  values  of  the  coefficients  a^ ,  a^ ,  q  remain  to  be  determined;  they 
represent  the  strengths  of  the  various  source  terms  to  which  the  aperture  is 
equivalent  as  R  ^  0. 
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2.3  The  scattering  problem 

Let  the  flexural  disturbance 

r  =  r  =  r  e^^^i  /cR  1 ,  (2.12) 

propagating  parallel  to  the  -axis  be  incident  on  the  aperture,  where  is 
constant.  and  the  accompanying  pressure  jump  across  the  plate  [p^ ] ,  say, 

satisfy  the  homogeneous  form  of  equation  (2.11)  (no  source  terms). 


Denote  the  scattered  flexural  displacement  and  pressure  by 
p  (x)  ,  respectively.  Equation  (2,11)  determines  the  functional  forms  of 
and  p  in  an  "outer  region"  at  distances  |x|  >  R  from  the  aperture  in  terms  of 
the  coefficients  a^ ,  a^ ,  q.  Write 


=J  fg .  k-  (k^.O.kg) 

~00 

/  ^  •  2  /  ^f  „i{k»x+7(k) |x,  I 

p^(x)  =  -ip^a;^sgn(x2)J  dk^ 


dk3  ,  X3  <  0,  J 


(2.13a,  b) 


where  f  (k)  is  the  Fourier  space  transform  of  with  respect  to  the  planar 
coordinates  x^  ,  Xg  ,  and  7(k)  =  (k  =lk|)  is  positive  for  real  k  <  , 

and  positive  imaginary  for  k  >  .  The  integral  representation  of  the 

scattered  pressure  is  a  solution  with  "outgoing"  wave  behavior  of  the  acoustic 
equation  (2.3),  and  is  applicable  in  the  outer  region  where  convection  by  the 
bias  flow  jet  can  be  ignored;  the  integrand  is  chosen  to  ensure  that  and  p^ 
satisfy  condition  (2.4). 

The  Fourier  amplitude  f ^ (k)  is  determined  by  substituting  from  (2.13) 
into  the  left  hand  side  of  (2.11).  The  net  displacement  on  x^  =  0  can  then  be 
cast  in  the  form: 


1  r“ 

"  (27r)2Ja 


{aj,  -a2k^+q(Bk^  -xau)  )  )e^ 
D(k,w) 


dkidk3  +  , 


(2.14) 


where 


D(k,w)  =  Bk^  -  mo)^  -  2ip^w^/7(k) 


(2.15) 
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is  the  bending  wave  dispersion  function  [33]  .  In  the  absence  of  dissipation 
D(k,cj)  has  two  real  zeros  that  determine  the  wavenumbers  k  of  undamped 
flexural  motions  of  the  homogeneous  plate ,  and  the  integration  contours  in 
(2.14)  must  be  indented  to  pass  around  the  corresponding  poles  of  the 
integrand  in  directions  that  are  consistent  with  the  radiation  condition. 
This  can  de  done  by  temporarily  assigning  to  u>  a  small  positive  imaginary 
part,  thereby  displacing  the  poles  off  the  real  k-axes  [36]. 


The  values  of  the  source  strengths  a^ ,  a^ ,  q  are  found  from  a 
consideration  of  conditions  in  the  immediate  vicinity  of  the  aperture.  This 
is  facilitated  by  first  re-casting  (2.14)  in  the  form 


1  r“ 

r(x,,=.3)  -  q«(*)  +  ^ 


{a(j  -a2k^+2iqp,,w^/7(k)  }e 
D(k,t^) 


ik»x 


dk,  dk,  +  r  e^'^^i. 

13  o 


(2.16) 


In  this  "outer"  formula  for  C ,  the  first  term  on  the  right  hand  side  is  the 
singular  component  of  f  associated  with  the  volume  displacement  through  the 
aperture.  The  remaining  terms  (which  are  non-singular  at  x  =  0)  make  finite 
contributions  to  the  plate  displacement  ,  and  two  linear  algebraic  equations 
involving  the  coefficients  a^ ,  a^ ,  q  are  obtained  by  substituting  for  into 
the  edge  conditions  (2.5).  Since  only  the  axisymmetric  part  of  the  scattered 
field  is  being  considered,  it  is  necessary  to  retain  only  the  axisymmetric 
component  f^Jp(/cr)  of  the  incident  wave  Bessel  function  expansion  [37] 

^  («r)  +  i"J^  (/cr)cos(no]  ,  x^  =  rcos^  . 

Conditions  (2.5)  are  then  found  to  give,  respectively. 


=  {kRJ^(kR)  -(l-a)J^(/cR))(/cryR) 

^02^0  +  ^2^2  +  ^2^  =  -K'j,(/cR)r„  . 

where  A„.  ,  A,.  ,  Q.  (j  =  1,  2)  are  defined  in  the  appendix. 

0  J  J 

A  third  relation  involving  ,  a^ ,  q  is  obtained  by  examining  the  motion 
in  an  "inner"  neighborhood  of  the  aperture.  The  pressure  here  can  be  written, 

p(x)  =  p^w^JT^x^  +  P’(x),  (2.18) 


(2.17) 
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where  the  first  term  on  the  right  is  defined  such  that  the  corresponding 
component  of  (-1/p^ )5p/ax2  at  =  ±0  is  equal  to  the  normal  acceleration  of 
the  plate  produced  by  the  incident  wave  When  the  characteristic 

bending  wavelength  is  much  larger  than  the  aperture  diameter  (K^R  «  1),  the 
motion  of  the  plate  at  distances  from  the  aperture  smaller  than  the 
characteristic  flexural  wavelength  will  be  identical  with  that  of  a  rigid 
plane  executing  translational  oscillations  of  amplitude  in  the  normal 
direction.  Since  also  k^R  <  1,  the  fluid  motion  may  similarly  be  regarded  as 
locally  incompressible,  and  the  additional  pressure  p’ (x)  in  (2.18)  calculated 
as  for  an  incompressible  fluid  subject  to  the  rigid  surface  condition 
d-p’ /dx^  =  0  on  x^  =0. 

The  volume  displacement  through  the  aperture  can  then  be  expressed  in  the 

form, 

q  =  [P„]/Po‘^^  •  (2.19) 

In  this  expression  is  the  Rayleigh  conductivity  of  an  aperture  in  a  rigid 
plate  taking  account  of  the  mean  bias  flow  [38];  it  has  the  dimensions  of 
length  and  typically  -  0(R).  p„  denotes  the  pressure  at  distances  lx|  from 
the  aperture  satisfying  R  <  |x|  <  .  where  the  flow  may  still  be  regarded 

as  incompressible.  This  interval  in  |x|  will  exist  provided  the  frequency  is 
sufficiently  small,  and  it  then  defines  an  "overlap"  region  in  which  both  the 
incompressible  representation  (2.18)  and  the  outer  solution  (2.13b)  are  both 
applicable.  This  implies  that  [p_]  can  be  calculated  by  considering  the 
limiting  value  of  the  integral  in  (2.13b)  in  the  overlap  region.  This 
limiting  value  contains  the  term 

-p^w2q/27r|xl  =  [p„  ]/2  |x|  <  [p„  ]  ,  (2.20) 

which  is  discarded,  and  we  then  obtain  (including  also  the  contribution  [Pj ]  • 
of  the  incident  disturbance) , 


a 

[pj  “  [Pii  -  + 


[ap  -a2k^+2iPotJ^  q/7(k)  ]kdk-> 

7(k)D(k,a))  J' 


(2.21) 


For  consistency  with  the  matching  procedure  the  remaining  terms  in  q  on  the 
right  hand  side  should  also  be  discarded,  since  they  are  easily  shown  to  be  at 
most  0(K  Kj^)[p^],  which  is  of  smaller  order  than  the  term  discarded  in  (2.20). 
We  shall  do  this,  but  it  should  be  noted  that  this  is  actually  equivalent  to 
neglecting  a  small  effect  of  fluid  compressibility  responsible  for  energy 
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losses  due  to  acoustic  scattering.  Then  substituting  into  the  right  hand  side 
of  (2.19),  we  obtain  the  third  equation  for  a^ ,  a^ ,  q: 


^03^0  +  A23a3  +  q  =  \  [Vj]/Py  . 


(2.22) 


where  and  are  defined  in  the  Appendix. 


The  solutions  of  equations  (2.17),  (2.22)  are  now  expanded  in  powers  of 
K  R  («^  1)  .  Only  the  axisynunetric  components  of  the  scattered  fields  are  being 
retained  in  the  present  approximation,  and  it  can  readily  be  verified  that  the 
expansion  must  therefore  stop  at  terms  of  order  (K^R)  .  The  expansion  can  be 

effected  by  first  nondimens ionali zing  the  expressions  (A.l)  -  (A. 3)  for  the 
coefficients  A^^  ,  A^^  >  found  that  a^  and  a^  are  each  0{(K^R)  }  and 

that  q  -  O(K^R).  To  this  order  of  approximation  equations  (2.17)  and  (2.22) 
yield 


a^  =  -7rR^B/c^r„  '  2RF^[Pi] 


a  =  -  - ttR  Bk  ^ 

2  (1-a) 


q  = 


p  w 


[Pt  ]  > 


,2 


(2.23) 


In  these  expressions 


{BV^-  mw2}r  +  [p] 


kV  =  so  that  (2.11)  becomes 

2^1  o  21 

=  -{^rR^BV^fj  +  2RKjj[Pj]}6(x)  +  •V2'5(x) 


+  -!^[p,  ]{BV^-  mw2}«5(x).  (2.24) 

Po’^ 

In  this  form,  the  equation  describes  the  scattering  of  an  arbitrary  incident 
displacement  and  pressure  field  p^  whose  characteristic  length  scale  is 
much  larger  than  the  aperture  radius. 

2.4  Scattering  of  a  bending  wave 

The  physical  significance  of  the  different  source  terms  on  the  right  of 
equation  (2.24)  can  be  exhibited  by  consideration  of  the  problem  in  which  the 
incident  disturbance  ^ is  a  plane  bending  wave,  for  which, 
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^  sgn(x,)ip„a>^ro^i{«x,-H7(«)|xJ}  ^  ^  ,  (2.25) 

Pi  7(«) 

where  the  wavenumber  k  >  may  be  taken  to  be  the  positive  real  zero  of  the 
dispersion  function  D(k,w)  defined  by  (2.15).  This  is  an  evanescent  wave  that 
propagates  subsonically  relative  to  the  fluid;  the  accompanying  fluid  motion 
is  confined  to  boundary  layer  regions  of  width  1/ (/c  -k^  )  ^  on  either  side  of 

the  plate  within  which  the  pressure  p^  decays  exponentially,  although  the 
decay  may  be  slow  when  k  -  (i.e. ,  when  w  >  ) . 

The  interaction  at  the  aperture  will  result  in  an  exchange  of  energy 
between  the  wave  field  and  the  bias  flow  jet,  whose  magnitude  can  be 
determined  from  a  consideration  of  the  energy  equation.  To  do  this,  consider 
an  homogeneous,  time - dependent ,  fluid  loaded  plate,  driven  by  a  distribution 
F(Xj^,X2,t)  of  normal  stresses  (taken  to  be  positive  in  the  x^ -direction)  . 

This  equation  is  derived  from  (2.1)  by  first  replacing  w  by  id/dt,  and 
inserting  F  on  the  right  hand  side.  Multiplying  by  f  s  d^/dt  and  rearranging, 
we  find 

+  |B(v2c)2j  +  div[BfV(v25:)  -  B(v2f)Vfj  +  r[p]  =  TF.  (2.26) 

'j'M0  first  term  on  the  left  is  the  time  rate  of  change  of  the  structural  energy 
density;  the  divergence  describes  the  fltix  of  energy  within  the  plate,  and 
f[p]  is  the  energy  flux  from  the  plate  into  the  fluid.  The  power  supplied  or 
dissipated  by  the  stress  field  F  is  determined  by  their  rate  of  working  fF. 

In  equation  (2.10)  the  role  of  the  stress  distribution  F  is  taken  by  the 
limiting  form  of  the  sources  as  e  ^  0,  where  the  composite  displacement  f  is 
given  by  (2.7).  To  evaluate  the  mean  energy  flux  from  the  aperture  the  time 
factor  is  restored  and  the  real  parts  taken  of  the  resulting  expressions 

for  all  perturbation  quantities.  The  sources  are  then  multiplied  by 

a$:/3t  =  Re{-iwr(x^  ,X3)e'^“''}  s  Re{ -iwf  ^  (x^  ,X3  )e‘^“''  -  iwr3  (^i -Xg  )  . 

averaged  over  a  wave  period  2n/o>,  and  integrated  with  respect  to  (x^^  ’X3)  • 

Note  that  the  first  four  source  terms  on  the  right  of  (2.10)  are  multiplied  by 
the  component  ar+/3t  of  d^/dt,  and  the  final  source  by  d^_/dt.  As  c  ^  +0  this 
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is  equivalent  to  multiplying  the  first  two  sources  in  (2.11)  by  3(*^/5t  and  the 
final  source  by  If  denotes  the  incident  bending  wave  power 

dissipated  at  the  aperture,  we  then  find 

n  =  -  .  R(/c2-K2)^/2]Im{IC  }  .  (2.27) 

A  (/c2-/c2)  ° 

Thus,  flexural  wave  power  is  dissipated  at  the  aperture  provided  the 
imaginary  part  of  the  conductivity  is  negative  (for  /cR  1).  In  the 
absence  of  the  bias  flow,  and  when  viscous  dissipation  is  ignored,  is  real 
and  equal  to  2R  for  a  plate  of  negligible  thickness  [38],  In  this  case, 
however,  energy  must  actually  be  extracted  from  the  bending  wave  by  scattering 
into  sound.  This  loss  cannot  appear  in  (2.27)  because  the  left  hand  side  of 
the  energy  equation  already  accounts  for  energy  interchanges  between  the 
structural  and  acoustic  fields. 


The  aperture  conductivity  becomes  complex  when  account  is  taken  of  the 
bias  flow.  For  a  high  Reynolds  number,  low  Mach  number  jet  it  can  be 
expressed  in  terms  of  the  aperture  Strouhal  number  =  |w|R/U^  in  the  form 


Kjj  =  2R{r(S^)  -  iA(S^)), 


(2.28) 


where  T  and  A  are  real  valued,  and  is  the  mean  convection  velocity  of 
vorticity  in  the  jet  shear  layer  (approximately  equal  to  one  half  of  the  jet 
velocity).  The  following  expressions  for  T  and  A  are  derived  in  [21]  and  used 
in  the  following  discussion: 


TT^  [1+S^  ]I^  (S^.  )^  +  4[S^cosh(S^)  -  sinh(Sj.  )]cosh(Sj.  )K^  (S^.  )^e^^t 

(S,.  )^  +  Acosh^  (S^.  )K^  (S^,  )^e^^t] 


r(s^)  = 


27rI^(S^)K^(S^)e2St 

A(S  )  =  — - - - „ 

SjTr^I^CS^)^  +  4cosh2(S^)K^(S^)2e^\] 


(2.29) 


where  I  and  K,  are  modified  Bessel  functions  [37].  r(S  )  and  A(S  )  are  both 
positive  for  >  0,  and  their  variations  with  are  illustrated  in  Figure  2. 
Predictions  based  on  these  formulae  of  the  attenuation  of  sound  by  a  rigid 
bias  flow  perforated  screen  have  been  found  to  be  in  good  agreement  with 
experiment  [17] . 
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The  first  term  in  the  square  brackets  of  (2.27)  is  supplied  by  the  volume 
displacement  (final)  source  in  equation  (2.24);  it  accounts  for  the  absorption 
of  energy  by  the  mean  flow.  The  second  term,  which  is  of  opposite  sign,  comes 
from  the  first  source  on  the  right  of  (2.24)  (i.e.,  from  a^ ) ,  and  describes 
the  piTodticdoTi  of  bending  waves  at  the  aperture  by  the  extraction  of  energy 
from  the  mean  jet  flow.  However,  /cR  ^  1,  and  the  net  effect  is  always  energy 
dissipation,  at  least  for  the  range  of  wavenumbers  and  frequencies  for  which 
the  present  approximations  are  applicable. 

The  absorption  cross-section  of  the  aperture  is  equal  to  n^/2RIIj  ,  where 
IIj  is  the  power  in  the  incident  bending  wave  per  unit  distance  parallel  to  the 
plate  and  wavefront.  It  can  be  expressed  in  the  form  [39] 


n 


I 


zSD.  , 
— (k.w) , 

OK 


(2.30) 


where  D(k,w)  is  the  bending  wave  dispersion  function  (2.15).  The  energy  in 
the  incident  wave  is  shared  between  the  kinetic  and  elastic  energies  of  the 
plate  and  the  energy  of  the  evanescent  fluid  motions  on  either  side  of  the 
plate  [33].  Combining  equations  (2.27)  and  (2.30)  we  find 

8(ep„R/P3h)^^2  A(S,)[1  - 

2BIl^  (MS^)1/2  (/c/K^)[5(/c/K^)''  -  4(K/cyK2)2  .jj  ’ 


where  M  =  U  /c  is  the  convection  Mach  number  of  the  bias  flow  jet,  p  =  m/h 

C  '  O 

is  the  mass  density  of  the  material  of  the  plate,  and  6  is  a  fluid  loading 
parameter  defined  by 

.  -  ^ _ ? — 

^s^l2p^cZ(l-a2) 

where  E  is  Young’s  modulus  for  the  plate. 

The  dependence  of  n^/2Rnj  on  Strouhal  number  =  uR/U^  is  illustrated  in 
Figure  3  for  an  aperture  in  a  steel  plate  in  water  (e  »  0.135)  when  the  bias 
flow  convection  Mach  number  M  =  0.003  (U^  =4.5  m/s).  The  absorption  peaks 
near  =  1  where  A(S^.)  s  -Im{Kj^}  assumes  its  maximum  value.  It  may  be 
remarked  that  the  corresponding  loss  due  to  acoustic  scattering  (which  has 
been  neglected)  is  O(k^R)  relative  to  (2.31). 


(2.32) 
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3.  THE  BIAS -FLOW  PERFORATED  ELASTIC  PLATE 
3.1  The  generalized  bending  wave  equation 

Consider  next  the  case  of  an  elastic  plate  perforated  with  a  distribution 
of  small  circular  apertures.  Let  the  apertures  have  equal  radius  R,  and  be 
distributed  with  their  centers  at  the  random  points  =  (x^^  ’®>^n3  ^ 
plate,  with  N  per  unit  area.  By  setting 

we  can  define  the  following  analog  of  equation  (2.7) 

fCx^.Xj)  =  []„  H(f„-£)  +  r.  H(-f^-0.  £  ^  +0.  (3.1) 

When  the  open  area  ratio  a  =  NttR^  «  1,  the  distance  (-  1/VN)  between 
neighboring  apertures  is  much  greater  than  R.  Then,  provided  k^^R.  1> 

the  length  scales  of  variation  of  all  disturbances  incident  on  an  aperture 
(including  those  -  0(1/VN)  due  to  scattering  by  neighboring  apertures)  will  be 
large  relative  to  R,  and  the  dominant  component  of  the  scattered  acoustic  and 
flexural  waves  will  have  radial  symmetry  with  respect  to  that  aperture.  To 
the  order  of  approximation  considered  in  §2,  it  can  then  be  deduced  that 

{BV^-  +  [p]  =  2RKj^[pJ)6(x-xJ 


(l+u) 

(1-u) 


.^R2BV2f^.v2  6(x-x^) 


+  ^^[p„]{BV^-  mw^}6(x-x  )1. 

n  0)^  n  z  J 


(3.2) 


In  this  equation  and  [p^]  respectively  denote  the  flexural  displacement  and 
pressure  jump  at  x  produced  by  any  incident  sound  or  flexural  waves  together 
with  collective  contributions  caused  by  scattering  by  the  remaining  apertures. 

When  the  characteristic  wavelengths  of  the  sound  and  structural  motions 
are  large  compared  to  the  distance  between  neighboring  apertures  it  is 
convenient  to  simplify  (3.2)  by  averaging  with  respect  to  aperture  positions. 
This  enables  attention  to  be  focussed  on  the  overall  effect  of  the  apertures 
rather  than  local  details  at  each  aperture,  and  henceforth  ("  and  p  will  be 
taken  to  denote  the  averaged  displacement  and  pressure  respectively.  This 
averaging  procedure  is  equivalent  to  that  used  to  study  wave  propagation  in 
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dusty  gases  and  bubbly  media  [40  -  45] ,  and  in  the  classical  Lorentz- Lorenz 
theory  of  a  dielectric  [46].  To  determine  the  average  contribution  from  the 
sources  in  (3.2)  at  x  =  ,  say,  we  first  obtain  appropriate  expressions  for 

r  and  [p  1  by  means  of  the  following  argument.  For  a  large  plate  of  area  A 
the  probability  that  x^  lies  within  the  area  element  dx^  =  dx^^^  dx^^  is  equal 
to  dx  /A,  and  there  are  NA  apertures  in  all.  In  a  first  approximation,  for 
a  <1,  the  apertures  may  be  regarded  as  independently  distributed,  and  the 
average  on  the  right  of  (3.2)  can  be  evaluated  by  first  averaging  over  all 
X  ^  JL  .  The  quantities  r  and  [p  ]  must  then  be  set  equal  to  the 
corresponding  averaged  incident  fields,  following  which  we  can  average  with 
respect  to  x^ .  When  the  total  number  of  apertures  is  large,  these  incident 
field  averages  cannot  differ  significantly  from  the  overall  average 
displacement  and  pressure  jump,  since  the  latter  are  hardly  affected  by  the 
presence  or  absence  of  the  aperture  at  x^  when  NA  Thus  in  the  leading 

approximation  the  contribution  from  x^  as  NA  -*•  “  is  given  by  writing 


[pj  =  [p(x^^  .0,x^3  )],  for  a  <  1. 

Making  this  substitution  in  (3.2)  and  performing  the  average  with  respect 

to  the  X  ,  we  obtain  the  generalized  bending  wave  equation 
n 


ib  ■  ^  (l  *  2"^ [l  ■  ■  ”“')]] 'Pi  - 


In  this  equation  the  term  in  a  accounts  for  an  effective  decrease  in  the 
stiffness  of  the  plate  caused  by  the  presence  of  the  apertures:  this  is 
independent  of  fluid  loading,  and  an  identical  term  is  obtained  if  the 
calculation  is  repeated  for  a  plate  in  vacuo.  The  terms  involving  the 
conductivity  represent  the  influence  of  the  fluid  at  the  apertures.  In  the 
absence  of  bias  flow  all  of  the  modified  coefficients  in  the  equation  are 
real,  and  equation  (3.3)  then  has  undamped  plane  wave  solutions  that 
correspond  to  the  undamped  bending  wave  modes  of  an  ideal,  fluid  loaded, 
homogeneous  plate  (the  attenuation  caused  by  acoustic  scattering  at  the 
apertures  is  a  higher  order  effect,  and  is  not  included  at  the  present  order 
of  approximation) . 
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3.2  Dissipation  of  sound 

To  illustrate  the  implications  of  equation  (3.3)  consider  the  absorption 
of  sound  incident  on  a  bias  flow  elastic  screen.  For  simplicity  attention  is 
confined  to  plane  acoustic  pressure  waves  of  the  form  e'^'^o^z  incident 
norimlly  on  the  screen  from  >  0  (see  Figure  4) .  The  pressure  can  be  taken 
in  the  form 


p  =  e*^^o^2  +  Re ^^0^2  ,  >  0, 

=  3e"^**'o^2  ,  <  0, 


(3.4) 


where  R  and  3  are  appropriate  reflection  and  transmission  coefficients, 
respectively.  The  (generalized)  displacement  C  of  the  screen  is  related  to  p 
by  $•  =  (1/p^w^  )3p/5x2  (x^  -»  ±0).  This  relation  and  equation  (3.3)  supply  two 
equations,  from  which  R  and  3  are  found  to  be  given  by 


R  =  1  -  3 


nuc  /p 

o'  ^  C 


nwc 


.■  2i(i  +  |.(r-iA)[i  + 


(3.5) 


Define  an  absorption  coefficient  611^/11^  ,  where  denotes  the  acoustic 
power  incident  on  the  screen  and  611^  is  the  dissipated  power.  Then 

-  1  -  |R|^  -  |3|^ 


^MS  aTi  +  2^ 

TT  t  1_  P  hj 


(3.6) 


The  solid  curves  in  Figure  5  depict  the  variation  of  611^/11^  with  aperture 
Strouhal  number  =  wR/U^  for  sound  incident  on  a  perforated  aluminium 
plate  in  air,  when  the  bias  flow  convection  Mach  number  M  =  0.05  and  when  R/h 
=  5.  According  to  this  figure  the  maximum  possible  attenuation  of  611^/11^  just 
less  than  0.5  (3  dB)  occurs  at  an  open  area  ratio  a  =  M  over  a  range  of 
Strouhal  numbers  centered  on  S^.  »  0.4.  At  very  low  frequencies  the  impedance 
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offered  to  the  sound  by  the  plate  becomes  negligible,  the  bias  flow  jet 
velocities  are  unmodulated  by  the  sound,  and  there  is  no  transfer  of  acoustic 
energy  to  hydrodynamic  motions.  In  this  low  frequency  limit  (S^  <  1) 


which  implies  that  significant  damping  can  occur  provided  wh/c^  >  p^/p^  ■  The 
results  shown  in  the  figure  are  for  R/h  =  5.  The  present  theory  is  applicable 
only  when  R  is  large  enough  to  ensure  that  the  Reynolds  number  of  the  bias 
flow  jets  is  large.  At  smaller  values  of  R/h  the  predicted  absorption  is 
similar  to  that  shown  in  the  figure,  but  the  widths  of  the  bell- shaped  solid 
curves  increase,  and  the  maxima  shift  to  slightly  lower  frequencies. 


The  corresponding  attenuation  produced  by  a  bias  flow  rigid  screen  at 
arbitrary  frequency  is  obtained  by  setting  /j^R/p^h  =  0  in  (3.6);  the 
predictions  in  this  case  are  depicted  by  the  broken  curves  in  Figure  5.  At 
low  frequencies  the  rigid  screen  absorption  coefficient  approaches  the 
limiting  value  given  by  (3.7)  when  p^R/p^h.  ■*  0  namely,  611^ /II^  =  2aM/(M+a)^ 

(S^  1)  [13]  ,  which  takes  its  maximum  value  (of  |)  when  a  =  M,  as  for  the 

elastic  screen.  The  elastic  screen  predictions  of  Figure  5  differ  from  those 
for  a  rigid  screen  when  the  frequency  is  low  enough  that  the  screen  may  be 
regarded  as  acoustically  transparent.  This  behavior  is  typical  of  situations 
in  which  the  fluid  mass  loading  of  the  plate  is  small  (e.g.,  for  metallic 
screens  in  air).  For  large  mass  loading  (which  in  practice  occurs  when  p^/p^ 
>  10’^  ,  for  steel  in  water,  say)  equation  (3.6)  predicts  that  the  absorption 
of  acoustic  energy  by  the  elastic  screen  is  negligible  at  all  frequencies,  in 
contrast  to  the  rigid  screen  predictions,  which  are  similar  to  those 
illustrated  by  the  broken  curves  in  Figure  5. 
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3.3  Damping  of  flexural  waves 

The  damping  of  a  flexural  wave  can  be  examined  similarly.  Let  the 
flexural  wave  be  the  plane  bending  wave  f  with  the  associated 

pressure  field  defined  as  in  equation  (2.25)  where  now,  however,  the  bending 
wavenumber  /c  is  a  zero  of  the  generalized  dispersion  function 

D(/c,cc)) 


2aa 

(T^- 


mccj 


7(«) 


^1  +  2NRKjj  1 


2aa 

(1-ct). 


B/c'’ 


which  is  obtained  by  substituting  for  f  and  [p]  in  equation  (3.3). 


""']]]■ 

(3.8) 


In  the  absence  of  the  bias  flow  (K^  ^  2R)  ,  D(/c,C4;)  has  precisely  two  equal 
and  opposite  real  zeros.  They  satisfy  \k\  >  (where  7(/c)  =  i (/c^ -/c^ ) and 
correspond  to  undamped  flexural  waves  on  the  perforated  plate;  they  are 
generally  close  to  the  corresponding  zeros  of  the  dispersion  function  D(k,(i;) 
of  the  homogeneous  plate  (see  (2.15)).  To  the  present  order  of  approximation 
there  is  no  damping  due  to  scattering  of  structural  energy  into  sound.  This 
is  analogous  to  the  corresponding  approximation  in  the  theory  of 
surface-acoustic  waves  propagating  over  a  nominally  plane  rigid  wall  with 
small  surface  irregularities,  where  the  scattering  by  surface  roughness 
elements  of  characteristic  dimension  R  causes  the  wave  to  be  damped  only  over 
very  large  propagation  distances  of  order  1/{R^k:^}  [47,  48]. 


The  absorption  of  bending  wave  energy  by  the  bias  flow  jets  causes  the 
formerly  real  bending  wavenumber  k  to  acquire  a  small  positive  imaginary  part 
that  accounts  for  the  progressive  decay  of  the  wave.  The  reduction  in  the 
flexural  wave  power  over  a  propagation  distance  is  equal  to 

20Im{/c)6x^iog^Q  (e)  dB,  and  the  power  absorbed  per  wavelength  of  propagation  is 
therefore  (e)Im{K}/Re{/c)  «  54. 6Im{K}/Re{Ac}  dB.  The  dependence  of  this 

absorption  on  the  bias  flow  Strouhal  number  is  illustrated  in 

Figures  6  and  7  for  a  perforated  steel  plate  in  water  and  several  different 
open  area  ratios  a  for  an  aperture  radius  R  =  h,  3h  respectively.  In  both 
cases  the  bias  flow  vorticity  convection  velocity  =  3  m/s  (M  ==  0.002). 
Increasing  the  open  area  ratio  increases  the  attenuation,  whereas  increasing 
the  size  of  the  apertures  reduces  the  maximum  possible  attenuation.  These 
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results  have  been  obtained  from  a  numerical  solution  of  the  dispersion 
equation  D(/c,w)  =  0,  by  iterating  about  the  real  zero  k  in  the  absence  of  the 
bias  flow  jets.  The  absorption  may  also  be  estimated  from  the  absorption 
cross-section  formula  (2.31).  The  fraction  of  the  incident  power  absorbed 
in  a  distance  in  the  direction  of  propagation  is  =  2Im{/c}6Xj^  ; 

using  (2.31)  it  follows  that 

Im{Ac}  ^  8a6A(St){l  -  R(«^-«g)^^^} _ 

Re{«;)  ttMS^  (k/K^  )2  [5(«/K^  )'*  -  4(«;«;^/K2)2  -  1] 

where  on  the  right  hand  side  k  denotes  the  real  value  of  the  wavenumber  in  the 
absence  of  dissipation.  This  approximate  formula  tends  to  over  predict  the 
attenuation  by  several  dB  near  the  maxima  of  Figures  6  and  7.  It  shows, 
however,  that  the  change  in  the  attenuation  with  R/h  (for  fixed  values  of  the 
remaining  parameters)  is  caused  by  the  consequent  change  in  the  real  part  of 
the  bending  wavenumber  k  relative  to  the  vacuum  wavenumber  (c.f.,  Figure  8 
below).  Note  also  that,  implicit  in  the  present  analysis  is  the  requirement 
that  R  >  h,  and  this  might  cast  some  doubt  on  the  validity  of  the  results 
shown  for  R/h  -  1  in  Figure  6.  However,  the  damping  is  governed  by  the 
imaginary  part  (-A)  of  the  Rayleigh  conductivity  ,  which  is  non-zero  because 
of  the  production  of  vorticity  in  the  apertures.  The  value  of  A  given  in 
(2.29)  would  be  expected  to  be  applicable  provided  the  length  scale  of 

the  vorticity  is  large  compared  to  the  plate  thickness  h,  and  the  predictions 
of  Figure  6  should  therefore  be  valid  at  least  for  <  1 . 

To  assess  the  likely  significance  of  the  attenuations  produced  by  a  bias 
flow  screen  we  plot  in  Figure  8  the  corresponding  absorption  of  bending  waves 
on  a  coated  steel  plate  in  water.  For  the  purpose  of  illustration  the  mass  of 
the  coating  is  neglected,  but  the  coating  is  assumed  to  make  the  effective 
stiffness  of  the  plate  complex,  such  that  the  bending  stiffness  B  of  the 
uncoated  plate  is  replaced  by  B(1  -  ii])  ,  where  rj  (>  0)  is  a  loss  factor  whose 
value  is  determined  by  the  properties  of  the  coating,  and  is  typically  of  the 
same  order  as  shown  in  the  figure  [32].  The  abscissa  ^ ’ 

where  ==  c^(m/B)^^^  is  the  coincidence  frequency.  As  w  increases  towards 
both  the  bending  and  acoustic  wavelengths  become  comparable  to  the  thickness 
of  the  plate  and  the  thin  plate  equation  (2.1),  on  which  the  present  analysis 
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is  based,  ceases  to  be  applicable.  A  comparison  with  figures  6  and  7 
indicates  that  the  attenuations  achieved  by  the  bias  flow  jets  can  be 
comparable  with  those  predicted  for  the  coated  plate. 

Much  smaller  attenuations  are  predicted  for  bending  waves  propagating  on 
a  bias  flow  screen  in  air.  Figure  9  depicts  typical  predictions  for  an 
aluminium  plate  with  three  different  open  area  ratios  when  R  =  5h  and  the  mean 
bias  flow  jet  velocity  U  =  2U^  =  3.4  m/s.  In  possible  applications  the 
thickness  h  of  the  plate  will  usually  be  small  (of  the  order  of  0.01  cms) ,  and 
the  bias  flow  Reynolds  number  UD/v  (D  =  2R  being  the  aperture  diameter,  and  v 
the  kinematic  viscosity)  must  be  large  enough  to  ensure  that  motion  in  the 
apertures  is  not  impeded  by  viscous  action.  For  the  cases  illustrated  in  the 
figure  UD/i/  =  250,  which  is  probably  close  to  the  lower  limit  for  which  the 
theory  is  likely  to  be  relevant.  The  Reynolds  number  could  be  raised  by 
operating  at  a  higher  jet  velocity,  but  this  would  degrade  the  absorption 
which,  according  to  (3.9),  is  inversely  proportional  to  the  vorticity 
convection  Mach  number  M. 
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4.  CONCLUSION 

Vorticity  is  produced  by  fluid  motion  relative  to  a  solid  surface,  tbe 
rate  of  production  being  greatest  in  regions  where  the  pressure  and  velocity 
in  the  primary  flow  change  rapidly,  such  as  at  corners  and  sharp  edges.  The 
kinetic  energy  of  the  motion  induced  by  this  vorticity  is  derived  from  the 
primary  flow,  and  vorticity  generation  accordingly  transfers  energy  from  that 
flow  to  generally  smaller  scale  vortex  motions.  In  this  paper  we  have 
examined  the  production  of  vorticity  at  the  sharp  edges  of  circular  apertures 
in  a  thin  elastic  plate  by  long  wavelength  sound  and  vibrational  motions  of 
the  plate  The  vorticity  is  diffused  from  the  edges  by  viscosity,  resulting 
in  the  dissipation  of  the  sound  and  structural  motions.  In  the  absence  of 
mean  flow  through  the  apertures  the  dissipation  is  caused  by  the  nonlinear 
convection  of  vorticity  from  the  apertures  and  by  viscous  damping,  both  of 
which  are  weak  because  the  growth  of  vorticity  tends  to  be  inhibited  by 
periodic  reversals  in  the  primary  motions,  which  produce  vorticity  of 
fluctuating  sign.  The  damping  is  significantly  increased,  however,  when  a 
high  Reynolds  number  mean  flow  is  maintained  through  the  apertures.  Viscous 
effects  are  now  important  only  in  the  immediate  vicinity  of  the  aperture 
edges,  where  vorticity  is  released  into  the  flow;  but  the  vorticity  is  swept 
away  from  the  plate  by  the  flow,  and  its  kinetic  energy  is  permanently  lost  to 
the  incident  acoustic  and  vibrational  motions. 

Our  analysis  of  this  mechanism  for  a  bias  flow,  perforated  elastic  plate 
has  confirmed  that  significant  attenuations  of  sound  and  flexural  vibrations 
are  possible  provided  a  sufficiently  large  fluctuating  pressure  gradient  can 
be  established  across  the  screen  to  cause  unsteady  production  of  vorticity  in 
the  apertures.  When  the  fluid  loading  is  large  (e.g.,  for  metallic  plates  in 
water)  an  unrestrained  elastic  screen  is  effectively  acoustically 
"transparent" ,  and  the  damping  of  sound  by  vorticity  production  tends  to  be 
negligible.  Significant  acoustic  attenuations  can  occur  for  a  lightly  loaded 
plate  (in  air,  say)  provided  the  frequency  w  satisfies  wh/c^  >  P^/p^- 

The  damping  of  bending  waves  (i.e.,  of  resonant  structural  waves  that 
propagate  subsonically  relative  to  the  fluid)  increases  with  the  fluid 
loading,  since  the  surface  pressure  fluctuations  produced  by  the  passage  of  a 
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bending  wave  of  given  surface  displacement  tends  to  increase  with  fluid 
density.  The  damping  is  greatest  for  aperture  Strouhal  numbers  S^.  =  wR/U^ 
between  about  0.03  and  3,  and  the  attenuation  experienced  by  bending  waves  on 
a  perforated  steel  plate  in  water  (of  the  order  of  4  -  6  dB  per  wavelength  of 
propagation)  can  then  exceed  or  be  comparable  with  that  possible  when  motions 
of  the  plate  are  heavily  damped  by  the  application  of  surface  coatings.  In 
air  the  predicted  attenuations  are  very  much  more  modest,  ranging  between  0.5 
-  1  dB  per  wavelength. 
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APPENDIX 


The  coefficients  A  ,  A  ,  Q 
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The  integrals  in  these  expressions  are  taken  along  the  positive  real  axis 
indented  to  pass  below  real  singularities  of  the  integrands  (equivalent  to 
replacing  w  by  w  +  ic  and  taking  the  limit  c  ■*  +0)  . 
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LIST  OF  PRINCIPAL  SYMBOLS 

source  coefficients  in  equation  (2,11)  (j  =1,  2) 
see  (2.17),  (2.22)  and  appendix 

bending  stiffness 
speed  of  sound 

uniform  plate  dispersion  function  (2.15) 
perforated  plate  dispersion  function  (3.8) 

Young’s  modulus 

see  equation  (2.6) 

plate  thickness 

Heaviside  unit  function 

(k^,0,kg),  wavenumber  in  plane  of  plate 

(ma;^/B)^^^,  vacuum  bending  wavenumber 

2R(r  -  iA) ,  Rayleigh  conductivity  of  aperture,  see  (2.28) 

p^h,  mass  of  plate  per  unit  area 

U  /c  ,  Mach  number 

number  of  apertures  per  unit  area 

pressure 

incident  bending  wave  pressure,  (2.25) 
scattered  pressure 

displacement  flux  through  the  aperture 
see  (2.17)  and  appendix 

aperture  radius 

cjR/U^  ,  Strouhal  number 

=*  |U,  vorticity  convection  velocity 

mean  bias  flow  velocity 

rectangular  coordinate  (j  =  1,  2,  3) 

(x^ ,0,X3 ) 


a  NttR^  ,  open  area  ratio 

7(k)  ^|^.|) 

e  fluid  loading  parameter,  (2.32) 
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f  displacement 

flexural  displacement  of  plate 
f  perturbation  displacement  in  the  aperture 

flexural  wave  amplitude 
C  scattered  flexural  disturbance 

^  s 

AC  flexural  wavenumber 

AC  (j)/c  ,  acoustic  wavenumber 

0  'o’ 

power  dissipated  at  the  aperture,  (2.27) 
IIj  incident  bending  wave  power,  (2.30) 

fluid  density 

p^  density  of  plate  material 

a  Poisson’s  ratio 

cj  radian  frequency 

H  reflection  coefficient,  (3.4) 

3  transmission  coefficient,  (3.4) 
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Figure  1.  Configuration  of  the  plate,  aperture  and  bias  flow 
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Figure  4.  Reflection  and  transmission  of  sound  by  a  bias -flow  perforated 
plate. 
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CHAPTER  2 


ENERGY  CONSERVATION  AND  THE  DAMPING  OF  FLEXURAL  WAVES 

BY  VORTICITY  PRODUCTION 
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SUMMARY 

An  analysis  is  made  of  the  transfer  of  energy  between  structural/acoustic 
vibrations  and  fluid  kinetic  energy  in  flows  at  very  small  Mach  number.  A 
general  energy  balance  equation  is  discussed  for  a  vibrating  rigid  body  in 
either  incompressible  or  low  Mach  number  compressible  mean  flows.  This 
equation  can  be  used  to  calculate  the  growth  or  decay  of  structural  or 
acoustic  oscillations,  and  to  locate  regions  of  the  flow  where  energy 
exchanges  are  significant.  A  similar  general  treatment  for  vibrating  elastic 
bodies  does  not  seem  to  be  possible,  except  in  simple  cases  involving  thin 
elastic  plates  in  parallel  flow.  Two  such  problems  are  discussed,  involving 
the  dissipation  of  structural  vibrations  by  vorticity  production  (i)  at  the 
trailing  edge  of  a  large  elastic  plate,  and  (ii)  in  the  circular  apertures  of 
a  perforated  elastic  plate  in  a  two-sided  grazing  mean  flow.  The  interaction 
of  boundary  layer  turbulence  with  the  apertures  of  a  perforated  plate  can  be  a 
particularly  intense  source  of  sound,  but  in  applications  where  the 
characteristic  frequencies  are  small,  a  grazing  flow  perforated  screen  has 
been  shown  to  be  an  efficient  sink  of  acoustic  energy.  In  this  paper 
predictions  are  given  for  the  damping  of  bending  waves  by  the  same  mechanism. 
When  the  fluid  loading  is  large,  such  as  for  a  steel  plate  in  water,  these 
predictions  indicate  that  the  damping  of  resonant  bending  waves  can  exceed 
that  normally  achieved  by  coating  the  plate  with  elastomeric  damping 
materials,  at  least  over  a  restricted  range  of  frequencies. 
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1 .  INTRODUCTION 

Vorticity  is  produced  by  fluid  motion  relative  to  a  solid  surface,  the  rate 
of  production  being  greatest  in  regions  where  the  pressure  and  velocity  in  the 
primary  flow  change  rapidly,  such  as  at  corners  and  sharp  edges.  The  kinetic 
energy  of  the  motion  induced  by  this  vorticity  is  derived  from  the  primary 
flow,  and  vorticity  generation  accordingly  transfers  energy  from  that  flow  to 
generally  smaller  scale  vortex  motions.  When  the  primary  motion  is  produced 
by  a  sound  wave  incident  on,  say,  a  sharp  edge,  vorticity  diffuses  from  the 
edge  by  viscous  action,  causing  the  sound  to  be  damped  (e.g.,  see  [1  -  3]). 
When  the  fluid  is  in  a  mean  state  of  rest  the  dissipation  is  caused  by  the 
nonlinear  convection  of  vorticity  from  the  edge  and  subsequent  viscous 
damping,  both  of  which  are  weak  because  the  growth  of  substantial  levels  of 
vorticity  is  inhibited  by  the  periodic  nature  of  the  acoustic  field,  which 
leads  to  the  generation  of  vorticity  of  fluctuating  sign.  The  damping  is 
significantly  increased,  however,  in  the  presence  of  a  high  Reynolds  number 
flow  [4  -  16] .  Viscous  effects  are  now  important  only  in  the  immediate 
vicinity  of  the  edge,  where  vorticity  is  released  into  the  flow;  the  shed 
vorticity  is  swept  away  from  the  edge  by  the  flow,  and  its  kinetic  energy  is 
permanently  lost  to  the  incident  acoustic  wave. 

Practical  devices  for  exploiting  this  mechanism  of  sound  control  usually 
involve  bias  or  grazing  flow  perforated  screens.  In  the  bias  flow  case  a  mean 
pressure  difference  is  maintained  across  the  screen  producing  a  steady  flow 
through  the  apertures.  Damping  is  caused  by  the  modulation  of  vorticity 
production  at  the  aperture  edges  by  the  impinging  sound  [6,  7,  12,  14  -  16]. 

A  grazing  flow  screen  works  in  a  similar  way:  unsteady  motion  produced  by  the 
sound  in  the  apertures  generates  vorticity  which  is  swept  away  by  the 
tangential  flow  past  the  screen.  In  practice  high  acoustic  intensities  are 
usually  accompanied  by  significant  structural  vibrations.  Since  the  near 
field  (non-acoustic)  pressure  fluctuations  produced  by  a  vibrating  body  can 
modulate  vorticity  production  at  an  edge,  it  seems  likely  that  a  neighboring 
perforated  structure  might  also  be  used  to  suppress  structural  vibrations. 

This  possibility  has  been  examined  theoretically  by  the  author,  who  considered 
the  damping  of  bending  waves  (i)  incident  on  the  trailing  edge  of  a  large 
plate  in  a  mean  flow  [17],  and  (ii)  propagating  on  a  bias-flow  perforated 
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elastic  screen  [18].  For  a  bias  flow  screen  subject  heavy  fluid  loading  (a 

steel  plate  in  water,  say),  it  was  predicted  that  damping  is  significant  over 

a  range  of  frequencies  centered  on  a  Strouhal  number  wR/U  -  0(1) ,  where  co 
denotes  radian  frequency,  R  the  aperture  radius,  and  U  the  mean  flow  velocity 

through  an  aperture.  These  frequencies  are  typically  less  than  about  O.lco^, 

where  a?  is  the  coincidence  frequency  of  the  plate  [19,  20].  The  predicted 
damping  can  be  comparable  with  that  normally  achieved  by  heavily  coating  the 
vibrating  plate  with  a  conventional  elastomeric  material,  whose  use  frequently 
entails  substantial  and  undesirable  increases  in  the  structural  mass, 
especially  at  low  frequencies. 

In  this  paper  equations  are  derived  in  §2  which  describe  the  transfer  of 
energy  from  acoustic  and/or  structural  vibrations  to  the  kinetic  energy  of  the 
fluid  in  low  Mach  number  flows,  first  for  vibrations  of  arbitrary  amplitude  of 
rigid  and  elastic  bodies,  and  second  for  small  amplitude  motions  of  elastic 
plates  in  parallel  mean  flow.  These  equations  generalize  results  given 
previously  in  [12].  Illustrative  applications  of  the  theory  are  made  in  §3, 
and  in  §4  the  analysis  of  [18]  for  a  bias  flow  perforated  plate  is  extended  to 
investigate  the  damping  of  bending  waves  on  a  grazing  flow  perforated  screen. 
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2.  THE  ENERGY  EQUATION  IN  FLOW  AT  INFINITESIMAL  MACH  NUMBER 
2.1  Rigid  body  vibrating  in  incompressible  flow 

Equations  describing  the  damping  of  structural  vibrations  by  vorticity 
production  are  easily  established  in  the  simplest  possible  case  of  a  rigid 
body  in  an  incompressible  fluid.  Consider  a  rigid  body  with  surface  S 
immersed  in  incompressible  fluid  of  infinite  extent  and  uniform  mean  density 
p  and  shear  viscosity  .  Let  the  body  execute  translational  oscillations  at 
velocity  u°  (whose  mean  value  is  assumed  to  vanish)  with  respect  to  a  fixed 
coordinate  system  (x^jX^jX^),  and  let  V  denote  the  fluid  region  bounded 
internally  by  S  and  externally  by  a  large  control  surface  S  (Figure  1).  The 
fluid  at  large  distances  from  the  body  is  in  steady  motion  at  uniform  velocity 
U^ .  The  velocity  v  at  any  point  in  the  fluid  can  then  be  cast  in  the  form 

V  =  U  +  u  +  V  (1) 

CJ 

where  U  and  u  are  the  ir rotational  velocities  due  to  the  steady  mean  flow  and 
the  translational  oscillations  of  the  body  that  would  exist  if  the  fluid  were 
ideal,  and  v^  is  the  additional  velocity  required  to  satisfy  the  no-slip 
condition  on  S  and  to  account  for  the  velocity  induced  by  vorticity  in  the 
flow.  This  vorticity  includes  "gusts"  convected  in  the  incident  stream,  and 
that  produced  by  vortex  shedding  from  the  body.  The  velocities  can  be  defined 
more  precisely  as  follows.  Let  denote  the  velocity  potential  of  the  flow 

that  would  be  produced  in  an  ideal  fluid  by  translational  motion  of  the  body 
at  unit  speed  in  the  i- direction.  Then 


U  =  U°VX  u=u°V<^;,  ] 

j  j  j  j  [■  (2) 

X.  =  X.  -  <t>*  ,  J 

J  J  J 

where  U»n  =  0,  u*n  =  u°«n  on  S,  and  n  is  the  unit  normal  on  S  directed  into 
the  fluid  (see  Figure  1) .  The  remaining  velocity  v^  is  determined  in  terms  of 
the  vorticity  w  by 


U  +  u  +  v  =u°,  v«n=0onS, 

U  *  (J 

curlv  =  0?  in  V. 

U) 


(3) 
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At  large  distances  from  S,  1 |  -  0(l/|x|^),  and  this  and  conditions  (3)  are 
sufficient  to  determine  v  [21]. 

CJ 


The  kinetic  energy  of  the  fluid  (U+u+v^)^d^x  can  be  partitioned 

(with  the  assistance  of  the  divergence  theorem)  in  the  form 


Iv 


V  )^d^x 

(J  ' 


Js+s 


V  )*dS 

W  ' 


I. 


\p  u^d^x. 

2 


The  surface  integral  may  be  discarded.  The  integration  over  S  vanishes 
identically  because  (U+v^)*n  =  0  on  S .  On  S,  p^u°^j  ~ 

-p  (x  P  /|x|^)U^,  where  P.  =  X.u°»dS  is  the  impulse  of  the  oscillatory 
motion  [21].  The  integration  over  S  accordingly  vanishes  when  S  is  large  and 
geometrically  symmetric;  more  generally  its  average  value  is  null,  since  u° 
has  zero  mean.  The  final  volume  integral  is  the  ("reversible")  kinetic  energy 
of  the  local  irrotational  motion  associated  with  the  added  mass  of  the  body; 
it  is  equal  to  u°u° ,  where  A.^  is  the  added  mass  tensor  for  translational 

motion  [21] . 


The  transfer  of  energy  from  the  body  to  the  fluid  is  balanced  by  an 
increase  in  the  kinetic  energy  of  the  fluid  and  the  dissipation  by  viscosity, 
and  can  be  evaluated  from  the  rate  of  working  of  the  net  surface  force  F.  Let 
V  denote  fluid  velocity  relative  to  the  body,  i.e., 

rei 

V  =  V  -  u°  =  U  -  u°VX.  +  v  .  (4) 

rel  J  J  ^ 

The  net  force  F  exerted  on  the  fluid  can  be  written  [22,  23] 


F. 


at 


’o  Jv 


WAV  •VX.  d®x  + 

rel  1 


wAVX.  "dS, 
1 


(5) 


where  the  volume  integral  in  (5)  corresponds  to  the  component  of  force 
produced  by  the  normal  (pressure)  stresses  on  S,  and  the  surface  integral  is 
the  net  contribution  from  the  skin  friction.  The  rate  at  which  energy  is 
supplied  to  the  fluid  is  equal  to  F.u°. 


If 


(6) 
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denotes  the  viscous  energy  dissipation  rate,  Q  plus  the  rate  of  change  of  the 
fluid  kinetic  energy  can  be  equated  to  to  yield  the  energy  equation 


WAV  •VX.  d  X  +  r]u. 

rei  j  J 


^  wAVX.*dS, 
s  J 


where  E  =  (U+v^ )^d^x  is  the  kinetic  energy  of  the  mean  flow  and  the 

vortex  field.  The  contributions  from  the  added  mass  have  canceled  from  both 
sides  of  this  equation,  which  therefore  expresses  the  rate  of  dissipation  and 
production  of  non-reversible  kinetic  energy  to  the  rate  of  working  of  the 
surface  stresses  on  S.  For  a  stationary  body  (u°  =  0)  the  kinetic  energy  is 
absorbed  entirely  by  viscous  action  in  the  body  of  the  fluid.  There  are 
frequently  circumstances,  however,  in  which  vortex  shedding  from  a  vibrating 
body  becomes  highly  organized  and  correlated  over  large  distances  on  S. 

Kinetic  energy  of  the  flow  can  then  be  extracted  to  drive  unstable 
oscillations  of  the  body  [24,  25]. 

Let  us  repeat  the  above  analysis  in  a  reference  frame  in  which  the  fluid 
is  at  rest  at  infinity.  The  kinetic  energy  is  now 

|v  +  J, 

The  second  integral  is  the  reversible,  added  mass  component  of  the  energy, 

which  is  equal  to  |a.  .  (u° -U° ) (u° -U° ) .  The  force  between  the  fluid  and  body  is 

still  given  by  equation  (5),  but  the  velocity  of  translation  of  the  body  is 

equal  to  -l)®+u°  ,  so  that  the  rate  of  working  of  the  surface  force  is 

F  (u°-U°).  Evaluation  of  this  using  (5),  leads  to  the  following  alternative 
j  j  j 

energy  equation 


+  Q  =  “AVX..dS 

+  wAv^^^'Ud^x  -  r?U°Jg  wAVXj'dS, 


where  E  =fipv^d^xis  the  kinetic  energy  of  the  vortex  field  alone 
(hereafter  referred  to  as  the  "vortex  kinetic  energy").  Equations  (7)  and  (8) 
correspond  to  alternative  views  of  the  motion.  If  the  body  does  not  vibrate 
(u°  =  0) ,  E  +  J^Qdt  is  constant,  and  mean  flow  kinetic  energy  is  transformed 
into  vortex  energy  and  thence  dissipated  by  viscosity. 
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By  subtracting  equations  (7)  and  (8) ,  the  variable  mean  flow  energy 
Eg  =  E  -  is  found  to  satisfy 

=  -  p  r  WAV  -Ud^x  +  r?U°r  wAVX.'dS,  (9a) 

^OJV  rel  jJS  J 

=  F  U°  (^^) 

i  i 

Since  the  component  of  -F  in  the  direction  of  the  uniform  mean  flow  If  is  the 
drag  experienced  by  the  body,  equation  (9)  equates  the  rate  of  change  of  the 
mean  flow  kinetic  energy  to  minus  the  rate  of  working  of  the  drag  force.  For 
a  body  at  rest  in  a  uniform  stream,  the  kinetic  energy  lost  by  the  mean  flow 
in  this  way  is  transformed  into  vortical  energy  and  dissipated  by  viscosity. 
For  a  vibrating  body  the  drag  force  can  be  negative]  the  body  experiences  a 
force  opposite  to  the  mean  flow  direction  and  proceeds  to  "swim"  against  the 
stream.  The  energy  transferred  to  the  fluid  goes  into  the  wake  and  into 
accelerating  the  mean  flow  [26] . 

2.2  An  arbitrary  elastic  body  in  low  Mach  number  compressible  flow 

In  a  compressible  fluid  there  can  be  a  transfer  of  energy  from  structural 
motions  into  both  sound  and  the  kinetic  energy  of  the  flow.  This  can  be 
important  if  sound  waves  are  incident  on  the  body,  or  when  it  is  necessary  to 
consider  the  production  of  aerodynamic  sound  by  surface  motions  and  vorticity. 
We  shall  consider  low  Mach  number  motions  and  write  the  velocity  in  the  form 

v=U+u+v  +  V<p,  where  n»u  =  n*u°  ,  n»Vcp  =  0  on  S .  (10) 

(j 

U  u  and  v  are  defined  as  before,  on  the  assumption  that  the  motion  is 

’  CJ 

incompressible,  and  the  velocity  potential  cp  is  taken  to  account  for  the  whole 
effect  of  compressibility.  For  an  elastic  body  the  surface  velocity  u°  will 
generally  vary  with  position  on  S .  In  a  nonuniform  mean  flow,  the  velocity 
field  U  can  be  assumed  to  be  incompressible  when  M^  1,  where  M  =  U°/c^  is 
the  free  stream  Mach  number,  c^  being  the  speed  of  sound  (which  can  also  be 
regarded  as  uniform  throughout  the  fluid  when  M^  <S  1). 
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An  equation  describing  the  production  of  the  kinetic  energy 
E  =  r  i/>  (U+v  )^d^x  of  the  mean  flow  and  vorticity  can  be  derived  from  the 
momentum  equation  taken  in  the  form 


(U+v  )  +  V  p  + 


\p  ^ 
2 


dcp'] 

p  —  +  p  —  I 

°dt  °dtJ 


p^wAv  -  fjcurlw  , 


(11) 


where  (for  <  1)  the  mean  density  p^  is  constant,  and  i)  denotes  the  velocity 
potential  of  incompressible  flow  produced  by  oscillatory  motion  of  the  body  in 
an  ideal  fluid  (so  that  u  =  VV>)  .  Dissipation  terms  proportional  to  f?divv  and 
those  associated  with  the  bulk  viscosity  of  the  fluid  have  been  neglected. 


The  scalar  product  of  equation  (11)  with  D+v^  can  be  written 

+  div[(U+vj[p  +  Ipy  +  +  pf^]]  =  p^a,Av.(VV.+Vcp) 

-  r/(U+v  )»curlw.  (12) 

(j 

The  argument  of  the  divergence  on  the  left  hand  side  may  be  interpreted  as  the 
flux  of  the  kinetic  energy  of  the  mean  flow  and  vorticity  when  convection  by 
sound  and  unsteady  surface  motions  is  neglected.  In  irrotational  regions  it 
reduces  exactly  to  the  energy  flux  vector  (U+v^  )d‘i>/5t ,  where  $  is  the 
velocity  potential  with  contributions  from  unsteady  surface  motion  and  sound 
excluded  [27] . 

The  rate  of  production  of  mean  flow  and  vortex  kinetic  energy  is  obtained 
by  integrating  equation  (12)  over  the  volume  V.  The  motion  of  the  surface  S 
of  the  body  is  accounted  for  by  recalling  that 


Jv  I'’. +  I 


ip  (U+V  )^u°«dS 

2  '  u  ' 


If  Q  is  the  viscous  dissipation  calculated  as  for  an  incompressible  fluid  (as 
in  (6)  but  with  v  defined  by  (10)),  equation  (12)  then  yields 


aE 

at 


+  Q  =  p^l^  WA(U+V^)*(V^  +  V(p)d^x  -  |p„|s  (U+v^)^u°«dS 

+  77 Jg  [«A(U+v^)  -  2(u°»V)vj»dS.  (13) 
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When  the  viscous  term  on  the  right  hand  side  of  this  formula  is  discarded,  and 
the  surface  is  at  rest  (u°  =  0) ,  this  equation  reduces  to  one  given  previously 
in  [7,  12].  However,  for  a  moving  surface  the  second,  surface  integral  on  the 
right  hand  side  of  (13)  was  mistakenly  omitted  from  the  corresponding  equation 
(2.5)  of  [12], 

For  a  rigid  body  u°  is  constant  on  S,  ^  =  u°^^  ,  and  the  identity 
a/ax.  (|u2)  =  a/ax.  (o.u^  -  (wau)^  ,  v  =  u+v^  (div  v  =  O) 

can  be  used  to  show  that 

ip  r  (U+v  )^u°*dS  =  p  u°  r  wA(U+v  )*Vx  d^x. 

Thus,  when  small  dissipative  terms  of  order  tjdiv  v  are  again  discarded,  so 
that  fg(u°«V)v*dS  =  0,  (13)  becomes 

^  +  Q  -  +  ,u;J, 

where 

V  r=v-u°=U  -  u°VX.  +  V  +  V<p.  (15) 

'^rel  3  3 

Equation  (14)  is  the  generalization  to  low  Mach  number  compressible  flow  of 
equation  (7). 

The  corresponding  equations  for  the  vortex  kinetic  energy  is  obtained 
by  taking  the  scalar  product  of  (11)  with  v^ .  This  gives 

+  Q  =  pj^  o,Av^.{U+VV.+V<p}d3x  -  ipjg  v^u°*dS 

+  ''Is 

For  a  rigid  body  this  simplifies  to 


-  2(u°»V)vj»dS. 


(16) 


wAVX. ^dS 

j 

-  r/Jg  ct>AV(p*dS,  (1^) 
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— ^  +  Q  = 
5t 


-p  u?  L  WAV  ‘VX.d^x  +  ri\x°  \  waVX  "dS 
'^o  jJV  reL  j  JJS  J 

+  p  U°  L  WAV  ‘VX.d^x  -  r?U°  r  wAVX  -dS 
jJV  rel  J  '  jJS  J 


+  P 


o  Jv 


wAv^^j^»V(pd 


>1. 


ti?AV(p*dS . 


(17) 


When  this  is  subtracted  from  (14),  the  fate  of  production  of  mean  stream 
kinetic  energy  =  E  -  is  found  to  be  given  by  (9a),  with  v^^^  defined  as 
in  (15). 

A  particular  form  of  the  rigid  body  formulae  (14)  and  (17)  arises  when  S 
is  acoustically  compact,  i.e.,  when  the  characteristic  wavelength  A  of  sound 
waves  involved  in  the  motion  is  much  larger  than  the  typical  body  dimension 
say.  The  generation  of  sound  by  the  unsteady  flow  near  the  body  can  then 
usually  be  neglected,  and  the  principal  contribution  to  the  compressible 
velocity  potential  (p  is  from  sound  waves  incident  from  remote  sources.  In  a 
first  approximation  (correct  to  the  neglect  of  terms  -  0{ (i/A)  })  the  behavior 
of  (p  near  S  is  given  by 

<p(x,t)  *  (p(x^,t)  +  Xja(p(x^ ,t)/axj , 

where  x  may  be  taken  as  the  center  of  volume  of  S.  Equations  (14),  (15)  and 
(17)  then  imply  that  the  production  of  kinetic  energy  proceeds  as  in  an 
incompressible  fluid  provided  the  vibrational  velocity  u°  of  the  body  is 
increased  by  -dcpix^  ,t) /dx. ,  i.e.,  by  an  amount  which  is  equal  and  opposite  to 
the  acoustic  particle  velocity  of  the  sound  incident  on  S. 


2.3  Small  amplitude  vibrations  of  elastic  plates  in  parallel  mean  flow 

The  surface  velocity  u°  cannot  normally  be  assumed  to  be  independent  of 
position  on  S  when  the  body  is  elastic,  and  it  is  not  then  possible  to  reduce 
the  energy  equation  (13)  to  the  simplified  form  (14).  An  important  special 
case,  however,  occurs  when  all  fluctuating  quantities  are  small  and  the 
viscous  term  on  the  right  of  (13)  may  be  discarded.  The  surface  integral 
involving  the  shear  viscosity  r)  is  negligible  when  the  Reynolds  number 
U°ip  /t]  >  1,  where  Jl  is  the  length  scale  of  the  unsteady  motions  near  S.  In 
making  this  approximation  observe  that,  although  viscosity  is  ultimately 
responsible  for  the  diffusion  of  vorticity  into  flow,  this  fact  has  no  bearing 
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on  the  relative  magnitudes  of  the  integrals  in  (13) . 

Consider  interactions  and  vorticity  production  in  the  neighborhood  of  one 
or  more  thin  elastic  plates  in  high  Reynolds  number  parallel  mean  flow  at 
velocity  .  When  all  fluctuating  quantities,  including  the  vorticity  w,  can 
be  regarded  as  small  perturbations,  the  leading  order  approximation  to  the 
terms  on  the  right  of  (13)  (omitting  contributions  of  zero  mean)  is 

a)AU°*(VV>  +  V<p)d^x  -  (U°»v^)VV-dS. 

These  integrals  can  be  combined  by  introducing  the  bound  vorticity  w,  say,  on 
the  surfaces  of  the  plates.  This  is  done  when  there  is  only  one  plate  as 
follows.  Let  the  equation  f(x,t)  =  0  define  a  surface  in  the  fluid  that  just 
encloses  the  plate,  where  $•  >  0  in  the  exterior  fluid.  Define  the  extended 
vorticity  Q  by 

Q  =  curl{H(f)v^}  =  H(r)w  +  VH(r)Av^ 

=  !!($■)“  + 

where  H(x)  is  the  Heaviside  step  function,  and  x^  is  a  local  coordinate  which 
vanishes  on  the  plate  and  is  orientated  in  the  direction  of  the  unit  normal  n. 
The  bound  vorticity  w  =  nAv  6(x  ).  Recalling  that  n»V(p  =  0  on  S ,  (13)  is  now 
reduced  to 

—  +  Q  =  p  r  QaU°»(V0  +  V(p)d^x,  (19) 

at  °J 

where  the  integration  is  taken  over  the  whole  of  space.  It  is  evident  that 
this  result  is  applicable  to  situations  involving  multiple  parallel  plates. 

Similarly,  since  If  is  constant  and  n*v^  =  0  on  S ,  the  vortex  and  mean 
Stream  energy  equations  become,  in  the  same  approximation, 

— ”  +  Q  =  p  L  WAV  ‘Ifd^x  =  -p  r  QAlf ‘v  d^x,  (20) 

'^ojv  u  '^oj  u 

— ®  =  p  r  QaU°*(VV'  +  V<p  +  V  )d^x.  (21) 

at  °J  " 

Equations  (19)  -  (21)  remain  valid  to  second  order  for  stationary  surfaces. 

In  that  case,  however,  Vip  =  0,  and  the  integrations  need  be  taken  only  over 
the  free  field  vorticity. 
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3.  APPLICATIONS  TO  RIGID  AND  ELASTIC  PLATES 
3 . 1  The  rigid  airfoil 

A  rigid,  two  dimensional,  flat  airfoil  of  chord  2a  is  set  at  zero  angle 
of  attack  to  an  incompressible,  high  Reynolds  number  mean  flow  of  velocity  U 
in  the  -direction.  In  the  undisturbed  state  the  airfoil  occupies  the  region 
I  <  a,  x^  =0,  of  the  rectangular  coordinate  system  (x^^.x^.x^).  The 
airfoil  executes  small  amplitude  time -harmonic ,  translational  oscillations  at 
velocity  u°  =  (0 , u^ , 0)e'^"^  in  the  x^-direction  (where  the  real  part  is  to  be 
taken,  and  the  harmonic  time  factor  is  henceforth  suppressed),  as  indicated  in 
Figure  2(a).  Viscous  forces  are  neglected  except  at  the  trailing  edge,  where 
vorticity  is  shed  into  a  wake  of  infinitesimal  thickness,  in  accordance  with 
the  Kutta  condition  of  thin  airfoil  theory.  Then  w  =  (0,0,0;^),  where  [28] 

=  7^H(x^-a)6(x2)e^'^^i,  7^  =  4iu2 /{H^^\/ca)+iH^^\/ca) }  ,  (22) 

are  Hankel  functions  [29],  and  k  =  c.;/U  is  the  wavenumber  of  the  wake 
vorticity.  In  high  Reynolds  number  flow  the  power  n  =  3E/3t  +  Q  fed  into  the 
mean  flow  and  wake  per  unit  span  of  the  airfoil  is  given  by  the  first  integral 
in  (7)  with  j  =  2  and  the  integration  over  x^  omitted.  To  evaluate  the 
integral  correct  to  second  order  in  the  amplitude  u^  of  the  oscillations,  we 
can  set  v  ,  =  (U,0,0)  and 

rei 

=  Re{ - iVz2 -a2 } ,  z  =  x^+ix^  (23) 

[21].  After  taking  the  real  parts  of  all  fluctuating  quantities  and  averaging 
over  the  period  In/w  of  the  motion,  we  find 

H  =  TrCj^p^aUlu^  1^  ,  (24) 

where  C  h  C  (/ca)  is  the  real  part  of  the  Theodorsen  function 

C(x)  =  iH<^\x)/{H'^>(x)+iH^^^(x)}.  (25) 

C  is  positive  for  all  real  frequencies,  so  that  the  motion  is  always  stable 
inasmuch  as  energy  is  always  drawn  from  the  oscillating  body  and  ceded  to  the 
fluid  [24] . 
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A  similar  calculation  shows  that  the  power  delivered  to  the  wake  per  unit 
span  is 


(/ca)  |H'^’(/ca)+iH[^\/ca)  |2 


(26) 


The  difference  Eg  =  II  -  represents  the  power  absorbed  directly  by  the  mean 
stream.  This  is  positive  and  arises  from  the  working  against  the  mean  stream 
of  the  suction  force  at  the  leading  edge  (which  is  not  balanced  by  an  equal 
and  opposite  force  at  the  trailing  edge  because  of  the  application  of  the 
Kutta  condition) .  The  airfoil  accordingly  has  a  tendency  to  "swim"  against 
the  mean  stream  under  the  action  of  a  propulsive  force  equal  to  Hg/U  (see  [26] 
and  the  references  cited  therein  for  an  extensive  discussion  of  problems  of 
this  type  relating  to  aquatic  propulsion,  etc) . 


When  the  characteristic  wavelength  of  the  vortex  wake  is  much  smaller 
than  the  airfoil  chord,  i.e.,  when  the  reduced  frequency  ks.  -»  it  follows  by 
use  of  the  asymptotic  formulae  for  the  Hankel  functions  in  (25)  and  (26)  that, 

n  =  n  =  jnp  aU|u  1^  as  /ca  (27) 

SO  that  the  power  absorbed  by  the  flow  is  ultimately  shared  equally  by  the 
mean  flow  and  wake.  At  low  frequencies,  however,  -  0(/ca)n,  and  the  energy 
flows  primarily  into  the  mean  stream. 


A  similar  calculation  can  be  performed  for  the  case  of  a  sound  wave 
incident  on  the  airfoil  in  mean  flow  of  Mach  number  M  =  U/c^  <  1,  when  the 
chord  is  acoustically  compact  (wavelength  of  the  sound  >  airfoil  chord) .  In 
these  circumstances  it  is  permissible  to  neglect  the  convection  of  sound  by 
the  flow,  and  we  may  consider  an  incident,  time»harmonic  plane  acoustic  wave 
whose  pressure  is  given  by 


p  =  +  x^sin^)^  (28) 

where  p  is  constant,  k  =  la/c  is  the  acoustic  wavenumber  and  6  is  the 

”  I  O  O 

propagation  direction  measured  from  the  positive  x^^-axis  (see  Figure  2(b)). 
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Assume  the  airfoil  to  be  fixed.  For  a  compact  airfoil  K^a.  is  small,  and 
the  velocity  potential  of  the  acoustic  field  near  the  airfoil  can  be 
approximated  by 


IP  r 

cp  =  -  — 1  +  i/C  X  cos^  +i/c  X  s 

p  wl  °  ^  °  ^ 


in^j  , 


where  is  given  by  (23).  The  wake  vorticity  is  again  given  by  equations 
(22)  with  u^  replaced  by  -p^ sin^/p^c^ .  Using  these  results  in  (14)  and  (17) 
(in  which  only  the  volume  integrals  involving  <p  are  retained),  we  find, 
correct  to  second  order, 


n 

-  =  7rC_  (/ca)Msin^  , 

TT 


n„  2Msin^ 

w  _ 

(/ca)|H^^\/ca)+iH[“(/ca)|2 


(30) 


where  =  | p^ | ^ asin^/p^ c^  is  the  acoustic  power  incident  on  unit  span  of  the 
plate.  At  high  reduced  frequencies  the  absorbed  sound  power  is  shared  equally 
between  the  wake  and  mean  stream,  but  as  /ca  0,  the  wake  power  becomes 
negligible  relative  to  that  ceded  to  the  mean  stream.  As  in  the  case  of  the 
vibrating  plate,  the  sound  induces  a  leading  edge  suction  force  and  a  tendency 
for  the  plate  to  be  propelled  against  the  stream. 


The  interaction  of  sound  with  the  trailing  edge  of  a  stationary  rigid 
airfoil  whose  chord  is  much  larger  than  the  acoustic  wavelength  (/c^a  >  1)  can 
also  be  investigated  in  this  way  [12,  30].  When  the  airfoil  is  modeled  by  a 
semi- inf inite  rigid  plate  in  flow  of  infinitesimal  mean  flow  Mach  number,  and 
the  incident  wave  is  defined  as  in  (28) ,  the  acoustic  power  dissipated  per 
unit  span  of  the  trailing  edge  is  [12] 


n - cos  (-)  , 

p  U)  ^ 


M  «  1, 


In  this  limit  it  is  again  found  that  ,  i.e.,  that  the  dissipated 

acoustic  power  is  shared  equally  by  the  wake  and  mean  flow,  but  this  relation 
is  not  maintained  at  higher  Mach  numbers  (see  [12]  for  details). 
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3.2  Damping  of  bending  waves  on  a  large  plate 

Vortex  shedding  from  the  trailing  edge  of  a  thin  elastic  plate  is 
probably  the  simplest  model  problem  for  studying  the  exchanges  of  energy  which 
occur  when  a  vibrating  elastic  structure  interacts  with  mean  flow.  Small 
amplitude  flexural  motions  of  the  plate  are  governed  by  the  linearized  bending 
wave  equation  [20],  and  the  canonical  structural  vibration  problem  is  that  in 
which  a  bending  wave  impinges  normally  on  the  edge.  Part  of  the  incident  wave 
energy  is  reflected,  part  is  transformed  into  kinetic  energy  of  the  mean  flow 
and  wake  and,  for  a  compressible  fluid,  part  is  scattered  into  sound.  A 
linearized  analysis  of  these  events  is  discussed  in  [17]  for  a  plate  that  is 
immersed  in  a  mean  flow  of  infinitesimal  Mach  number  flow. 

Let  a  thin,  semi-infinite  elastic  plate  of  bending  stiffness  B  and  mass  m 
per  unit  area  occupy  the  half -plane  x^  <0,  x^  =  0,.  in  fluid  in  uniform  motion 
at  speed  U  (where  M  =  U/c^  4  1)  in  the  positive  x^ -direction  (Figure  3). 
Time-harmonic  flexural  displacements  f(x^)e‘^“''  (w  >  0)  of  the  plate  (measured 
in  the  X2-direction)  satisfy  the  bending  wave  equation  [20], 

{BaV^xJ  -  maj^)f  =  p(x^,-0)  -  p(x^,+0),  Xj^  <  0,  (32) 


where  p(x^,X2)  denotes  pressure,  and  the  terms  on  the  right  hand  side  give  the 
net  normal  stress  on  the  plate  (here  and  henceforth  the  harmonic  time  factor 
is  suppressed)  .  The  displacement  and  the  pressure  are  also  related  by 
the  x^ -component  of  the  fluid  momentum  equation  {w+iUa/ax^ f  =  dp/dx^,  x^  = 
±0. 


A  steady  state  bending  wave  C  =  propagating  towards  the  edge  has 

subsonic  phase  velocity  (<  c^)  and  the  structural  motion  is  accompanied 

by  a  traveling  pressure  field  that  decays  exponentially  with  distance  from  the 

plate 


p  =  -sgn(x  )p  {(w-UK^)^/VK2-/c2)C^e^'‘+^l‘^''+'''o  '^2  I  . 


(33) 


provided  the  Mach  number  is  sufficiently  small  that  convection  of  sound  by  the 
flow  can  be  neglected 


It  follows  from  (32)  that  is  the  positive  real  root  of  the  dispersion 
equation 
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D(k,w)  =  Bk''  (w-Uk)^/Vk2 -/c2  =  0.  (34) 

A  reflected  bending  wave  f  emerges  from  the  localized, 

complex  fluid- structure  interaction  at  the  edge  and  propagates  without 
attenuation  to  x^  =  -®.  The  reflection  coefficient  R  is  a  function  of 
frequency,  Mach  number  M,  and  the  mechanical  constraints  imposed  on  the  motion 
of  the  plate  at  the  edge;  /c_  is  the  positive  root  of  (34)  when  the  sign  of  U 
is  reversed.  In  addition,  sound  is  generated,  and  vorticity  is  shed  into  a 
wake  extending  downstream  from  the  edge.  The  strength  of  the  shed  vorticity 
can  be  estimated  by  application  of  the  unsteady  Kutta  condition.  The  analysis 
of  this  problem  in  [17]  is  based  on  linear  perturbation  theory,  in  which  the 
wake  is  treated  as  a  vortex  sheet. 

The  flexural  wave  power  11^  ,  say,  incident  on  unit  length  of  the  edge 
(including  that  conveyed  in  the  evanescent  fluid  motions  on  both  sides  of  the 
plate)  can  be  calculated  from  the  formula  11^  =  |  1  ^aD/3k,  where  the 

derivative  is  evaluated  at  k  =  [31],  The  total  flexural  wave  power 

dissipated  per  unit  length  of  edge  consists  of  a  portion  scattered  into 
sound  and  a  component  Hj,  which  increases  the  kinetic  energy  of  the  mean  flow 
and  wake.  The  total  dissipated  power  is  therefore  given  by 


(35) 


Using  formulae  given  in  [17],  the  fractional  power  dissipated  at  the 
edge,  A  =  {II  +11  } /H  ,  has  been  plotted  in  Figure  4(a)  for  a  steel  plate  in 
water  for  several  values  of  the  Mach  number  when  the  edge  of  the  plate  is 
"free"  (S^f/Sx^  =  =  0  at  x^  =  -0  [20] ) .  In  this  figure  the  frequency 

is  normalized  by  the  coincidence  frequency  o>^  =  c^  (m/B)  ^  ,  above  which  the  in 
vacuo  bending  wave  phase  speed  exceeds  c^ . 


For  M  =  0  all  of  the  dissipation  of  bending  wave  power  can  be  attributed 
to  sound  radiated  from  the  edge.  When  M  0  there  is  a  significant  increase 
in  the  bending  wave  power  loss  at  low  frequencies  caused  by  the  transfer  of 
energy  to  the  mean  flow  and  wake,  with  negligible  changes  in  the  radiated 
sound  power  (at  least  for  M  «  1) .  This  is  clear  from  Figure  4(b),  which 
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compares  the  total  power  loss  A  for  a  steel  plate  in  water  (solid  curve)  with 
the  acoustic  power  (dashed  curve)  when  M  =  0.01.  Inspection  of  Figures  4(a) 
and  (b)  confirms  that  the  sound  power  is  hardly  influenced  by  the  flow. 
However,  at  high  frequencies  dissipation  is  caused  entirely  by  the  radiation 
of  sound;  this  is  because,  as  w  ^  ,  bending  wave  energy  is  contained 

principally  in  the  evanescent  motions  on  either  side  of  the  plate,  and 
propagates  with  negligible  plate  motion  at  a  velocity  which  is  only  slightly 
less  than  the  speed  of  sound. 

The  dotted  curve  in  Figure  4(b)  shows  the  fraction  of  the  dissipated 
energy  which  appears  as  vortex  kinetic  energy.  This  indicates  that,  over  the 
range  of  frequencies  within  which  acoustic  dissipation  is  negligible,  the 
dissipated  energy  is  shared  equally  between  the  mean  flow  (Eg)  and  the  vortex 
motions  (E^^);  this  is  not  generally  the  case,  but  depends  on  the  nature  of 
mechanical  constraints  at  the  edge.  The  momentum  of  the  mean  flow  increases 
at  a  rate  equal  to  Hg/U,  so  that  the  plate  experiences  an  equal  and  opposite 
thrust  and  a  tendency  to  swim  against  the  mean  stream. 
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4.  DAMPING  OF  FLEXURAL  WAVES  ON  A  GRAZING  FLOW  PERFORATED  SCREEN 


Perforated  screens  have  been  used  to  attenuate  sound  in  low  Mach  number 
flows.  The  flow  is  usually  turbulent,  and  typically  the  sound  is  generated  by 
the  interaction  of  the  mean  flow  and  turbulence  with  solid  bodies,  such  as  the 
cross -tubes  in  a  flow-through  heat  exchanger.  Ver  [10,  13]  has  described  the 
use  of  grazing  flow  perforated  plates  in  the  tube  bank  cavity  of  a  nuclear 
reactor  heat  exchanger  that  successfully  eliminated  acoustic  resonances  which 
had  not  responded  to  more  conventional  means  of  attenuation.  Acoustic 
attenuation  by  this  means  is  useful  only  when  the  dominant  frequencies  are 
small.  At  higher  frequencies  boundary  layer  turbulence  (which  is  enhanced  by 
flow  over  the  perforations)  can  interact  strongly  with  the  apertures  to 
produce  intense  sound  (and  vibration)  [32] .  Analytical  studies  [33  -  35] 
suggest  that  boundary  layer  generated  noise  is  increased  by  several  orders  of 
magnitude  by  the  presence  of  the  apertures. 

4.1  Modified  bending  wave  equation 


The  attenuation  of  low  frequency  sound  and  structural  vibrations  by  a 
grazing  flow  perforated  plate  can  be  estimated  by  a  simple  extension  of  the 
method  developed  in  [18]  for  a  bias  flow  screen.  A  generalized  form  of  the 
bending  wave  equation  (32)  was  derived  in  [18]  for  a  plate  with  circular 
apertures  when  the  length  scale  of  the  sound  and  plate  motion  is  large 
compared  to  both  the  aperture  radius  R  and  the  distance  between  neighboring 
apertures.  For  an  infinite,  thin  elastic  plate,  with  N  apertures  per  unit 
area,  and  occupying  the  plane  x^  =  0  in  the  undisturbed  state,  long 
wavelength,  small  amplitude  motions  (proportional  to  e  )  are  governed  by  the 
equation 


2aa 

(1-a)- 


BV^ 


mo) 


^  (  r  1  ff-i  2aa 

)c  +  (l  H.  2HRK,[l  -  — 


Bv: 


=  0. 


(36) 

where  a  =  NttR^  is  the  fractional  open  area,  a  is  Poisson’s  ratio  for  the 
plate,  =  (S^/ax^+a^/SXj)^ ,  and  is  the  Rayleigh  conductivity. 
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The  conductivity  determines  the  volume  flux  Q,  say,  through  an  aperture 
caused  by  a  long  wavelength  perturbation  [p]  =  p(x^ ,+0 )  -pCx^.-OjX^)  in  the 
pressure  difference  across  the  plate  according  to  the  relation 

Q  =  Kjj[p]/iPo“- 

This  definition  is  a  generalization  of  that  originally  proposed  by  Rayleigh 
[36,  37]  in  terms  of  the  electrical  analogue  Q  =  where  §  is  the 

perturbation  velocity  potential,  and  the  motion  is  imagined  to  be  driven  by 
the  potential  difference  [$]  across  the  aperture.  In  more  general  situations, 
however,  where  the  flow  in  an  aperture  can  contain  vorticity,  local  variations 
in  the  velocity  potential  need  not  be  continuous  through  the  aperture  and 
Rayleigh’s  definition  must  be  expressed  in  terms  of  the  pressure  (which  is 
continuous) . 

It  was  argued  in  [18]  that  is  well  approximated  by  its  value  for  a 
rigid  plate  provided  the  length  scale  of  the  plate  motion  is  much  larger  than 
the  aperture  radius  R.  has  the  dimensions  of  length  and  is  equal  to  2R  in 

a  stationary  ideal  fluid,  but  becomes  complex  valued  and  frequency  dependent 
in  the  presence  of  flow,  because  of  the  unsteady  shedding  of  vorticity  from 
the  aperture  edges . 

We  consider  the  application  of  equation  (36)  to  the  case  in  which  the  low 
Mach  number  mean  flow  is  tangential  to  the  plate  at  uniform  speed  U  in  the 
x^ -direction  on  both  sides  (Figure  5),  and  where  the  Reynolds  number  based  on 
aperture  radius  is  large.  In  the  linearized  approximation  the  pressure  p  and 
velocity  potential  $  are  then  related  by  p  =  Pg(iw  •  Vd/dx^)^ .  The 
conductivity  is  a  function  of  the  Strouhal  number  S  =  coR/U,  and  can  be  set  in 
the  form 

=  2R{r(S)  -  iA(S)}.  (38) 

The  real  and  imaginary  parts  F  and  A  are  not  known  in  analytic  form,  but 
have  been  computed  by  Scott  [38]  for  high  Reynolds  number  flow.  Scott  applied 
linear  perturbation  theory  and  the  unsteady  Kutta  condition  to  determine  the 
strength  of  vorticity  shed  from  the  upstream  semi-circular  arc  of  the  aperture 
as  a  result  of  the  unsteady  aperture  flow;  the  vorticity  was  assumed  to  occupy 
a  vortex  sheet  in  the  plane  of  the  aperture  and  to  convect  downstream  at  the 
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mean  velocity  U. 

The  calculated  dependencies  of  T  and  A  on  Strouhal  number  are  depicted  in 

Figure  6.  F  and  A  are  approximately  periodic  (with  period  =  tt)  when  wR/U 

exceeds  about  3.  The  imaginary  component  A  governs  the  damping  of  aperture 

motions  by  vorticity  production  and  acoustic  radiation  losses .  The  latter  can 

be  neglected  when  k^R  is  small,  so  that  the  value  of  A  then  determines  the 

exchange  of  perturbation  energy  with  the  mean  flow;  aperture  motions  are 

damped  by  vorticity  production  provided  A  >  0;  when  A  <  0,  the  unsteady  motion 

in  the  aperture  is  enhanced  at  the  expense  of  mean  stream  kinetic  energy. 

'j']^i3  ©xtraction  of  energy  from  the  flow  is  caused  by  the  interaction  of  shed 

vorticity  with  the  semi-circular  "leading  edge"  on  which  it  impinges  after 

convecting  across  the  aperture.  The  aperture  motion  should  not  be  regarded  as 

unstable,  however,  since  the  transfer  of  energy  from  the  mean  flow  occurs  only 

when  an  incident  disturbance  contains  a  component  with  frequency  within  an 

interval  where  A  <0.  This  means  that,  in  the  absence  of  an  incident 
R 

pressure  differential  [p]  (and  on  the  basis  linear  theory),  large  amplitude, 
unsteady  motions  within  the  aperture  cannot  arise  spontaneously  as  a  result  of 
a  random  fluctuation  in  the  mean  flow;  mathematically  this  conclusion  follows 
from  the  fact  that  the  Rayleigh  conductivity  for  this  problem  can  be 

shown  to  be  a  regular  function  of  u)  in  the  upper  half  of  the  complex  frequency 
plane . 


4.2  Attenuation  of  resonant  bending  waves 


To  examine  the  damping  of  bending  waves  consider  a  time  harmonic 
displacement  f  >  q)  ;  from  (33)  we  have 

[p]  =  -2p^  {  (w-U/c)VV/c2-K2}f^e^''''i. 

The  flexural  wave  dispersion  equation  is  obtained  by  substituting  these 
expressions  for  C  and  [p]  into  equation  (36),  which  supplies 


D  (/c  ,  to)  = 


2ac7 


aa 

-cr). 


(1-C7) 


B/c^  -  mtj^ 


2p-  (co-U/c) 


^  +  2NRK„[l  -  -  "-")]]  -  0.  (39) 


When  there  is  no  mean  flow  -*  2R,  and  D(/c,w)  has  two  real  zeros  of 
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opposite  sign  that  correspond  to  undamped  flexural  waves  propagating  in  the 
positive  and  negative  -directions .  They  satisfy  |/c|  >  ,  so  that  V/c^-zc^  is 

real  and  positive.  In  the  presence  of  flow  the  bending  wavenumber  acquires  a 
small  imaginary  component  that  accounts  for  the  exchange  of  wave  energy  with 
the  flow  via  vorticity  production.  The  wave  power  dissipated  in  propagating  a 
distance  6s  is  equal  to  20 1  Im{/c)  |  dsiog^^  (e)  dB,  and  the  power  absorbed  per 
wavelength  of  propagation  is  therefore  (e)Im{K;}/Re{K}  » 

54.6Im{K)/Re{K;)  dB.  Figures  7  -  9  illustrate  typical  predictions  of  the 
absorption  for  a  perforated  steel  plate  in  water  at  different  fractional  open 
areas  a  and  aperture  radii  (for  R/h  =  3  and  6,  where  h  is  the  thickness  of  the 
plate)  and  mean  flow  velocities.  The  solid  and  dashed  curves  are  respectively 
for  waves  propagating  in  the  ±x^ -direction,  i.e.,  with  and  against  the  mean 
stream.  The  case  of  a  lightly  fluid  loaded  plate  is  illustrated  in  Figure  10 
(aluminium  in  air:  R/h  =10,  U  =  10  m/s).  The  ability  of  the  bending  wave 
motions  to  force  fluid  through  the  apertures  is  now  much  reduced,  and  the 
maximum  losses  are  correspondingly  very  much  smaller  than  for  a  plate  in 
water.  Waves  propagating  against  the  mean  flow  are  more  strongly  damped  than 
those  propagating  in  the  mean  flow  direction.  Increases  in  flow  velocity  or 
aperture  radius  reduce  the  maximum  possible  attenuation.  In  all  case,  waves 
can  experience  small  negative  damping  when  coR/U  >  3. 

To  assess  the  significance  of  these  predictions ,  it  should  be  noted  that 
the  largest  attenuations  usually  achieved  by  conventional  means,  i.e.,  by 
coating  a  nominally  homogeneous  plate  with  one  or  more  layers  of  an  absorptive 
material,  are  of  the  order  of  4  dB  per  wavelength  of  propagation  [19].  The 
present  results  are  accordingly  comparable  with  the  best  obtainable  by 
coating,  at  least  for  Strouhal  numbers  in  the  interval  0.7  -  2.  However, 
conventional  coating  materials  tend  to  be  effective  over  a  much  broader  range 
of  frequencies. 
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5.  CONCLUSION 

Energy  is  transferred  between  structural  and  acoustic  vibrations  and  the 
kinetic  energy  of  essentially  incompressible  fluid  motions  by  vorticity 
production  at  the  surface  of  a  solid.  In  many  instances  this  can  lead  to 
unstable  and  catastrophic  structural  motions  and  the  generation  of  intense 
acoustic  noise,  but  there  exist  regimes  in  which  vibrational  and/or  acoustic 
modes  may  be  dissipated  through  the  generation  of  vorticity.  This  energy 
transfer  mechanism  has  been  discussed  for  elementary  rigid  body  vibrations  of 
a  solid  in  low  Mach  number  flow,  and  the  effects  of  structural  and  acoustic 
damping  have  been  illustrated  by  reference  to  the  classical  problem  of  the 
vibrating  airfoil.  At  low  Mach  numbers  the  dissipated  energy  appears  first  as 
an  increase  in  the  kinetic  energy  of  the  mean  flow  and  wake  before  being 
ultimately  dissipated  by  viscosity.  Modified  versions  of  these  formulae  are 
applicable  also  to  vibration  problems  involving  parallel  flat  plates  in  a  mean 
flow,  but  such  problems  are  usually  so  complicated  that  it  is  generally  easier 
to  make  direct  predictions  of  the  energy  distribution  from  the  analytical 
solution.  The  simplest  problem  of  this  kind  involves  the  generation  of 
vorticity  at  the  trailing  edge  of  a  large,  vibrating  elastic  plate  in  a 
nominally  uniform  mean  flow.  At  low  Mach  numbers  energy  is  dissipated  both  by 
the  generation  of  sound  and  by  the  conversion  of  plate  energy  into  kinetic 
energy  of  the  flow.  At  low  frequencies  (well  below  the  coincidence  frequency) 
most  of  the  dissipated  energy  is  ceded  to  the  flow,  but  as  the  frequency 
increases  towards  coincidence  the  damping  is  progressively  dominated  by 
acoustic  radiation. 

The  exploitation  of  vorticity  production  in  the  apertures  of  a  perforated 
elastic  plate  for  vibration  damping  has  been  examined  for  low  Mach  number, 
two-sided  grazing  mean  flow.  As  in  the  case  of  the  ’‘bias  flow"  screen 
discussed  in  [18],  the  efficiency  with  which  bending  waves  are  damped 
increases  with  the  fluid  loading.  The  attenuation  experienced  by  waves  on  a 
p03^forated  steel  plate  in  water  is  significant  for  aperture  Strouhal  numbers  S 
=  coR/U  in  the  range  0.7  -  2,  and  can  exceed  that  usually  obtained  when  an 
unperforated  plate  is  heavily  coated  with  damping  material. 
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Figure  1. 


A  rigid  body  in  oscillatory  motion  in  an  incompressible  fluid 
in  uniform  mean  flow  at  velocity  T3°. 
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^2 


Figure  2.  (a)  Rigid  airfoil  vibrating  in  a  mean  stream. 

(b)  Sound  incident  on  a  stationary  airfoil  in  a  mean  stream. 
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Figure  3.  Flexural  vibrations  near  the  trailing  edge  of  a  large  plate  in 
uniform  mean  flow. 
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Figure  5.  Schematic  grazing  flow  perforated  elastic  plate 
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Figure  7 


S  =  coR/U 


Flexural  wave  power  dissipated  per  wavelength  of  propagation  on  a 
perforated  steel  plate  in  water  for  different  values  of  the 
fractional  open  area  a,  and  for  R/h  =  aperture  radius/plate 

thickness  =3,  U  =  5  m/s :  - ,  wave  propagating  in  the  mean 

stream  direction; - ,  wave  propagating  against  the  mean 


stream. 


Power  absorbed  per  wavelength  (dB) 
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Figure  10.  Flexural  wave  power  dissipated  per  wavelength  of  propagation  on  a 
perforated  aluminium  plate  in  air  for  different  values  of  the 

fractional  open  area  a,  and  for  R/h  =10,  U  =  10  m/s:  - ,  wave 

propagating  in  the  mean  stream  direction; - ,  wave  propagating 

against  the  mean  stream. 
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CHAPTER  3 


THE  INFLUENCE  OF  TANGENTIAL  MEAN  FLOW  ON  THE 
RAYLEIGH  CONDUCTIVITY  OF  AN  APERTURE 
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SUMMARY 

An  investigation  is  made  of  the  influence  of  grazing  mean  flow  on  the 
Rayleigh  conductivity  of  an  aperture  in  a  thin  rigid  plane.  The  fluid  is 
assumed  to  be  incompressible  and  inviscid,  but  the  Kutta  condition  is  applied 
to  permit  the  generation  of  vorticity  at  the  edge  of  the  aperture  by  an 
applied  time-varying  pressure  field.  Numerical  results  are  given  for  a 
circular  aperture  in  the  two  cases  of  (i)  one-sided  mean  flow,  when  the 
aperture  is  spanned  by  a  plane  vortex  sheet  in  the  undisturbed  state,  and  (ii) 
two-sided  mean  flow,  when  the  mean  velocity  is  the  same  on  both  sides  of  the 
plane,  so  that  the  undisturbed  motion  is  irrotational .  In  both  cases  there 
exist  frequency  ranges  within  which  perturbation  energy  is  either  absorbed  or 
generated  by  the  mean  flow.  The  numerical  results  are  supplemented  by  an 
approximate  analytical  treatment  of  the  same  problem  for  a  rectangular 
aperture  of  large  aspect  ratio  (with  its  long  edge  transverse  to  the  mean  flow 
direction).  The  aperture  flux  for  one-sided  flow  is  shown  to  be  absolutely 
unstable,  and  may  in  principle  be  triggered  by  an  arbitrary,  small  disturbance 
in  the  mean  stream.  For  two-sided  flow  the  motion  is  conditionally  unstable, 
in  the  sense  that  perturbations  are  amplified  by  the  extraction  of  energy  from 
the  mean  flow  only  when  the  frequency  of  the  applied  pressure  lies  in  certain 
discrete  bands. 
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1 .  INTRODUCTION 

Consider  the  unsteady  flow  of  fluid  through  an  aperture  of  area  A  in  a 
rigid  plane  of  negligible  thickness  produced  by  a  uniform,  time-harmonic 
pressure  load  [p]  =  (p+-p. (w  >  0)  applied  across  the  plane.  The 
"upper"  and  "lower"  faces  of  the  plane  are  taken  to  coincide  respectively  with 
surfaces  z  =  ±0  of  the  coordinate  plane  z  =  0  of  the  rectangular  system 
(x,y,z),  and  p^ ,  p_  are  the  respective  amplitudes  of  the  applied  pressures 
above  and  below  the  plane.  The  induced  motion  through  the  aperture  is  assumed 
to  be  incompressible  and  to  be  governed  by  the  linearized  equations  of  motion. 
Let  Q  denote  volume  flux  through  the  plane  (in  the  positive  z-direction) , 
where  here  and  henceforth  the  exponential  time  factor  e'^“^  is  suppressed. 

Then  the  Rayleigh  conductivity  of  the  aperture  is  defined  by  (Rayleigh 
1870,  1945) 

=  iwPoQ/(p+ -P.)  .  (1-1) 

where  denotes  the  density  of  the  fluid. 

The  pressure  varies  continuously  through  the  aperture,  and  p^  are  its 
limiting  values  at  large  distances  VA)  on  either  side.  When  there  is  no 
mean  flow,  and  the  motion  at  large  distances  can  be  regarded  as  irrotational , 
the  definition  (1.1)  is  equivalent  to 

\  ,  (1-2) 

where  are  the  uniform  velocity  potentials  above  and  below  the  plane, 
determined  by  p^  =  -p^dip^/8t.  This  form  of  the  definition  is  analogous  to 
Ohm* s  law  in  electricity.  Conductivity  has  the  dimensions  of  length;  for 
example,  it  is  equal  to  2R  for  a  circular  aperture  of  radius  R  in  an  ideal, 
stationary  fluid. 

For  an  ideal,  inviscid  fluid,  Kj^  is  real  and  the  motion  through  the 
aperture  is  entirely  irrotational  and  conservative.  In  practice,  however, 
viscous  action  causes  dissipation,  and  the  dissipated  power  11,  say,  is  given 
in  terms  of  [p]  and  the  normal  component  v^  of  fluid  velocity  at  z  =  0  by 

n  =  -]!.  <[p]v„)dxdy. 
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where  the  angle  brackets  denote  a  time  average,  and  the  real  parts  of  [p] , 
are  to  be  taken  after  restoration  of  the  time  factor  e  It  follows  that 

n  =  -|lm{Kj^}  I  [p]  IVPqC'^.  (1-3) 

and  therefore  that  energy  is  dissipated  provided  the  imaginary  part  of  is 
negative . 

The  dissipation  is  frequently  enhanced  for  certain  frequencies  by  the 
presence  of  a  high  Reynolds  number  mean  tangential  flow  over  the  aperture  or  a 
"bias"  flow  through  the  plane  (Bechert  1979;  Crighton  1981),  and  perforated 
screens  have  been  deployed  in  mean  flow  environments  to  eliminate  harmful 
acoustic  resonances  (for  example,  in  heat  exchanger  cavities  (Ver  1990)).  For 
a  bias  flow  screen  with  circular  apertures,  the  dissipation  can  be  estimated 
from  (1.3)  by  use  of  an  approximate  formula  for  given  by  Howe  (1979),  which 
has  been  validated  experimentally  by  Hughes  &  Dowling  (1990) .  The  mean  flow 
is  maintained  by  a  steady  pressure  differential  across  the  screen.  If  U  is 
the  average  velocity  within  an  aperture,  and  UVA/i/  ^  1,  the  motion  will  be 
uninfluenced  by  viscosity  except  in  the  immediate  vicinity  of  the  aperture 
edge,  where  the  flow  separates  to  form  a  jet.  The  undisturbed  motion  is  no 
longer  irrotational ,  and  a  superimposed,  fluctuating  pressure  [p]  causes 
additional  vorticity  to  be  generated,  which  is  then  swept  away  in  the  jet. 

The  energy  required  to  produce  this  vorticity  is  supplied  by  the  applied 
pressure  [p] ,  and  dissipation  occurs  by  its  transfer  to  the  kinetic  energy  of 
the  jet. 

In  this  paper  we  consider  the  determination  of  for  an  aperture  in  the 
presence  of  a  tangential  mean  flow.  When  the  mean  velocities  are  at  speeds 
,  U_  above  and  below  the  plane  (see  Figure  1)  and  the  Reynolds  number  is 
large,  the  undisturbed  mean  shear  layer  in  the  aperture  can  be  approximated  by 
a  vortex  sheet.  The  sheet  is  deformed  by  a  small  amplitude  applied  pressure 
[p] ,  and  additional  vorticity  is  shed  from  the  upstream  edge  of  the  aperture. 
However,  because  the  vortex  sheet  is  unstable  (Lamb  1932),  perturbation  energy 
may  be  expected  to  be  either  created  or  absorbed  by  the  interaction,  depending 
on  the  frequency  u)  of  the  applied  pressure.  This  problem  was  examined  by  Howe 
(1980),  who  obtained  high  and  low  frequency  approximations  for  for  a 
circular  aperture  when  ,  so  that  the  undisturbed  motion  may  be  regarded 
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as  irrotational .  The  applied  pressure  was  assumed  to  induce  separation  from 
the  upstream  semi-circular  edge  of  the  aperture  in  the  form  of  a  thin  vortex 
sheet,  the  strength  of  which  was  set  by  the  Kutta  condition  (Crighton  1985), 
and  it  was  concluded  that  Im{Kj^ )  <  0  for  all  to.  We  shall  argue  below  (§2) 
that  Howe’s  (1980)  solution  is  erroneous,  because  proper  account  was  not  taken 
of  conditions  at  the  downstream  edge  of  the  aperture. 

The  tangential  flow  problem  is  formulated  in  §2 ,  and  solved  numerically 
for  a  circular  aperture  for  the  two  cases  of  (i)  one-sided  mean  flow,  when  the 
fluid  is  at  rest  on  one  side  of  the  plane  in  the  undisturbed  state,  and  (ii) 
two-sided  flow  at  the  same  velocity  on  both  sides.  The  problem  is  amenable  to 
approximate  analytical  treatment  for  a  rectangular  aperture  of  large  aspect 
ratio,  orientated  with  its  long  edges  perpendicular  to  the  mean  flow 
direction.  This  is  discussed  in  §3,  and  comparison  made  with  the  numerical 
results.  In  all  cases  there  exist  frequency  ranges  where  Im{Kjj}  >  0,  within 
which  perturbation  energy  is  extracted  from  the  mean  flow.  For  both  the 
circular  and  rectangular  apertures  we  show  in  case  (ii),  for  which  the 
undisturbed  motion  is  irrotational,  that  the  motion  is  conditionally  unstable 
(§4),  i.e.,  that  energy  is  extracted  from  the  mean  flow  only  when  [p]  5^  0;  in 
case  (i),  on  the  other  hand,  the  motion  is  absolutely  unstable,  and  an 
arbitrarily  small  disturbance  can  trigger  the  release  of  mean  flow  kinetic 
energy . 

A  brief  discussion  is  given  in  the  Appendix  of  the  generalized  Rayleigh 
conductivity  Kjj(k,w),  which  is  defined  as  in  (1.1),  but  for  an  applied 
pressure  [p]  =  where  p^  is  constant  and  k  is  an  arbitrary 

"trace -wavenumber"  parallel  to  plane.  Rayleigh’s  (1870)  definition 
corresponds  to  Kj^(0,w).  The  reverse  flow  reciprocal  theorem  (Howe  1975; 
Mohring  1979),  extended  to  problems  involving  vortex  sheets,  can  be  invoked  to 
express  Kj^(k,a))  in  terms  of  the  vortex  sheet  displacement  in  the  reverse  flow 
problem  for  k  =  0. 
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2.  CIRCULAR  APERTURE  WITH  A  GRAZING  MEAN  FLOW 
2.1  The  vortex  sheet  model 

Consider  incompressible,  uniform  mean  flows  over  the  upper  and  lower 
surfaces  of  the  rigid  plane  z  =  0  which  is  pierced  by  a  circular  aperture  of 
radius  R  whose  center  is  at  the  coordinate  origin.  Let  the  mean  velocities  be 
at  speeds  in  the  x-direction  respectively  in  z  ^  0,  as  in  Figure  1.  The 
Reynolds  number  is  assumed  to  be  sufficiently  large  that  viscosity  can  be 
neglected  except  for  its  role  in  generating  vorticity  at  sharp  edges.  In  the 
steady  state,  when  ^  U_ ,  the  aperture  is  spanned  by  a  vortex  sheet. 

When  the  plane  is  subject  to  a  small  amplitude,  time  harmonic  pressure 
load  [p]  =  P^.-P.  >  the  velocity  potentials  (x)  of  the  unsteady  motion  above 
and  below  the  plane  may  be  expressed  in  the  form  (Lamb  1932) 


(x) 


^  [■“  v±(x’  ,y’) 
27r  Joo  |x-x’  I 


dx*  dy* 


X* 


(x’ ,y’ ,0) , 


(2.1) 


where  v^  is  the  z-component  of  velocity  as  z  ±0,  which  vanishes  on  the  plane 
except  at  the  aperture.  In  the  linearized  approximation,  the  z-component  of 
displacement  r(x,y)  of  fluid  particles  initially  lying  just  above  or  below  the 
plane  (z  =  ±0)  vanishes  outside  the  aperture,  and  equals  the  displacement  of 
the  vortex  sheet  within  the  aperture.  Thus 

v^(x,y)  =  (d/dt  +  U^3/3x)f  =  -i(w  +  iU^S/Sx)^,  (2.2) 


and  integration  by  parts  permits  (2.1)  to  be  cast  in  the  form 


<P.,(x)  = 


iif 


w  +  iU 


r(x’ ,y’ ) 
|x-x’ I 


dx  ’  dy  ’  , 


(2.3) 


where  the  final  integration  may  be  restricted  to  the  area  S  occupied  by 
aperture  provided  -  0(1/|X|^),  a  <  1,  near  the  aperture  edge,  where  X 
denotes  distance  from  the  edge. 


The  volume  flux  Q  through  the  aperture  is  needed  to  calculate  the 
conductivity  from  the  definition  (1.1).  This  is  given  by 
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•i(w  +  iU_^a/ax)$:  dxdy  = 


f (x,y)dxdy. 


(2.4) 


The  equation  determining  r(x,y)  in  terms  of  [p]  is  obtained  from  the 
condition  that  the  pressures  just  above  and  below  the  vortex  sheet  must  be 
equal.  The  linearized  Bernoulli  equation  supplies  the  perturbation  pressure 
in  the  form 

p^  +  ipjj  (w+iU^a/ax)(p^  ,  (z  ^0). 

Continuity  of  pressure  across  the  vortex  sheet  therefore  requires  that 


r(x’ ,y’ ) 
27r|x-x’  I 


dx’dy’ 


(P+-P. )/Po  > 


for  z  =  0,  (x^+y^)^^^  <  R,  (2.5) 


Integrating  with  respect  to  the  second  order  differential  operator  on  the  left 
hand  side,  we  find  (for  constant  (p^-p.  )/P(j) 


Is 


Z(r  ,v’)dVdr,’ 

V{(M’)2+(r?-r7’)2} 


=  1  +  a^(r?)e^^i^  +  a^(ri)e^‘^2^  ,  <  1.  (2.6) 


where  Z  =  {Pp  w^R/?r(p_^ -p_ )  }{■  is  the  displacement  in  dimensionless  form, 

^  =  x/R,  r?  =  y/R,  etc  (so  that  the  integration  is  over  the  unit  disk 

<  1),  and  and  are  the  dimensionless  Kelvin- Helmholtz 
wavenumbers  of  instability  waves  on  the  vortex  sheet  (Lamb  1932) 

a)R(l+i)  wR(l-i) 

a  =  - ,  (7 - .  (2.7) 

^  U^+  iU.  ^  U^-  iU. 

The  amplitudes  a^(ri)  ,  a^in)  of  these  waves  are  constant  along  any  line 
drawn  on  the  sheet  in  the  streamwise  direction,  but  vary  with  transverse 
location  t]  .  They  must  be  chosen  (in  accordance  with  the  Kutta  condition)  to 
ensure  that  the  vortex  sheet  leaves  the  upstream  semi-circular  arc  of  the 
aperture  edge  tangentially,  i.e.,  such  that  Z(^,ri)  =  3Z(^,»j)/3^  =  0  for 
=1,  ^  <  0.  Evidently,  are  even  functions  of  rj . 

Potential  theory  (Kellogg  1954;  Sneddon  1966)  implies  that  f  will  generally 
exhibit  a  mild  singularity  of  order  1/VR-r  at  the  downstream  edge  r  = 

(x^+y^)^^^  =  R,  X  >  0,  of  the  aperture.  The  singularity  is  integrable,  however 
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(yielding  a  finite  value  for  the  volume  flux  Q  (2.4)),  and  is  a  physically 
acceptable,  linear  theory  representation  of  the  large  amplitude  and  violent 
motions  produced  by  vortex  impingement  on  an  edge  (Rockwell  1983).  A 
correction  that  incorporates  the  influence  of  displacement  thickness  waves  on 
the  downstream  surfaces  of  the  plane  (modeling  vorticity  convected  in  the 
boundary  layers  after  ejection  from  the  aperture)  can  be  introduced  to 
eliminate  the  singularity,  but  this  cannot  be  done  in  a  unique  manner  and  is 
not  expected  to  modify  significantly  the  conclusions  of  the  present  simpler 
approach  (Howe  1981b) . 

The  solution  of  the  integral  equation  (2.6)  determines  the  volume  flux  Q 
via  (2.4).  The  definition  (1.1)  then  yields 

Kjj  =  TTRjg  Z(^,ri)d^dr,.  (2.8) 

In  the  absence  of  flow  the  conductivity  is  real  and  equal  to  2R.  In  the  more 
general  case  it  is  convenient  to  set 

Kjj  =  2R(rj^  -  i\),  (2.9) 


where  F  and  A„  are  real  valued  and  depend  on  the  frequency  to. 

R  R 

2.2  Numerical  results 

Numerical  solutions  of  equation  (2.6)  have  been  obtained  by  discretizing 
the  integral.  The  aperture  was  overlaid  by  a  mesh  of  square  integration 
elements  within  each  of  which  Z,  and  were  assumed  to  be  constant,  and 
the  Kutta  condition  imposed  by  requiring  Z  to  vanish  in  those  elements 
overlapping  the  upstream  edge  (where  the  flow  separates)  and  in  each  of  the 
contiguous  elements  just  downstream.  The  application  of  (2.6)  in  each  element 
of  the  mesh  supplies  a  sufficient  number  of  equations  to  determine  the 
discretized  values  of  Z,  and  a2  •  Further  details  of  the  procedure  are 
given  by  Scott  (1995).  In  the  degenerate  case  in  which  the  mean  velocities  on 
each  side  of  the  plane  are  the  same  and  equal  to  U  (so  that  ^2  ^  wR/U)  , 

the  right  hand  side  of  equation  (2.6)  was  taken  in  the  modified  form 

1  +  |^aj^(r7)  +  (*7)?je^‘^^ ,  where  a  =  wR/U. 

The  real  and  imaginary  parts  of  the  conductivity  (2.9)  are  functions  of 
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the  Strouhal  number  uR/U,^  and  the  velocity  ratio  U_/U^.  Their  calculated 
values  are  given  in  Tables  1  and  2  for  the  two  cases  of  (i)  one-sided  flow: 

=  U,  U_  s  0,  and  (ii)  two-sided  flow:  U_|.  =  U_  =  U.  In  applications  these 
tabulated  results  can  be  approximated  by  cubic  splines,  and  this  has  been  done 
in  plotting  and  for  these  two  cases  in  Figure  2. 

For  one-sided  flow  (Case  (i))  r„  =■  4  when  wR/U  becomes  large  and  the 
imaginary  part  ^  0.  \  is  negative  for  1.9  <  wR/U  <3.9.  According  to 

equation  (1.4),  a  negative  amount  of  work  is  then  done  by  the  applied  pressure 
[p]  =  P+ -P  .  i-e.,  energy  is  extracted  from  the  mean  flow.  In  Case  (ii)  (U^  = 
U_ )  the  motion  is  irrotational  when  [p]  =  0.  However,  there  still  exist 
periodically  spaced  frequency  bands  (of  period  n  when  wR/U  >  3)  within  which 
A  <0,  and  where  perturbations  grow  at  the  expense  of  the  mean  flow.  This 
extraction  of  mean  flow  energy  is  not  associated  with  vortex  sheet 
instability,  but  should  in  this  case  be  regarded  as  a  consequence  of  the 
interaction  of  the  vorticity  with  the  semi-circular  "leading  edge"  on  which  it 
impinges  after  convecting  across  the  aperture.  We  shall  show  in  §4  that  cases 
(i)  and  (ii)  are  fundamentally  different,  inasmuch  as  the  aperture  motions  are 
respectively  absolutely  and  conditionally  unstable. 
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Conductivity  of  a  circular  aperture  for  one-sided  grazing  flow 

“  U.  U_  =  0 


wR/U 

coR/U 

Tr 

0.00 

0.000 

0.000 

1.50 

1.949 

0.623 

0.10 

-0.030 

0.005 

1.75 

2.097 

0.244 

0.20 

-0.119 

0.041 

2.00 

2.233 

-0.131 

0.30 

-0.253 

0.145 

2.25 

2.409 

-0.540 

0.40 

-0.395 

0.360 

2.50 

2.708 

-1.006 

0.50 

-0.464 

0.710 

2.75 

3.277 

-1.459 

0.60 

-0.350 

1.140 

3.00 

4.184 

-1.555 

0.70 

-0.016 

1.503 

3.25 

4.887 

-0.969 

0.80 

0.439 

1.673 

3.50 

4.897 

-0.300 

0.90 

0.872 

1.652 

3.75 

4.628 

-0.000 

1.00 

1.216 

1.515 

4.00 

4.398 

0.059 

1.10 

1.466 

1.333 

4.25 

4.258 

0.033 

1.20 

1.644 

1.142 

4.50 

4.189 

0.003 

1.30 

1.774 

0.959 

4.75 

4.175 

0.011 

1.40 

1.872 

0.786 

5.00 

4.182 

0.115 

Table  1 
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2.3  Comparison  with  previous  work 


The  two-sided  flow  problem  where  =  U.  =  U  has  been  investigated  by  Howe 
(1980) ,  who  derived  the  following  integral  equation  for  the  dimensionless 
normal  velocity  v’  =  -i(a  +  id/d^)Z  of  the  vortex  sheet 


I. 


V’  (^’  ,r?’)d^’dr?’ 


<  1, 


(2.10) 


where  C  is  a  known  constant.  The  Kutta  condition  was  imposed  by  choosing  a(ri) 
to  ensure  that  v’  is  finite  over  the  upstream  edge  +r]^  =1,  C  <  0  of  the 

aperture.  Consequently,  v  -  0(l/A^Tr)  at  the  downstream  edge,  on  which  the 
shed  vorticity  impinges  after  convecting  across  the  aperture.  In  fact, 
equation  (2.10)  is  invalid  because  the  range  of  integration  in  (2.10)  cannot 
be  restricted  to  the  region  S  of  the  aperture  unless  account  is  taken  of  a 
singular  contribution  to  v’  from  the  downstream  edge.  When  v  -  0(l/VR-r)  at 
this  edge,  the  displacement  f  will  have  a  simple  discontinuity  there,  with 
limiting  value  ,  say,  as  r  R-0  (x  >  0)  ,  and  the  singular  velocity  at  the 
edge  would  equal  -U(*jj5(r“R)  .  This  should  appear  as  a  separate  contribution  in 
(2.10)  and  must  be  included  in  the  evaluation  of  the  volume  flux  Q.  A  finite 
limiting  behavior  of  f  at  the  downstream  edge  is  possible  only  when  =  U_ . 

since  in  general  the  edge  singularity  is  much  stronger  than  a  simple 
discontinuity,  being  proportional  to  1/VR-r ,  although  the  aperture  flux 
remains  finite. 
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3.  THE  RECTANGULAR  APERTURE 

It  is  of  interest  to  compare  the  numerical  predictions  for  a  circular 
aperture  with  corresponding  results  that  can  be  derived  in  analytical  form  for 
a  rectangular  aperture  of  large  aspect  ratio. 


Take  the  mean  flow  in  the  x-direction  with  velocities  U^  above  and  below 
the  plate,  as  before,  and  let  the  rectangular  aperture  occupy  the  region  |x|  < 
s.  Iyl  2^>  where  the  aspect  ratio  b/2s  >  1  (Figure  3).  In  the  notation  of 
§2,  the  vortex  sheet  displacement  f  satisfies  the  integro-differential 
equation  (2.5),  in  which  S  now  denotes  the  rectangular  area.  An  approximate 
solution  of  this  equation  can  be  obtained  when  b/2s  is  large  by  assuming  the 
motion  of  the  vortex  sheet  to  be  two-dimensional,  i.e.,  by  neglecting  the 
dependence  of  r(x’,y’)  on  the  spanwise  coordinate  y’ .  The  integration  with 
respect  to  y’  over  (-|b,|b)  in  (2.5)  may  then  be  performed  explicitly.  If  we 
also  integrate  over  -|b  <  y  <  |b,  the  resulting  equation  may  be  regarded  as 
satisfying  the  condition  of  pressure  continuity  across  the  vortex  sheet  in  a 
spanwise  averaged  sense. 


To  do  this  observe  that,  when  x  and  x’  both  lie  within  the  aperture,  so 
that  |x-x’ I  <  2s , 


II 


.b/2 


dy  *  dy 


-b/2  )^+(y-y’ 


=  2b  sinh'^  (b/|x-x’ I  )  -  2{Vb2-i-(x-x’  )2  -  |x-x’|}. 


=  2b^ln{2b/e |x-x’ I }  + 


|x-x’ I  |x-x’ I  2 

~b  4^2" 


...  ).  (3. 


1) 


where  e  is  the  exponential  constant.  When  b/2s  >  1  only  the  first  term  in  the 
brace  brackets  of  this  expansion  need  be  retained.  Making  this  approximation 
in  (2.5),  and  integrating  with  respect  to  the  differential  operator,  we  obtain 


2b 


1  + 


Kl  <  1. 


(3.2) 


where 


Z  =  {2p(,w2s/7r(p^-p.)}f ,  ?  =  x/s,  C  =  x’/s 
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and  o  and  o  are  nondimensional  Kelvin-Helmholtz  wavenumbers  defined  as  in 
1  2 

(2.7)  with  R  replaced  by  s.  The  coefficients  and  are  constants  whose 
values  are  fixed  by  the  Kutta  condition  Z  =  3Z/3^  =  0  at  ^  =  -1. 

The  solution  of  (3.2)  satisfying  these  conditions  has  been  given  by  Howe 
(1981b).  Z  varies  like  (^+1)^^^  near  the  upstream  edge  ^  =  -1+0  of  the 
aperture,  and  is  singular  like  l/(l-0^^^  at  the  downstream  edge  (C  ^  1-0) 
the  vortex  sheet  impinges  on  the  solid  surface.  The  volume  flux 
Q  =  {7rb(p^ -p_  )/2ippW}J^^  Z(Od?. 

and  the  Rayleigh  conductivity  is  therefore  equal  to  |b/^^Z(^)dC.  Using  the 
solution  given  by  Howe  (1981b)  we  find 

K  _ _ — _  (3.3 

^  2{F(a^ ,a^)  +  $} ’ 


where 


-o^Jq(a3)[Jp(g,)-2W(a^)]  +  Jp  (^^  )  [  Jq  (^2  ^ ^ 

CTiW(a2)[Jo(ai)-2W(c7^)]  -  C7^yi{a^)[3 -2Vi{a^)] 


^  =  ln(4b/es) , 


(3.4) 


W(x)  =  ix[Jg(x)  -  iJ^(x)],  and  and  are  Bessel  functions. 


The  conductivity  in  the  absence  of  flow  is  also  given  by  (3.3)  by 
formally  setting  F  =  0.  For  that  case  we  can  write 

Kj^/2R^  =  (jr^b/32s)^'Vln{4b/sc},  (3.5) 

where  R^  =  (2sb/7r)^^^  is  the  radius  of  the  circle  of  area  2sb  (that  of  the 
rectangle).  According  to  Rayleigh  (1945)  the  circular  aperture  has  the 
minimum  conductivity  of  all  apertures  of  the  same  area,  and  indeed  the  right 
hand  side  of  (3.5)  has  a  minimum  value  of  about  1.1,  and  does  not  exceed  1.2 
when  2.5  <  b/ s  <  12.8. 


Figure  4(i)  illustrates  the  dependence  of  the  real  and  imaginary  parts  of 
K^(w)/2R  =  on  a  =  ws/U  for  one  sided  flow,  =  U,  U_  =0  when  the 

aspect  ratio  b/2s  =  10.  These  curves  are  qualitatively  the  same  as  the 
corresponding  plots  in  Figure  2(i)  for  the  circular  aperture.  It  follows  from 
equations  (3.3),  (3.4)  for  the  rectangular  aperture  that 
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(7r^b/32s) 


1/2 


®  ln(4b/se^) 


=  0,  a 


r  =  -3a2(7r^b/32s)^'2_  (7r^b/32s) a 

R  ^ 


3  /_3i 


,1/2 


0. 


(3.6) 


Thus,  at  high  frequencies  the  conductivity  is  real  and  exceeds  that  in  the 
absence  of  the  vortex  sheet  (provided  4b/se^  >  1) . 


For  two-sided  flow  with  U_^  =  U_  =  U,  there  is  no  vortex  sheet  in  the 
steady  state  and  a ^  ^  ^  a  =  ws/U.  The  corresponding  limiting  form  of 

function  F  of  (3.4)  is 


F(a,a) 


iaJ^(J^-iJ^)  -  [Jp-2ia(J^-iJ^)][J„-ia(J„+iJ,)] 

a[JQj^+a{j2-h(Jjj-2iJ^)M] 


(3.7) 


where  the  Bessel  functions  are  evaluated  at  a  =  ws/U.  The  variations  of 
and  A  with  a  are  shown  in  Figure  4(ii)  for  b/2s  =  5.  Once  again,  the  results 
are  qualitatively  similar  to  those  in  Figure  2(ii)  for  the  circular  aperture. 
However,  whereas  Figure  2(ii)  indicates  that  the  oscillatory  behavior  of  Kj^ 
ultimately  dies  out,  and  that  ->  1,  Aj^  -»  0  as  u  ^  oo,  the  limiting  behaviors 
for  the  rectangular  aperture  remain  oscillatory,  with  period  ;r ,  since 
F  =  4/{5sin(2a)-h4icos(2a) -3}  when  a  >  2.  ks  a  -*  0  we  find 


6a^  (7r^b/32s)^'2. 


(3.8) 
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4.  ABSOLUTE  AND  CONDITIONAL  INSTABILITY 

The  volume  flux  Q  is  determined  by  the  applied  pressure  [p]  s  p^ -p_  and 
the  conductivity  .  All  of  these  quantities  are  generally  functions  of  the 
frequency  cj,  and  the  real  form  of  (1.1)  is 

00 

p^dQ(t)/dt  =  -J  Kj^(w)  [p(w)  dw.  (4.1) 

-00 

The  integration  is  nominally  taken  over  all  real  values  of  w.  However,  if  the 
integrand  has  singularities  in  the  upper  half  of  the  complex  w-plane,  the 
integration  contour  must  pass  above  them  to  ensure  that  causality  is 
satisfied.  Since  [p(t)]  is  arbitrary,  and  may  be  assumed  to  vanish  prior  to 
some  initial  instant,  any  singularities  in  Im(a))  >  0  must  be  associated  with 
Kj^(w).  When  they  are  present  the  system  is  unstable,  and  an  unsteady  flux 
through  the  aperture  can  develop  spontaneously.  Also,  because  Q(t)  and  [p(t)] 
are  real  valued  quantities,  it  follows  that  Kj^(-w  )  =  Kj^(w),  where  the 
asterisk  denotes  complex  conjugate. 

Consider  the  stable  case  in  which  (w)  =  2R{rj^  (w)  -  iA^,  (w) }  is  regular  in 
Im{w)  >  0.  For  real  u  we  have 

rjj(-c^)  =  rj^(w),  Aj^(-w)  =  -Aj^(w).  (4.2) 

In  general  (w)  =  constant  as  |w|  “  and  Aj^(w)  -♦  0.  The  function 

/(w)  =  r  (w)-r  -iA^(w)  is  regular  in  Im(w)  >  0  and  vanishes  as  |w|  -> 

Cauchy’s  theorem  (Carrier,  et  al  1966)  applied  in  the  usual  way  to  a  closed 
contour  consisting  of  the  real  axis  and  a  semi -circle  in  the  upper  half -plane 
whose  radius  grows  without  limit,  then  implies  that 

/(co)  =  (l/27ri)J^^/(Od?/(^-^^)  >  0, 

where  the  integration  is  along  the  real  axis.  The  real  and  imaginary  parts  of 
this  expression  as  Im(w)  +0  yield  the  Kramers -Kronlg  formulae 


r  (cj)  -  r  - 

R  ^  ^  00 


2  r“  \A^(X)dX 
TT  Jo  A^-CO^ 


Aj^(w) 


^  £  [rR(A)-r„]dA 

TT  Jo  A2-w2 


(w  real) 


where  the  integrals  are  principal  values. 


(4.3) 
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When  the  motion  is  stable  these  equations  can  be  used  to  determine  , 
say,  when  is  known.  If  vanishes  at  zero  frequency  (as  it  would 

when  vortex  shedding  completely  blocks  the  aperture  flux) ,  the  constant  F^  is 
found  by  setting  to  =  0  in  the  first  of  (4.3). 

Alternatively,  when  r^(co)  and  \  are  known  for  real  w,  equations  (4.3) 
can  be  used  to  investigate  the  stability  of  the  flow.  If  the  equations  are 
not  satisfied  (w)  must  be  singular  in  Im(ci))  >  0,  and  the  motion  is 
absolutely  unstable .  This  is  the  case  for  one-sided  grazing  flow  over  the 
circular  and  rectangular  apertures  of  §§3,  4,  and  is  presumably  a  consequence 
of  vortex  sheet  instability.  According  to  Figures  2(ii)  and  4(ii),  <  0  for 

discrete  ranges  of  real  frequencies  for  two-sided  flow  when  =  U_  =  U. 
Equation  (1.4)  indicates  that  energy  is  extracted  from  the  mean  flow  at  these 
frequencies.  However,  the  Kramers -Kronig  formulae  are  satisfied  by  the 
numerical  data  for  two-sided  flow  in  Figures  2(ii)  and  4(ii) .  The  dots  in 
these  figures  represent  values  of  F^^  computed  from  the  first  of  the 
Kramers -Kronig  equations  (4.3)  using  the  respective  values  of  shown  in  the 
figures  (taking  Aj^(A)  =  0  for  A  >  20)  .  The  motion  is  therefore  conditionally 
unstable:  energy  is  extracted  from  the  mean  flow  only  when  the  applied 
pressure  [p]  contains  a  component  with  frequency  within  an  interval  where  A^^  < 
0. 

When  the  aperture  is  backed  by  a  cavity  (in  which  U_  =  0)  whose 
dimensions  are  large  compared  to  the  aperture  diameter,  the  instability  of  the 
shear  layer  can  excite  and  sustain  cavity  acoustic  resonances.  This  can  occur 
only  when  the  frequency  lies  within  the  range  where  Aj^(w)  <  0.  For  the 
circular  aperture  it  is  usual  to  express  the  frequency  at  which  a  resonance  is 
excited  in  terms  of  the  Strouhal  number  /D/U,  where  /  —  co/2n  and  D  is  the 
diameter.  According  to  Figure  2(i)  flow  excited  cavity  resonances  must  occur 
for  0.6  <  /D/U  <  1.15.  This  range  is  consistent  with  experimentally  observed 
Strouhal  numbers  for  laminar,  high  Reynolds  number  flow  over  cavity  openings 
in  cases  where  the  shear  layer  is  well  approximated  by  a  vortex  sheet,  i.e., 
when  the  shear  layer  thickness  is  very  much  smaller  than  D  (see,  e.g.,  DeMetz 
and  Farabee  1977). 
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5.  CONCLUSION 

The  Rayleigh  conductivity  of  an  aperture  in  a  thin  plane  determines 
the  volume  flux  through  the  aperture  generated  by  a  small  amplitude,  time 
harmonic  pressure  load.  When  the  motion  is  incompressible  and  irrotational 
is  real  and  constant,  and  depends  only  on  the  shape  of  the  aperture.  The 
conductivity  is  complex  for  a  real  fluid.  If  there  is  a  tangential  flow  over 
the  plane,  and  the  Reynolds  number  based  on  aperture  diameter  is  large,  the 
pressure  load  can  be  assumed  to  generate  vorticity  at  the  aperture  edge,  which 
causes  perturbation  energy  to  be  absorbed  or  generated  by  the  mean  flow, 
according  as  Im{Kjj}  ^  0. 

In  this  paper  Kj^  has  been  calculated  numerically  for  tangential  flow  over 
a  circular  aperture  of  radius  R  by  modeling  the  shear  layer  as  a  vortex  sheet. 
For  one-sided  flow  at  speed  U,  energy  is  extracted  from  the  mean  flow  when  the 
Strouhal  number  wR/U  is  between  about  1.9  and  3.9.  The  aperture  motion  is 
absolutely  unstable,  in  the  sense  that  a  small  perturbation  in  the  mean  stream 
can  trigger  motion  of  the  vortex  sheet  that  becomes  unbounded  on  the  basis  of 
linear  theory.  When  the  mean  velocity  is  the  same  on  both  sides  of  the  plane 
the  flow  is  irrotational  in  the  steady  state.  Vorticity  is  generated  by  an 
applied  pressure  (in  accordance  with  the  Kutta  condition)  and  there  are 
discrete,  quasi -periodic  frequency  intervals  within  which  perturbation  energy 
is  alternately  absorbed  and  generated  by  the  mean  flow.  However,  this  motion 
is  conditionally  unstable,  in  that  perturbation  energy  grows  at  the  expense  of 
the  mean  flow  only  when  the  applied  pressure  is  non-zero  and  of  the 
appropriate  frequency. 

The  numerical  results  for  a  circular  aperture  are  in  qualitative 
agreement  with  analytical  predictions  for  a  rectangular  aperture  of  large 
aspect  ratio,  whose  long  edge  is  transverse  to  the  mean  flow  direction. 
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CHAPTER  4 


DAMPING  OF  SOUND  AND  VIBRATION  BY  FLOW  NONLINEARITY  IN 
THE  APERTURES  OF  A  PERFORATED  ELASTIC  SCREEN 
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SUMMARY 

An  analysis  is  made  of  the  influence  of  flow  nonlinearity  in  the 
apertures  of  a  perforated  elastic  plate  on  the  damping  of  sound  and  flexural 
vibrations.  Fluid  is  forced  through  the  perforations  by  the  pressure 
differential  established  across  the  plate  by  the  incident  disturbance.  The 
Reynolds  number  is  assumed  to  be  sufficiently  large  that  separation  occurs, 
and  the  reciprocating  aperture  flows  form  "jets"  on  alternate  sides  of  the 
plate.  The  growth  of  these  jets  is  modeled  by  means  of  a  nonlinear  equation 
proposed  by  Cummings  (1986).  This  equation  is  solved  simultaneously  with  a 
generalized  bending  wave  equation  derived  by  the  author  (Howe  1995a)  which 
governs  motions  of  a  perforated  elastic  plate  whose  lengths  scales  are  large 
compared  to  the  aperture  spacing.  It  is  shown  that  significant  attenuations 
of  large  amplitude  acoustic  waves  can  occur  except  when  the  frequency  is  so 
small  that  the  plate  is  acoustically  transparent.  Bending  waves  are  also 
damped  provided  the  amplitude  of  the  plate  surface  velocity  is  not  too  large 
and  the  frequency  is  small  enough  to  ensure  the  formation  of  substantial  jots 
in  the  apertures.  Numerical  results  are  given  for  large  amplitude  sound  waves 
incident  on  a  perforated  screen  in  air,  and  for  bending  waves  propagating  over 
aluminium  and  steel  screens  immersed  in  either  air  or  water. 
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1.  INTRODUCTION 

When  a  sound  wave  impinges  on  a  solid  surface  in  the  absence  of  mean  flow 
dissipation  usually  occurs  through  the  direct  action  of  viscous  and  molecular 
thermal  diffusion.  At  very  high  acoustic  amplitudes,  however,  free  vorticity 
is  formed  at  sharp  edges ,  and  dissipation  takes  place  by  the  transfer  of 
energy  to  the  kinetic  energy  of  the  vorticity,  which  is  convected  away  from 
the  interaction  region  by  self-induction.  This  nonlinear  mechanism  can  be 
particularly  important  when  sound  generates  motions  in  small  perforates  or 
apertures;  when  the  acoustic  Reynolds  number  is  large  (w  being  the 

acoustic  frequency,  R  the  aperture  radius,  and  i/  the  kinematic  viscosity)  it 
is  independent  of  the  viscosity,  except  insofar  as  viscosity  must  be  called 
into  play  to  permit  vorticity  to  be  released  into  the  fluid  (Sivian  1935; 
Thurston,  Hargrove  &  Cook  1957;  Ingard  &  Ising  1967;  Zinn  1970;  Mailing  1973; 
Cummings  1984,  1986;  Salikuddin  1990;  Salikuddin  &  Brown  1990).  In  addition, 
aperture  nonlinearity  typically  causes  sinusoidally  varying  incident  sound  to 
produce  periodic  reflected  and  transmitted  waves  containing  multiple  harmonics 
of  the  incident  wave  frequency. 

This  mechanism  of  attenuation  by  kinetic  energy  production  has  long  been 
exploited  to  increase  the  damping  of  moderate  amplitude  sound  by  the  simple 
expedient  of  introducing  a  high  Reynolds  number  mean  flow  (Barthel  1958; 
Bechert,  Michel  &  Pfizenmaier  1977;  Bechert  1979;  Howe  1979,  1980a,  b; 

Rienstra  1981;  Cargill  1982;  Hughes  &  Dowling  1990;  Guo  1991).  For  example, 
the  absorption  of  acoustic  energy  by  a  perforated  screen  is  enhanced  when  a 
"bias  flow"  is  maintained  through  the  apertures.  Viscous  effects  are 
important  only  in  the  immediate  vicinity  of  the  surface,  where  vorticity  is 
produced  by  the  acoustic  oscillations;  the  vorticity  is  subsequently  swept 
away  by  the  mean  flow  and  its  kinetic  energy  is  permanently  lost  to  the 
incident  sound.  When  the  fluid  loading  is  small,  such  as  for  a  steel  plate  in 
air,  it  is  usually  permissible  to  estimate  the  damping  by  assuming  the  screen 
to  be  rigid,  but  for  plates  immersed  in  liquids  the  loading  is  often  large 
enough  that  high  acoustic  intensities  are  accompanied  by  significant 
structural  vibrations.  Vorticity  production  is  then  modulated  both  by  the 
sound  and  the  surface  motions.  Thus,  a  bias  flow  may  be  expected  to  cause 
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enhanced  absorption  of  structural  vibrations  as  well  as  sound. 

An  experimental  investigation  is  currently  in  progress  at  Boston 
University  to  quantify  the  influence  of  grazing  and  bias  flows  on  the 
absorption  of  structural  vibrations  by  a  perforated  plate.  In  this  paper  the 
importance  of  nonlinear  aperture  motions  on  damping  is  examined  in  the  absence 
of  mean  flow.  An  objective  is  to  assess  the  likely  contribution  of  aperture 
nonlinearity  to  the  observed  attenuations.  A  generalized  bending  wave 
equation  governing  the  propagation  of  small  amplitude  flexural  motions  of  a 
thin,  fluid  loaded,  perforated  elastic  plate  has  previously  been  applied  (Howe 
1995a,  b)  to  determine  the  influence  of  mean  flow  (either  a  bias  flow,  or  a 
tangential  mean  flow)  on  damping.  This  equation  is  also  applicable  in  the 
absence  of  mean  flow  and  for  nonlinear  aperture  motions,  provided  the  plate 
motion  is  linear.  Cummings  (1986)  has  proposed  a  simple  and  convenient  system 
of  nonlinear  equations  for  fluctuating  forced  motions  in  small  circular 
apertures  in  the  absence  of  mean  flow,  predictions  of  which  have  been 
validated  by  comparison  with  experiment.  These  equations  are  here  adapted  and 
used  in  conjunction  with  the  generalized  bending  wave  equation  to  investigate 
the  interaction  of  long  wavelength  sound  and  flexural  vibrations  with 
apertures  in  a  thin  elastic  plate.  In  particular,  we  consider  the  problems  of 
the  transmission  of  sound  through  a  perforated  plate,  and  the  damping  of 
bending  waves. 

The  governing  equations  are  reviewed  in  §2.  The  sound  transmission 
problem  is  discussed  in  §3,  where  a  comparison  is  made  with  corresponding 
predictions  for  a  rigid  screen.  In  §4  we  investigate  the  interaction  of  a 
bending  wave  on  a  fluid  loaded  elastic  plate  with  a  single,  isolated  aperture, 
and  make  extension  to  bending  waves  propagating  over  a  homogeneously 
perforated  plate. 
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2.  THE  GOVERNING  EQUATIONS 


2.1  Nonlinear  theory  of  xmsteady  high  Reynolds  n\imber  flow  through  a  circular 
aperture 

Consider  a  rigid  wall  of  thickness  h  containing  a  circular  cylindrical 
aperture  of  radius  R.  The  wall  is  immersed  in  fluid  of  uniform  mean  density 
p  and  sound  speed  c  ,  and  a  time-periodic  pressure  differential  [p]  = 
P2(t)-p^(t)  (of  radian  frequency  w  >  0)  is  applied  across  the  wall  as 
illustrated  schematically  in  Figure  1.  The  mean  value  of  [p]  is  assumed  to 
vanish,  so  that  there  is  no  mean  flow  through  the  aperture;  the  aperture 
radius  R  and  wall  thickness  h  are  assumed  to  be  much  smaller  than  the 
characteristic  acoustic  wavelength  2n/K^  ,  where  =  w/c^  is  the  acoustic 
wavenumber.  During  the  interval  of  time  in  which  p.^  >  p^  ,  fluid  is  forced 
through  the  aperture  from  left  to  right  in  the  figure.  The  Reynolds  number  is 
assumed  to  be  sufficiently  large  that  the  flow  through  the  aperture  separates 
and  forms  a  jet. 


To  motivate  the  following  discussion,  the  motion  within  the  jet  will 
first  be  assumed  to  be  irrotational ,  and  the  pressure  on  the  free  surface  of 
the  let  is  taken  to  be  p„ .  Let  V  (t)  denote  the  mean  axial  velocity  of  the 
jet  where  it  enters  the  aperture  (at  0  in  the  figure).  Bernoulli’s  equation 
applied  to  the  axial  streamline  within  the  jet  between  0  and  the  point  marked 
J  just  within  the  potential  flow  region  at  the  head  of  the  jet,  where  the 
velocity  is  and  the  pressure  p^ ,  implies  that 


d 

at 


+  -v^ 

2  J 


2  A 


Pa  ~P2 
f>o 


(2.1) 


where  cp  and  cp  are  the  values  of  the  velocity  potential  at  0  and  J.  The  flow 
^  A  J 

into  the  aperture  may  be  regarded  as  entirely  irrotational,  so  that  we  also 
have 


— {<p  -(p  )  +  xVj 

2  A 


Pi -Pa 
Po 


(2.2) 


where  the  motion  at  large  distances  (>  R)  from  the  aperture  in  the 
irrotational  domain  is  assumed  to  be  small  enough  to  be  governed  by  linearized 


equations,  i.e.,  »  -p^dcp^/dt. 
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Rayleigh’s  formula  (Rayleigh  1945)  can  be  used  to  write 
where  =  0.82R  is  the  "end-correction"  of  the  circular  opening  into  the 
irrotational  region.  Also, 

dcpj/dt  -  d(p^/dt  =  av(C,t)/at  d^, 


where  v(z,t)  is  the  velocity  on  the  axis  of  the  jet  at  distance  z  from  the 
inlet  plane.  The  value  of  the  integral  in  this  formula  is  denoted  by 
2’(t)aV^/at,  where  i’ (t)  is  a  suitable,  time  dependent  length. 


Inserting  these  expressions  into  equations  (2.1),  (2.2)  and  adding,  we 


find 


,(t)^A  ^  iv2  =  ll£l 
^  ^at  2  J  p 


(2.3) 


and 


where  i(t)  =  +  .^’(t).  Equation  (2.3)  describes  the  motion  of  a  slug  of 

fluid  of  variable  effective  length  i(t)  and  variable  mass  p^A^£it)  (where 
Ag  =  ttR^  is  the  area  of  the  aperture)  subject  to  a  driving  pressure  p^ -p, 
a  resistive  force  -- p  A  .  In  the  half-cycle  during  which  the  flow  through 
the  aperture  is  from  right  to  left  in  Figure  1,  equation  (2.3)  remains  valid 
provided  the  sign  of  the  nonlinear  resistive  term  is  reversed.  The  velocity 
Vj  can  be  eliminated  by  writing  =  V^(t)/r;,  where  rj  is  a  contraction 
coefficient.  By  this  means  we  arrive  at  the  following  equation  describing 
reciprocating,  periodic  flow  through  the  aperture  (Cummings  &  Eversman  1983; 
Cummings  1986) 

VaIVJ  -[p] 


2r,‘ 


(2.4) 


If  the  flow  were  entirely  irrotational  (no  jet  formation)  the  length  i(t) 
would  be  constant  and  equal  to  ^2  *  ^j^+h  =  0.82R  +  h  (Rayleigh 

1945).  When  separation  occurs  i(t)  can  be  expressed  in  terms  of  an  effective 
length  of  the  jet  defined  by 


(2.5) 


where  the  time  t  is  measured  from  the  beginning  of  the  half  cycle  during  which 
the  sign  of  V  (t)  is  constant.  From  an  analysis  of  data  derived  from  aperture 
flow  experiments  performed  by  Ingard  (1953)  ,  Thurston,  Hargrove  &  Cook  (1957) 
and  Ingard  &  Ising  (1967),  in  which  the  working  fluid  was  either  water  or  air. 
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Cuminlngs  (1986)  deduced  the  following  empirical  relation  between  the  lengths 


.2(t)  and  L, 


i(t)  =  A  + 


1  1+/9(L/D)  1-585  ’ 


P  = 


D  =  2R. 


(2.6) 


There  is  no  corresponding  formula  for  the  dependence  of  the  contraction 
ratio  ri  on  T .  This  would  be  expected  to  vary  significantly  only  near  the 
beginning  of  a  half  cycle  when,  however,  the  nonlinear  resistance  is  small. 
Experiments  for  steady  flow  suggests  that  »j  =  0.61,  but  in  a  careful 
comparison  of  predictions  of  equations  (2.4)  -  (2.6)  with  experiments 
involving  high  amplitude  sound  incident  on  small  circular  apertures,  Cummings 
(1986)  has  concluded  that  r;  may  be  regarded  as  constant  and  equal  to  0.75. 

2.2  The  generalized  bending  wave  equation 

Consider  next  a  thin  elastic  plate  of  bending  stiffness  B,  mass  density 
=  m/h  (m  being  the  mass  of  unit  area  of  the  plate)  and  Poisson’s  ratio  a, 
perforated  with  small  circular  apertures  of  equal  radius  R  which  are 
distributed  uniformly  with  N  per  unit  area.  In  the  undisturbed  state,  let  the 
median  plane  of  the  plate  occupy  the  plane  z  =  0  of  the  rectangular  coordinate 
system  (x,y,z).  The  fractional  open  area  a  =  NttR^  is  assumed  to  be  small,  so 
that  the  average  distance  d  (-  1/VN)  between  neighboring  apertures  is  very 
much  greater  than  R.  When  a  1  the  average  equation  governing  motion  of  the 
fluid-coupled  plate  whose  length  scale  is  large  relative  to  d  can  be  cast  in 
the  form  (Howe  1995a) 

where  s  (3^/3x^+a^/3y^ )^ ,  Tp .  respectively  denote  the  displacement  of 
the  plate  and  the  mean  fluid  displacement  in  the  apertures,  both  measured  in 
the  x^ -direction,  [p]  =  p(x,y,-f0,t)  -  p(x,y,-0,t)  is  the  pressure  jump  across 
the  plate,  and  all  quantities  are  averaged  over  a  region  of  the  plate 
containing  many  apertures.  The  flexural  motions  of  the  plate  are  small,  and 
governed  by  linear  theory,  but  the  fluid  motion  in  the  apertures  can  be  large, 
and  satisfy  in  a  first  approximation  (Howe  1995a)  equations  (2.4)  -  (2.6), 
wherein  =  V^(x,y,t)  =  3f^/3t. 
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3.  ABSORPTION  OF  SOUND  AT  A  PERFORATED  ELASTIC  SCREEN 
3 . I  Equations  of  motion 

The  equations  of  §2  are  first  applied  to  estimate  the  influence  of 
flexural  motions  on  the  transmission  of  sound  through  a  perforated  screen  (the 
average  aperture  spacing  d  being  small  compared  to  the  characteristic 
wavelength  of  the  sound).  Consider  the  simplest  case,  depicted  in  Figure  2, 
in  which  a  plane,  time-harmonic  acoustic  wave  of  pressure  p  =  PjCOs{w(t-z/c^ ) } 
(where  p^  is  constant  and  w  >  0)  is  incident  from  z  <  0  on  a  thin  elastic 
screen  of  thickness  h,  whose  median  plane  is  z  =  0.  The  transmitted  and 
reflected  acoustic  pressures  at  distances  |z|  >  d  from  the  screen  must  also  be 
plane  periodic  waves,  so  that 

p  =  PjCos{(j(t-z/c^ ) }  +  Pj^(t+z/c^),  z  <  0, 

=  Pj (t-z/c^) ,  z  >  0,  (3.1) 

where,  because  of  the  nonlinear  motions  in  the  apertures,  the  transmitted  and 
reflected  waves  p^(t+z/c^)  and  pj^(t+z/c^)  are  generally  not  simple  harmonic. 

Equation  (2.10)  supplies  the  functional  relations 

p^(t)  =  p^c^V(t),  Pjj(t)  =  PjCOs(a;t)  -  p^c^V(t),  (3.2) 

which  determine  the  reflected  and  transmitted  waves  in  terms  of  the  aggregate 
normal  velocity  V  of  the  perforated  plate,  and  yield 

[p]  =  2pjCos(wt)  -  2p^c^V(t). 

Substituting  this  expression  for  [p]  in  equations  (2.3),  (2.7),  and  expressing 
and  fp  in  terms  of  V  and  ,  the  equations  governing  the  fluid- structure 
interaction  reduce  to 


^  VJVJ  ,  _  2pi 


^(t)— + 

at 


+  2c  V(t)  =  — -  cos(wt), 

2r?2  0  p 


m—  +  2p  c  V  -  a(m  +  2Rp  )^^  =  2p  cos(wt), 

at  °  °  °  at  ^ 


(3.3) 


where  i(t)  is  determined  by  equation  (2.6),  in  which  L  is  the  discontinuous 
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function  defined  by  the  differential  form  of  (2.5): 

2T 

—  =  |V  I ,  L  =  0  when  V  =  0.  (3.4) 

at  '  a'  ’  A 

The  periodic  solution  of  these  equations  must  be  determined  numerically 
for  prescribed  values  of  the  incident  pressure  amplitude  p^  and  frequency  w. 
This  is  done  by  assuming  that  p^  vanishes  for  t  <  0,  that  ^  =  ^  t  =  0, 

and  integrating  for  t  >  0  until  the  solution  becomes  stably  periodic,  usually 
after  three  or  four  cycles. 

The  acoustic  power  incident  on  unit  area  of  the  screen  is 

-  P?/2p  c 

when  averaged  over  one  cycle  27r/w.  Similarly,  the  mean  reflected  and 
transmitted  sound  powers  per  unit  area  are  given  respectively  by 

2 

Jl^  =  J  ^pjCOsCwt)  -  p^c^V(t)]  dt,  I  v2(t)  dt,  (3.5) 

^  o  o 

where  the  integrations  are  over  one  acoustic  cycle. 

An  absorption  coefficient  may  also  be  defined  by 

A  =  1  -  f  V(t)[p  cos(a,t)  -  p  c  V(t)]  dt,  (3.6) 

Hj  7rp2  J  I  °  ° 

which  is  the  fractional  acoustic  power  absorbed  by  vorticity  production  at  the 
screen. 

3.2  Numerical  results 

It  is  convenient  to  introduce  dimensionless  velocities  and  time  as 
follows: 


V  =  V/(wR),  =  V^/(a)R);  T  =  cot/27r.  (3.7) 

Let  the  fluid  be  an  ideal  gas  of  mean  pressure  p^  =  p^c^/y,  where  7  is  the 
ratio  of  specific  heats.  The  incident  wave  amplitude  p^.  may  then  be 
normalized  with  respect  to  p^ ,  and  equations  (3.3)  -  (3.6)  become 
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M  5V.  TT  -  -  ,  47rV 

R  ai  ,,2  a'  a'  (^^r) 


47r(Pi/Po) 

7(/c^R)2 


COS  (27rT) 


3V  47rV  1 

-|_  1 

fRpo] 

1  '  ^  1 

2Rp„'| 

1  + 

aT  (k^R)' 

^  in  - 

c  m 

laT 

47r(Pi/Po)fRPo1 

7 (/c^ R)  2  I  m  J 


cos  (27rT) 


y  (3.8a,  b,  c) 


d 


;r|VJ  .  D  =  2R. 


^  r  V(T)  f^^^^^°^cos(27rT)  -  V(T)1  dT,  (3.9) 

(Pl/Po)^  -*0  ^  ^ 

where  k  =  w/c  is  the  acoustic  wavenumber  of  the  fundamental  mode,  and  i/R  is 

O  '  O 

defined  by  (2.6)  with  =  0.82R,  A^  =  0.82R  +  h.  In  this  form  the  equations 
clearly  exhibit  the  dependence  of  the  fluid- structure  interaction  on  the 
dimensionless  quantities 

Pl/Po  .  Po/Ps  >  “• 

Figures  3-5  depict  typical  numerical  predictions  for  an  aluminium  plate 
in  air  with  fractional  open  area  a  =  0.03  and  R  =  lOh.  The  thick  and  thin 
curves  in  Figure  3  represent  the  respective  contributions  of  the  components 
aV^  and  (l-a)Vp  of  the  normal  velocity  V  (i.e.,  the  respective  contributions 
from  the  aperture  fluxes  and  the  plate  motion)  for  k^R  =  0.1  and  different 
incident  wave  amplitudes.  Each  pair  of  curves  is  plotted  to  the  same  scale 
for  a  given  value  of  Pj/P^.  but  the  scales  for  different  values  of  Pj/p^  are 
not  the  same.  At  low  acoustic  amplitudes  the  motion  is  dominated  by  the 
aperture  flows,  which  are  essentially  simple  harmonic;  the  very  much  smaller 
amplitude  plate  motions  are  seen  to  be  in  phase  with  the  aperture  velocity. 

As  the  acoustic  amplitude  increases,  higher  harmonic  components  of  the 
fundamental  frequency  make  progressively  larger  contributions  to  the  aperture 
flux,  and  the  plate  motion  becomes  larger,  but  still  simple  harmonic. 

However,  the  phases  of  the  plate  and  aperture  motions  diverge  as  Pj/p^ 
increases,  the  difference  being  close  to  90°  when  Pj/P^  =  0.1. 

The  dependence  of  the  absorption  coefficient  A  on  the  acoustic  amplitude 
is  plotted  in  Figure  4  for  several  values  of  k^R  (the  range  of  values  of  Pj/p^ 


111 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


in  the  figure  corresponds  to  sound  pressures  between  about  110  -  170  dB  (re 
2x10"^  N/m^)).  The  curves  for  /c^R  >  0.05  are  essentially  identical  to  those 
for  a  rigid  perforated  plate  of  the  same  thickness  (which  can  be  calculated  by 
letting  p^R/m  ^  0  in  equations  (3.8)).  At  lower  frequencies,  however,  the 
elastic  plate  becomes  progressively  more  transparent  to  the  incident  sound 
and,  whereas  the  rigid  plate  damping  tends  to  a  finite  limit  for  each  value  of 
K  R,  the  elastic  plate  damping  ultimately  decreases  to  zero.  This  is  clear 
from  Figure  5,  which  compares  the  absorption  coefficients  for  elastic  and 
rigid  screens  as  a  function  of  frequency  for  two  values  of  p^/p^. 
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4.  NONLINEAR  DAMPING  OF  BENDING  WAVES 
4.1  Absorption  by  a  single  aperture 

A  flexural  wave  propagating  over  a  thin,  fluid  loaded  elastic  plate  of 
thickness  h  is  accompanied  by  a  traveling  pressure  field  which  takes  equal  and 
opposite  values  at  opposite  points  on  the  surfaces  of  the  plate.  If  the  wave 
encounters  a  small,  isolated  aperture,  this  pressure  differential  causes  fluid 
to  flow  through  the  aperture,  and  part  of  the  incident  energy  is  converted 
into  the  kinetic  energy  of  the  vorticity  produced  in  the  aperture.  An 
approximate  expression  for  the  wave  energy  absorbed  in  this  way  is  easily 
derived  as  follows. 


Let  the  bending  wave  propagate  in  the  x- direction  and  denote  the 
displacement  of  the  plate  by 

r  =  k,  >  0,  (4.1) 

where  is  constant  and  kh  ^  1  (wavelength  P  h) .  It  follows  from  equations 
(2.8)  and  (2.9)  that  the  pressure  jump  [p]  =  p(x,y,+0,t)  -  p(x,y,-0,t)  is 


The  bending  wavenumber  k  is  the  real  positive  zero  of  the  dispersion  equation 
obtained  by  substituting  from  (4.1)  and  (4.2)  into  equation  (2.7)  with  a  =  0 
and  set  equal  to  C,  i.e.,  of 

D(k,w)  =  Bk'’  -  mw^  -  ^  =  0.  (4.3) 

Q 

This  zero  satisfies  k  >  k  ,  and  characterizes  a  wave  that  propagates 

O 

subsonically  relative  to  the  fluid,  and  whose  influence  decays  exponentially 
with  distance  from  the  plate. 

Suppose  the  wave  impinges  on  a  small  circular  aperture  whose  center  is  at 
the  origin,  and  whose  radius  R  is  much  smaller  than  the  bending  wavelength  (kR 
I).  It  may  be  assumed,  without  loss  of  generality,  that  is  real,  so  that 
the  velocity  in  the  aperture,  ,  is  determined  by  equation  (2.4)  with  the 
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right  hand  side  replaced  by  (2cij^Zp/Vk2-/c2)cos (cot)  .  When  the  equation  is 
nondimensionalized  as  in  (3.7),  we  obtain 


R  dT 


47rZp 

o 


cos(27rT)  . 


(4.4) 


The  bending  wave  power  11^  dissipated  at  the  aperture  is  supplied  to  the 
aperture  flow  via  the  driving  pressure  [p] ,  and  is  easily  seen  to  be  equal  to 
-7rR^<V^  [p]>,  where  the  angle  brackets  denote  a  time  average,  i.e., 

Ha  =  -  |wR2J  VA(t)[p(t)]dt,  (4.5) 

where  the  integration  is  over  the  period  2%/u  of  the  wave .  The  energy  flux 
associated  with  the  incident  bending  wave  is  shared  between  that  propagated  by 
the  elastic  motions  in  the  plate  and  that  conveyed  in  the  evanescent  motions 
of  the  fluid.  The  power  flux  11,  say,  per  unit  length  of  wavefront  parallel  to 
the  plate  can  be  calculated  from  the  formula  II  =  (w| | ^/4)aD(k, w)/ak,  where  D 
is  the  dispersion  function  (4.3),  and  the  derivative  is  evaluated  at  the 
bending  wavenumber.  The  ratio  Ha/II  defines  the  absorption  cross-section  of 
the  aperture,  which  may  be  cast  in  the  following  nondimensional  form. 


_Ea 

2Rn 


k[5k'‘ -4(kK^  )2-iiiw2]  M 


L 


V,  (T)cos(27rT)dT. 


(4.6) 


where  M  =  wZ^/c^  is  the  amplitude  of  the  bending  wave  surface  velocity 
expressed  as  a  Mach  number. 

The  dependence  of  IIa/ZRII  on  frequency  predicted  by  this  formula  is  shown 
in  Figures  6  and  7  for  different  values  of  M  respectively  for  the  two  cases 
(i)  a  steel  plate  in  water  with  R/h  =  3,  and  (ii)  an  aluminium  plate  in  air 
with  R/h  =10.  In  these  figures  the  frequency  is  normalized  by  the 
coincidence  frequency  of  the  plate;  co^  =  c^(m/B)^^^  (above  which  bending  waves 
in  vacuo  have  phase  velocity  w/k  >  ) .  These  results  indicate  that,  however 

small  be  the  amplitude  M  of  the  incident  wave,  a  significant  fraction  of  the 
wave  energy  is  absorbed  at  sufficiently  small  frequencies,  where  the  curves 
peak. 


114 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


Equation  (4.6)  can  be  used  to  derive  an  approximate  formula  for  the 
damping  of  bending  waves  on  a  homogeneously  perforated  plate.  If  there  are  N 
apertures  per  unit  area  of  the  plate,  the  change  dH  in  the  bending  wave  power 
over  a  distance  6x  in  the  direction  of  propagation  of  the  wave  is 
approximately  equal  to  -2NRIIS5x,  i.e., 


■1  ”  -  2MI:  = 

n  ax 


k[5k^  -4(k/c^  )2-nictj2]  M 


I. 


V.  (T)cos(27rT)dT. 


(4.7) 


To  account  for  damping  the  amplitude  of  the  plate  motion  must  now  decrease 
with  increasing  x.  If  we  define 


-1  aZp  _  -1  an 

Zpdx  ~  m  ’ 


(4.8) 


then  fp  =  ,  i.e.,  Z^  =  Z^e'^i^,  where  Zp  =  at  x  =  0. 

(>0)  may  be  regarded  as  the  imaginary  component  of  a  complex  wavenumber  k  = 
k+ikj ,  where,  however,  k^  must  itself  depend  on  x,  because  of  the  nonlinear 
nature  of  the  damping. 


This  formula  for  k^  neglects  the  collective  influence  of  the  apertures  on 
the  wavenumber  k,  i.e.,  it  assumes  that  the  bending  wave  dispersion  equation 
is  (4.3),  as  for  an  unperforated  plate.  Similarly,  it  ignores  the  important 
fact  that  the  amplitude  of  the  motion  in  the  apertures  will  tend  to  be  smaller 
than  for  an  isolated  aperture,  because  the  pressure  difference  [p]  must  be 
very  much  smaller  for  the  perforated  plate,  even  if  the  amplitude  of  the  plate 
motion  is  unchanged.  Thus,  equation  (4.7)  would  be  expected  to  determine  an 
upper  bound  for  the  actual  damping  due  to  vorticity  production. 

It  is  convenient  to  express  the  attenuation  in  terms  of  the  fractional 
power  dissipated  per  wavelength  (2ff/k)  of  propagation,  which  is  equal  to 
407rlog^g(e)kj/k  (dB)  =  5.46kj/k  (dB)  .  Predictions  of  this  derived  from  (4.7), 
(4.8)  are  illustrated  in  Figure  8  for  a  thin  aluminium  screen  in  air  with  3% 
fractional  open  area  and  aperture  radius  R  =  lOh,  when  the  vibration  Mach 
number  M  =  10'^,  lO"'*  and  10"^.  These  results,  and  the  physical  constraints 
on  the  amplitudes  of  possible  plate  motions  in  practice,  suggest  that  the 
damping  is  significant  only  at  small  frequencies  -  M  (when  the  aperture 
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motion  is  nonlinear  even  when  M  0,  because  the  lengths  of  the  jets  become 
large  when  the  period  2i^/o)  . 

4.2  Equations  for  bending  waves  on  a  perforated  plate 

To  obtain  more  precise  estimates  of  the  damping,  and  extend  these  results 
to  underwater  applications  where  the  pressure  differential  [p]  produced  by 
surface  motion  might  be  expected  to  be  very  much  larger  than  in  air,  we  now 
proceed  to  the  formulate  a  self-consistent  analytical  model  of  flexural  wave 
propagation  over  a  homogeneously  perforated  plate.  Consider  motions  in  which 
the  displacement  of  the  plate  has  the  harmonic,  traveling  wave  form 

fp  =  0,  >  0,  (4.9) 

where  may  be  taken  to  be  a  real,  positive  constant.  The  wavenumber  k  must 
be  complex,  however,  and  for  propagation  in  the  positive  x-direction  we  write, 

K  =  k  +  ikj  ,  (4.10) 

where  both  k  and  k^  must  both  be  positive. 

Expand  the  mean  aperture  displacement  ^ ^  in  the  form 

=  Re[z^ei('‘''-“^)  +  I  j  ,  (4.11) 

^  n>2 

where  the  coefficients  Z"  (n  >  1)  are  generally  complex  quantities.  The  first 
of  these,  Z^  is  assumed  to  be  constant,  so  that  the  decay  of  this  term  as  the 
wave  propagates  is  accounted  for  by  the  imaginary  part  of  the  wavenumber  /c . 
This  hypothesis,  and  the  related  assumption  that  is  constant,  is  correct  in 
a  first  approximation;  it  will  be  seen,  however,  that  and  Z^  must  actually 
vary  slowly  with  x,  but  only  over  distances  large  compared  to  the  dissipation 
length  1/kj .  The  higher  order  terms  in  (4.11)  have  real  wavenumber  k,  and 
their  decay  is  accounted  for  by  assuming  the  coefficients  Z^  (n  >  2)  to  be 
slowly  varying  functions  of  x  on  a  scale  of  the  bending  wavelength  27r/k. 

The  pressure  fluctuations  on  either  side  of  the  plate  satisfy  the  linear 
wave  equation  (2.8) .  The  solution  of  this  equation  must  have  outgoing 
behavior  with  respect  to  the  plate,  and  may  be  expressed  in  terms  of  and 
via  the  boundary  conditions  (2.9)  where,  as  before,  the  median  section  of  the 


116 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


plate  is  assumed  to  coincide  with  the  plane  z  =  0  in  the  undisturbed  state. 

At  distances  from  the  plate  that  are  large  compared  to  the  aperture  spacing  d, 
we  find 


P  = 


-sgn(z)p^w^Re 


O 


] 


»ssn(z)p„u,^  o-f  7  „zn^n(l(kx-»t)-|2|VP^)l  (412) 

Vk2-/c^  L  A  j 

o  n>2 


The  first  term  on  the  right  hand  side  is  the  bending  wave  pressure  field. 
The  sum  represents  the  contribution  from  the  nonlinear  components  of  the 
aperture  motions.  We  shall  argue  that  this  latter  contribution  is  negligible 
at  distances  much  larger  than  d  from  the  screen,  because  "waves"  of 
wavenumbers  nk  and  frequencies  nw  (n  >  2)  induced  by  the  aperture  fluxes  do 
not  satisfy  the  dispersion  equation  (analogous  to  (4.3)  and  to  be  derived 
below)  for  flexural  waves  on  the  perforated  plate,  and  their  growth  to  finite 
amplitude  is  therefore  opposed  by  the  plate.  Thus,  the  pressure  jump  that 
drives  the  aperture  flows  is  taken  to  be 


[p] 


-2p  w^Re 

^  O 


o 


(4.13) 


Further  justification  of  this  approximation  is  given  in  §4.4,  where  an 
alternative  approach  is  outlined  that  is  strictly  applicable  when  the  damping 
due  to  vorticity  production  is  small. 


The  equation  describing  the  unsteady  motion  in  an  aperture  whose  center 
is  at  (x,y)  is  deduced  from  (2.4)  by  replacing  t  by  t  -  kx/w,  and  [p]  by 
(4.13).  In  nondimens ional  form  we  find 


i 

R  aT 


+  ^2^aI^A 


47rRe 


-[(l-a)Zp/R  +  aZ^/R]e 
.  RV«;2  -/c2 


•27riT 


]■ 


(4.14) 


where  T  =  (cot-kx)/2;r ,  and 


Z'e-ih*. 

A 


(4.15) 
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The  solution  of  equation  (4.14)  depends  on  the  complex  wavenumber  k, 
which  is  the  solution  of  the  appropriate  dispersion  equation  for  the 
perforated  plate.  The  latter  is  obtained  by  substituting  from  (4.9),  (4.11) 
and  (4.12)  into  the  linear,  generalized  bending  wave  equation  (2.7).  When 
attention  is  confined  that  component  of  the  motion  of  wavenumber  k  and  radian 
frequency  w  the  nonlinear  terms  from  (4.11)  and  (4.12)  may  be  discarded,  and 
(2.7)  yields 


fl  -  f 

2py  r. 

L  + 

(l-ct)m  ZpJ 

o 

1-a  ZpJ 

Equations  (4.14)  and  (4.16)  must  solved  simultaneously  for  and  the 

wavenumber  k.  The  definition  (4.11)  supplies 


(4.17) 


4.3  Numerical  results 

The  nonlinear  aperture  equation  (4.14)  is  solved  in  conjunction  with 
equation  (3.8c),  which  determines  the  variation  of  the  length  i(T)  via  the 
empirical  relation  (2.6).  The  solution,  which  must  be  periodic  with  period 
2n,  depends  on  both  the  local  amplitude  of  the  plate  motions  and  the 
wavenumber  k,  given  by  (4.16). 

The  solution  can  be  effected  by  the  following  iterative  procedure.  The 
frequency  w  and  magnitude  Z^  of  the  local  amplitude  of  the  plate  motion  are 
prescribed,  and  equation  (4.16)  is  first  solved  for  k  by  setting  Z^/Z^  =  0. 
This  initial  estimate  for  k  will  be  real,  and  corresponds  to  one  of  the  two 
equal  and  opposite  real  roots  of  the  dispersion  equation  (4.16)  when  Z^  =  0; 
the  positive  root  is  selected  when  the  wave  is  taken  to  propagate  in  the 
positive  x-direction.  The  presence  of  the  apertures  will  change  this  root, 
and  endow  it  with  a  positive  imaginary  part  k^.  that  accounts  for  the  damping 
of  the  wave.  Next,  (i)  equations  (4.14)  and  (3.8c)  are  solved,  and  the  result 
used  to  calculate  a  corrected  value  of  Z^/Z^  =  (4.17);  (ii)  the 

revised  value  of  Z^/Zp  is  inserted  into  (4.16)  and  the  wavenumber  k  is 
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re-calculated  by  iterating  about  its  value  for  Z^/Z^  —  0;  steps  (i)  and  (ii) 
are  repeated  until  convergence  is  achieved. 


If  the  local  plate  amplitude  is  expressed  as  a  Mach  number  M  =  wZ^/c^, 
this  naive  iteration  procedure  is  found  to  converge  rapidly  except  at  very 
small  values  of  M  (typically  less  than  10"^)  and  within  an  intermediate  range 
of  frequencies  w/u>  marginally  larger  than  M,  wherein  successive  iterates  of 
Z^/Zp  tend  to  fluctuate  wildly.  This  difficulty  is  avoided  by  using  at  stage 
(ii)  of  the  iteration  the  running  mean  <Z^/Zp>  instead  of  the  most  recently 
calculated  approximation  for  ,  where  at  the  nth  iteration 


<zi/Zp)  = 


I  (Z;'/Zp).,  n  >  1. 


When  n  becomes  large  successive  estimates  of  the  wavenumber  k  obtained  in 
stage  (i)  are  then  found  to  converge  mono tonic ally,  and  a  rapidly  converging 
estimate  for  k  can  be  obtained  by  making  use  of  a  simple  extrapolation  formula 
of  the  type 


ln[ (n-l)/(n-2) ]/ln[n/(n-l) ] . 


The  fractional  bending  wave  power  dissipated  per  wavelength  (27r/k)  of 
propagation  (=  54.6kp/k  dB)  calculated  by  this  procedure  is  illustrated  in 
Figure  9  for  the  aluminium  plate  in  air  considered  in  §4.1.  A  comparison  of 
these  predictions  with  those  shown  in  Figure  8  (based  on  the  absorption 
cross-section  of  an  isolated  aperture  in  an  unperforated  plate)  reveals  that 
the  major  differences  occur  at  lower  plate  Mach  numbers  M.  The  scattering 
cross-section  theory  tends  to  overestimate  the  damping  except  at  very  low 
frequencies  where  the  attenuation  is  small.  Predictions  of  equations  (4.14)  - 
(4.17)  for  a  perforated  steel  plate  in  water  (for  a  =  0.05  and  R/h  =  3)  are 
given  in  Figure  10.  Large  fluctuations  in  pressure  can  now  be  caused  by 
flexural  motions  of  the  plate  because  of  the  relatively  high  density  of  water, 
and  the  predicted  attenuations  are  correspondingly  much  higher  than  for  a 
screen  in  air. 


119 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


To  estimate  the  significance  of  these  predictions,  observe  that  for  a 
plate  of  thickness  h,  the  amplitude  of  the  plate  motion  can  be  expressed  in 
the  form 


Mh 

'P  (w/cd^  )  (w^h/c^  )  ’ 


(4.18) 


where  w^h/c^  =0.22  and  0.95  respectively  for  an  aluminium  plate  in  air  and  a 
steel  plate  in  water.  Thus,  must  typically  be  larger  than  M  to  ensure 

that  the  amplitude  of  the  plate  motion  is  small  compared  to  the  plate 
thickness.  Although  the  latter  condition  is  not  strictly  necessary  for  the 
applicability  of  the  thin  plate  equation  at  low  wavenumbers,  structural 
motions  are  in  practice  unlikely  to  attain  amplitudes  comparable  to  h. 
Accordingly,  it  is  only  at  frequencies  >  M,  to  the  right  of  the 

corresponding  peaks  in  the  Figures  9  and  10,  that  the  present  predictions  are 
relevant.  However,  surface  Mach  numbers  M  are  generally  much  smaller  than 
discussed  here,  and  correspondingly  at  the  finite  frequencies  of  interest  in 
applications  (typically  >  10’'’),  it  appears  from  these  results  that 

nonlinear  damping  will  be  negligible. 

4.4  Alternative  analysis  applicable  for  small  damping 

The  aperture  equation  (4.14)  was  derived  by  assuming  that  the  dominant 
contribution  to  the  local  pressure  jump  [p]  is  given  by  the  linear  theory 
approximation  (4.13),  associated  with  the  bending  wave  component  of  the 
pressure.  This  approximation  accords  with  intuition,  but  its  accuracy  is 
difficult  to  quantify  in  general.  An  alternative  formulation  is  possible  when 
the  damping  is  small,  however,  that  fully  justifies  the  approximation  in  that 
limit . 

Let  the  plate  and  aperture  displacements  be  taken  in  the  form  (c.f. 

(4.9),  (4.11)) 

Cp  =  Re[Zpe^(’^-“^)],  =  Ref  ^  ,  o)  >  0,  (4.19) 

where  the  wavenumber  k  is  required  to  be  real  and  positive  (for  propagation  in 
the  positive  x-direction) ,  and  where  all  of  the  coefficients  and  (n  >  1) 
are  assumed  to  be  slowly  varying  functions  of  x  on  a  scale  of  the  flexural 
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wavelength  2n/k.  The  damping  of  the  wave  is  now  accounted  for  by  the  decay  of 
the  coefficients  Z„  and  .  As  before,  it  may  be  assumed  that  is  real. 

P  A 

The  pressure  in  the  fluid  at  distances  >  d  from  the  plate  is 


-sgn(z)PoU)^ 

Vk^-K^ 

o 


(1 


-a)Z„cos(kx-a)t)e' +  p 


(4.20) 


where  the  first  term  on  the  right  is  the  contribution  from  the  plate  motion, 
and  p  is  the  collective  effect  of  the  apertures.  When  the  damping  is  small 
p^  may  be  regarded  locally  as  a  function  of  kx-wt  and  z,  in  which  case  it 
satisfies  the  following  reduced  form  of  equation  (2.8) 


3^ 

5x2 


=  0. 


(4.21) 


Now  K^/k^  <  1  for  a  subsonically  propagating  flexural  wave,  and  the 

O  ' 

solution  of  this  equation  with  outgoing  wave  behavior  can  therefore  be 
expressed  in  terms  of  the  limiting  values  as  z  -»  ±0  of  5p^/5z  s  -ap^5V^/5t  by 
the  methods  of  two-dimensional  potential  theory  (Carrier,  Krook  &  Pearson 
1966).  Correct  to  the  present  order  of  approximation  we  can  write  5V^/5t  = 

- (k/w)5V^/5x,  and  it  then  follows  that 


sgn(z)a/3QW  r”  Va  (^H  ^  ~  (t-kx/tj)  ]  dr _ 

7rVk2-/c2  { [r- (t-kx/w)  ]  2  +  z2 (k2/w2 ,i/c2)  } 


(4.22) 


The  equation  of  motion  of  fluid  in  the  aperture  at  (x,y)  is  formed  by 
replacing  t  by  t  -  kx/w  in  equation  (2.4),  and  calculating  [p]  from  (4.20)  and 
(4.22).  In  the  usual  nondimens ional  form  the  result  is 


£  V  IV  I  .  I  ^a(T’.)^ 

R  ST  ,,2  A  I  a'  RVk2-K2  T’  -T 


47r(l-a)  Zp 


cos  (ZttT) 


(4.23) 


where  T  =  (wt-kx)/27r.  The  principal  value  of  the  integral  represents  the 
pressure  jump  associated  with  the  aperture  flows. 

The  solution  of  this  equation  determines  from  the  relation  d^^/dt  = 

V, ,  from  which  we  find 
A 
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Z 


1 

A 


27riT 


dT. 


(4.24) 


Next  substitute  from  (4.19)  and  (4.20)  into  the  generalized  bending  wave 
equation  (2.7).  In  this  linear  equation  we  are  concerned  solely  with  that 
component  of  the  solution  of  wavenumber  k  and  frequency  w,  but  with  account 
taken  of  the  slow  variation  in  the  amplitudes  caused  by  damping. 

Hence,  we  can  set 

[pl  -^^Ra([(l-o)Z,  (4.25) 

o 


When  these  various  substitutions  are  made  in  equation  (2.7)  and  the 
separate  coefficients  of  cos(27rT)  and  sin(27rT)  in  the  resulting  equation  are 
equated  to  zero,  we  find 


fi  -  f 

^Imw^ 

[ 

1  +  1 

(l-a)m 

Z  J 

P 

0 

(l-a)ZpJ 

)  =0, 


(4.26) 


-1  aZp 
ZpSx 


2a(Zl./Z^)pyV^^[l 


4/c2)Bk^  - 

O  ^ 


RVk^-K^ ] 

O 


(l-a)m 


(4.27) 


where  Z^  ,  Z^^  respectively  denote  the  real  and  imaginary  parts  of  Z^ .  In 
deriving  these  equations  derivatives  of  Z^  of  higher  order  than  the  first  have 
been  neglected,  and  we  have  used  the  result  (l/Zp)9Zp/3x  =  (1/Z^  )9Z^/9x 

When  the  ratio  Z^/Z^  is  known,  equation  (4.26)  becomes  the  dispersion 
equation  that  determines  the  real  wavenumber  k  in  terms  of  w.  For  this  value 
of  k  equation  (4.27)  defines  an  imaginary  wavenumber  component  k^.  = 

- (l/Zp)9Zp/3x,  in  terms  of  which  the  fractional  wave  power  absorbed  per 
wavelength  is  54.6kp/k  dB,  as  before. 


Equations  (4.23)  -  (4.27)  are  to  be  solved  simultaneously  with  equation 
(3.8c),  and  this  is  accomplished  by  the  iterative  procedure  described  in  §4.3. 
The  details  are  simplified  by  observing  that,  when  \(T)  has  period  1,  the 
principal  value  integral  of  (4.23)  can  be  transformed  to  a  regular  integral 
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over  (0,1),  and  (4.23)  cast  in  the  form 

+  lv  \V  \  +  r'  <^a(T’)  -  VA(T-[T])}dT’^  hna-a)Z^  cos(2^T)  , 

R  3T  ^2  A'  A  I  RVk2-«;2  Jo  tan{7r(T’ -  (T- [T] ) ) }  RWk^-K^ 

(4.28) 

where  [T]  is  the  integer  part  of  T.  The  periodic  solution  of  this  equation  is 
obtained  by  integrating  from  T  =  0,  with  V^(0)  =  0.  During  the  first  period  0 
<  T  <  1  the  contribution  from  the  integral  is  ignored.  Integration  of  the 
equation  over  subsequent  periodic  intervals  (n,  n+1)  is  performed  by 
evaluating  the  integral  from  the  solution  in  (n-l,n).  This  procedure 
typically  yields  a  closely  periodic  solution  for  n  >  6. 

The  calculations  have  been  performed  for  two  of  the  cases  discussed  in 
§4.3  involving  an  aluminium  plate  in  air  and  a  steel  plate  in  water,  both  for 
a  surface  Mach  number  M  =  wZ^/c^  =  10  * .  The  predicted  attenuations  are 
indicated  by  the  dotted  curves  in  Figures  9  and  10,  and  are  seen  to  conform 
closely  to  the  corresponding  predictions  based  on  the  approximation  (4.13). 
This  agreement  occurs  only  for  the  small  attenuations  obtaining  in  these 
cases,  however.  Indeed,  equations  (4.26)  and  (4.27),  which  determine  k  and  k^ 
=  - (1/Zp )aZp/ax,  may  also  be  deduced  from  the  general  dispersion  equation 
(4.16),  by  expanding  the  solution  «  =  k  +  ik^  in  powers  of  Z^^ /Z^ .  The  zeroth 
and  first  order  approximations  are  respectively  equivalent  to  equations  (4.26) 
and  (4.27).  Equations  (4.26)  and  (4.27),  with  terms  of  0(q:^)  discarded,  may 
also  be  seen  to  reduce  to  the  scattering  cross-section  approximation  (4.3), 
(4.7)  (wherein  the  integral  determines  Z^.  /2R) . 
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5.  CONCLUSION 

In  this  paper  we  have  investigated  the  influence  of  flow  nonlinearity  in 
the  apertures  of  a  perforated  elastic  screen  on  the  attenuation  of  sound  and 
of  flexural  vibrations.  The  damping  is  a  result  of  the  direct  transfer  of 
acoustic  or  vibrational  energy  to  the  kinetic  energy  of  vorticity  generated  in 
the  apertures  by  the  incident  disturbance.  The  sound  or  flexural  wave 
produces  a  pressure  differential  across  the  plate,  causing  fluid  to  "jet" 
through  the  apertures.  Large  amplitude  sound  waves  can  experience  substantial 
attenuations  which,  however,  progressively  diminish  at  lower  frequencies,  when 
an  elastic  plate  ultimately  becomes  transparent  to  incident  sound.  Similarly, 
it  is  predicted  that  flexural  waves  can  be  significantly  damped,  by  several  dB 
per  wavelength  of  propagation.  Aperture  nonlinearity  is  particularly 
important  for  low  amplitude  bending  waves  at  low  frequencies  =  M  (M  being 

the  Mach  number  of  the  plate  surface  velocity) ,  when  the  wave  period  is 
sufficiently  long  that  nonlinear  processes  within  an  aperture  have  ample  time 
acquire  a  relatively  large  amplitude.  However,  the  frequencies  at  which  this 
occurs  are  likely  to  be  below  those  normally  of  interest  in  practice.  At 
higher  frequencies,  and  at  surface  wave  amplitudes  encountered  in 
applications,  the  predicted  structural  damping  is  usually  small,  and  it  is 
certainly  negligible  compared  with  that  which  can  be  achieved  (by  the  same 
vorticity  production  mechanism)  when  there  exists  a  grazing  mean  flow  past  the 
screen  or  bias  flow  through  the  apertures. 
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PRINCIPAL  SYMBOLS 

bending  stiffness  of  the  plate 
speed  of  sound  in  the  fluid 
dispersion  function  (4.3) 
plate  thickness 

real  part  of  flexural  wavenumber 
imaginary  part  of  flexural  wavenumber 
Rayleigh  conductivity 

aperture  flow  length  defined  by  (2.6) 

effective  jet  length 

p^h,  mass  per  unit  area  of  the  plate 

u)Z  /c  or  wZ  /c  ,  Mach  number  of  the  plate  displacement  velocity 

number  of  apertures  per  unit  area  of  the  plate 

pressure 

p(x,y ,+0) -p(x,y , -0) 
aperture  radius 
time 

cot/ZTT  or  (a;t-kx)/27r 
aggregate  plate  velocity,  (2.10) 
mean  aperture  velocity 
jet  velocity 

amplitude  of  the  nth  Fourier  coefficient  of  the  aperture 
motion 

amplitude  of  plate  displacement 


N;rR^  ,  fractional  open  area  of  the  plate 
absorption  coefficient  (3.6) 
contraction  ratio  of  jet 
(l-a)fp+af^ 

mean  displacement  of  fluid  in  an  aperture 
displacement  of  the  plate 
k+ikj ,  complex  wavenumber 


127 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


Po 
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co/c  ,  acoustic  wavenumber 

velocity  potential 

end  correction 

mean  fluid  density 

density  of  plate  material 

Poisson’s  ratio  of  the  plate  material 

scattering  cross-section  (4.6) 

acoustic  or  bending  wave  power 

radian  frequency 

c^(m/B)^^^,  coincidence  frequency 
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screen 

PjCos{w(t-z/c^)} 

PR(t  +  z/C^) 

-^A/XA/ 

z  =  0 


Figure  2.  Interaction  of  a  normally  incident  sound  wave  with  a  perforated, 
elastic  screen. 
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Aluminium/ Air,  R/h  =  10,  k^R  =  0.05,  a  =  0.03 


Figure  3.  Relative  contributions  of  the  apertures  and  plate  ((l-Q)Vp) 

to  the  effective  normal  velocity  V  of  the  perforated  plate. 
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Air/ Aluminium,  R/h  =  10 


Figure  7.  Absorption  cross  section  calculated  from  equation  (4.6)  for  a 
bending  wave  on  an  aluminium  plate  in  air  incident  on  a  small 
circular  aperture. 
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Estimate  based  on  absorption  cross-section 


Figure  8.  Estimates  of  the  damping  of  bending  waves  on  a  perforated  plate  in 
air  based  on  the  absorption  cross-section  formulae  (4.6)  -  (4.8). 
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Figure  9.  Damping  of  bending  waves  on  a  perforated  aluminium  plate  in  air 

for  a  =  0.03  and  R/h  =  10:  - ,  predictions  of  equations 

(4.14)  -  (4.17):  ••••,  small  damping  approximation  based  on 
equation  (4.23). 


Absorption  per  wavelength  (dB) 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


Water/Steel,  R/h  =  3,  a  =  0.05 


Figure  10.  Damping  of  bending  waves  on  a  perforated  steel  plate  in  water  for 

a  =  0.05  and  R/h  =  3:  - ,  predictions  of  equations 

(4.14)  -  (4.17);  •••• ,  small  damping  approximation  based  on 
equation  (4.23). 
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RECOMMENDATIONS  FOR  FURTHER  RESEARCH 
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1 .  INTRODUCTION 

The  B-IB  and  F15  aircraft  are  configured  with  twin  engine  nacelles,  and 
both  have  suffered  premature  failure  of  external  nozzle  engine  flaps  caused  by 
high  dynamic  pressures.  Model  tests  [1]  confirm  that  the  highest  dynamic 
pressures  occur  where  structural  damage  is  observed  on  full  scale  aircraft, 
principally  between  the  nozzles  of  the  twin  nacelle.  Seiner  et  al  [2] 
correlated  these  high  pressures  with  jet  screech  tones  produced  by  the 
interaction  of  turbulence  and  shock  waves  in  the  supersonic  jet  plumes.  In 
particular,  intense  surface  pressures  corresponded  to  a  dynamic  coupling  of 
the  neighboring  plumes,  a  common  feature  of  parallel  jets  whose  separation  is 
less  than  about  four  or  five  jet  diameters  [3  -  5].  More  recent  wind  tunnel 
tests  [6,  7]  have  revealed  that  the  supersonic  plume  resonance  is  actually 
important  only  at  low  flight  Mach  numbers,  typically  less  than  about  0.5.  At 
higher  flight  speeds  the  pressures  appear  to  be  attributable  to  large-scale 
vortex  structures  impinging  from  the  aircraft  forebody. 

In  practice  the  dominant  mechanism  is  determined  by  flight  profile  and 
aircraft  configuration.  The  jet  plume  resonance  can  be  eliminated  by 
inserting  tabs  into  the  nozzle  flow  or,  more  effectively,  by  means  of  a  small 
supersonic  jet  tube  within  the  nozzle  [7].  However,  there  is  currently  no 
effective  means  of  controlling  the  vortex  dominated  pressures. 

In  Chapter  2  we  have  advanced  and  developed  arguments  to  determine  the 
damping  of  structural  vibrations  by  vorticity  production.  The  idea  is  that, 
in  the  presence  of  a  mean  flow  a  vibrating  structure  will  generate  vorticity 
that  is  subsequently  convected  away  by  the  flow.  Energy  is  required  to  create 
vorticity  so  that,  provided  the  back  reaction  of  the  pressure  field  of  the 
vorticity  on  the  structure  is  small,  structural  vibrations  tend  to  be  damped. 
Strong  back  reactions  are  avoided  if  the  length  scale  of  the  vorticity  is 
small,  and  this  can  be  achieved  by  perforating  all  or  part  of  the  vibrating 
structure,  and  forcing  the  mean  flow  to  pass  either  through  or  over  the 
perforates,  in  which  case  the  vortex  length  scale  is  comparable  to  the 
perforation  scale. 

To  assess  the  feasibility  of  this  method  for  inhibiting  flap  vibration  we 
have  performed  a  series  of  basic,  proof  of  principle  tests  using  vibrating 
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perforated  plates  in  the  Harvard  University  flume.  The  flume  has  not  been 
used  for  some  time,  and  is  in  need  of  routine  maintenance.  In  spite  of  this 
we  have  obtained  encouraging  results  which  are  discussed  below  in  §3.  In  §2 
we  outline  an  approximate  theory  of  vibration  damping  based  on  our  work 
described  in  [8,  9]. 

2.  THEORY  OF  VIBRATION  DAMPING 

Consider  a  flat,  rigid,  rectangular  airfoil  of  large  aspect  ratio  and 
chord  2a  set  at  zero  angle  of  attack  to  a  low  Mach  number  mean  flow.  The  flow 
is  at  speed  U  in  the  x  direction  of  the  rectangular  coordinate  system  (x,y,z). 
The  airfoil  executes  small  amplitude,  time  harmonic  motions  normal  to  its 
su^rface  (in  the  y  direction,  see  Figure  1) ,  such  that  its  displacement  at  time 
t  is  5-Qexp(-iwt)  ,  where  the  real  parts  of  all  complex  quantities  are  to  be 
taken.  The  rate  11^  at  which  energy  is  supplied  to  the  flow  due  to  vortex 
shedding  at  the  trailing  edge  of  the  airfoil  is  given  approximately  by 

=  TrpaUw^CjjCwa/U)  Ko  1^  (1) 

per  unit  span,  wliere  p  is  tine  mean  density  of  the  fluid  and  is  the  real 
part  of  the  Theodor sen  function 

C(x)  =  iH[^\x)/{H^^\x)+iH^^\x)}.  (2) 

When  coa/U  >  1,  (x)  ~  0.5,  and 

Hq  -  TrpaUw^  ICq  I V2,  >  1.  (3) 

Consider  next  the  case  of  an  airfoil  perforated  with  a  uniform 
distribution  of  small  circular  apertures  of  radius  R  ^  a.  The  oscillatory 
motion  forces  fluid  through  the  apertures  and  leads  to  the  ejection  of 
additional  vorticity  into  the  flow.  Provided  the  aperture  Strouhal 
number  wR/U  is  less  than  about  3,  theoretical  estimates  [9,  10]  indicate  that 
energy  is  supplied  to  the  fluid  at  a  rate 

n;  =  pwA|  [(p]  1V2. 

where  [cp]  =  is  the  difference  in  the  unsteady  velocity  potentials  just 

above  and  below  the  aperture,  and  A  =  A(a;R/U)  is  the  imaginary  part  of  the 
Rayleigh  conductivity  of  the  aperture,  defined  by 
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=  2R(r-iA).  (5) 

In  the  absence  of  the  mean  flow,  is  real  and  equal  to  2R  (i .  e .  ,  T  =  1 ,  A  = 
0).  The  influence  of  tangential  flow  has  been  investigated  in  [9],  with 
predictions  for  T  and  A  shown  in  Figure  2.  In  the  low  frequency  range  0  < 
wR/U  <3,  A  is  positive  and  large,  and  airfoil  vibrations  will  tend  to  be 
damped  at  these  frequencies  by  vorticity  production  in  the  apertures. 

To  estimate  the  damping  when  the  distribution  of  apertures  over  the 
airfoil  is  uniform,  we  take  the  unsteady  potential  difference  across  the 
airfoil  to  be 

[cp]  =  iwf  ^  (a^ -x^  )^^^exp( -iwt)  ,  (6) 

where  is  the  amplitude  of  the  airfoil  motion,  /S’  is  a  constant,  and  the 
coordinate  origin  is  at  the  center  of  the  airfoil.  This  should  be  an  adequate 
approximation  for  large  values  of  wa/U,  except  for  points  very  close  to  the 
trailing  edge.  For  a  rigid  airfoil  /S’  =2,  but  its  value  is  actually  much 
smaller  because  of  the  reduction  in  [cp]  caused  by  "short  circuiting"  by  the 
apertures.  If  there  are  N  apertures  per  unit  area  of  the  airfoil,  and  a  = 

NttR^  is  the  fractional  open  area,  we  now  find  by  integration  of  (4)  over  the 
airfoil  chord  that  the  power  dissipated  per  unit  span  can  be  cast  in  the  form 

=  a/3pa)^a®A|f l^/R,  (7) 

where  /3  is  another  constant. 

Combining  (1)  and  (7),  for  an  airfoil  vibrating  with  amplitude  the 

net  power  dissipated  per  unit  span  by  the  vibrating  perforated  airfoil  becomes 

=  pw^aV[itC^  +  a/S(wa/U)(a/R)A]  (8) 

Suppose  next  that  the  11^  =  ,  i.e.,  that  the  vibratory  motions  of  the 

perforated  airfoil  at  amplitude  and  the  unperforated  airfoil  at  amplitude 
fp  are  maintained  by  driving  mechanisms  delivering  the  same  power.  Then 
equations  (1)  and  (8)  imply  that, 

Ifl/Col"  =  +  a;0(a,a/U)(a/R)A].  (9) 
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3.  PRELIMINARY  DAMPING  TESTS 

The  experiments  were  conducted  in  a  low  speed  water  flume  that  was 
capable  of  flow  speeds  up  to  one  meter  per  second.  Two  interchangeable, 
rectangular  steel  plates  were  used,  each  of  11"  span,  6"  chord  and  1/32" 
■thickness .  One  plate  was  unperforated  and  the  other  was  drilled  with  a 
uniform  distribution  of  circular  holes  of  diameter  0.25"  with  a  fractional 
open  area  a  =  0.1.  Each  plate  could  be  bolted  at  one  end  to  the  lower  end  of 
a  wooden  support,  as  indicated  schematically  in  Figure  3a  (where  the  mean  flow 
is  directed  into  the  paper) .  The  support  could  be  raised  or  lowered  to 
submerge  the  plate  to  a  water  depth  of  about  3  inches .  The  other  end  of  the 
plate  supported  a  3"  tall  aluminum  block,  on  the  top  of  which  was  mounted  a 
lead-eccentric  electric  motor  and  a  polyvinylidene  difluoride  (PVDF)  vibration 
sensor.  The  output  of  the  PVDF  was  displayed  on  a  digital  oscilloscope  which 
gave  a  direct  reading  of  the  voltage.  Steady  operation  of  the  motor  at  a 
fixed  speed  caused  the  plate  to  vibrate  at  a  known  frequency,  which  was 
measured  stroboscopically .  The  connection  between  the  block  and  plate  was 
rigid  so  that  motion  detected  by  the  sensor  also  gave  the  amplitude  of  the 
plate  motion.  The  arrangement  was  such  that  both  the  motor  and  vibration 
sensor  were  above  the  water  surface. 

The  flow  speed  in  the  flume  was  adjusted  to  about  0.6  m/s,  and  was 
measured  by  observing  the  average  travel  times  of  several  straws  placed  on  the 
water  surface  over  a  distance  of  2  yards.  The  electric  motor  was  set  to  a 
predetermined  speed,  and  the  output  (voltage)  from  the  vibration  sensor  was 
displayed  on  the  oscilloscope.  The  sensor  output  is  proportional  to  the 
vibration  amplitude  of  the  plate,  and  steady  state  amplitudes  were  recorded 
at  different  frequencies  for  both  the  unperforated  (Tq)  perforated  (T^) 

plates.  This  permitted  a  comparison  to  be  made  of  the  vibration  levels  when 
the  plates  were  submerged  to  the  same  depth,  the  motor  operated  at  the  same 
frequency,  and  when  the  delivered  mechanical  powers  were  the  same.  However, 
the  electric  motor  would  function  in  a  controlled  manner  only  for  a  restricted 
range  of  frequencies,  and  the  measurements  were  therefore  limited  to  the  three 
Strouhal  numbers  indicated  in  Table  1  (which  correspond  to  frequencies  between 
10  -  30  Hz).  The  measurements  were  repeated  on  several  different  occasions  to 
guarantee  the  consistency  of  the  data. 
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wR/U 

iTo/rj" 

20xlog|ro/rj  (dB) 

0.35 

2.4 

3.8 

0.5 

1.23 

0.9 

0.95 

1.86 

2.7 

TABLE  1 :  Measured  Damping 

The  measured  damping  (column  3  of  Table  1)  corresponds  to  the  prediction 
of  equation  (9),  and  is  plotted  as  open  squares  in  Figure  3b.  To  apply 
equation  (9)  to  these  measurements  we  first  recognize  that  when  wR/U  =  0(1) 
the  Strouhal  number  wa/U  based  on  the  airfoil  semi-chord  is  very  large,  so 
that  IIq  may  be  approximated  by  formula  (3)  (i.e.,  Cj^  =  0.5).  Then  (9)  becomes 

Ki/fol"  =  V[i  +  5(wR/u)a],  (10) 

where  6  =  (16/3?r^  )a)3(a/R)^  .  A  proper  estimate  of  the  coefficient  ^  should 
strictly  be  made  by  solving  a  boundary  value  problem  to  determine  the 
effective  velocity  potential  of  the  fluid  motion  produced  by  the  vibrating 
perforated  airfoil.  The  solution  of  this  problem  requires  the  numerical 
integration  of  the  equations  of  motion.  Experience  with  similar  problems 
concerned  with  the  acoustics  of  perforated  plates  [11]  indicates  that  ^  is 
likely  to  be  small. 

According  to  Figure  2  the  additional  damping  attributable  to  the 
apertures  should  increase  over  the  frequency  range  of  the  experiments.  This, 
is  evidently  not  the  case,  at  least  for  the  first  measurement  at  wR/U  =  0.35, 
which  exhibits  the  largest  damping  of  nearly  4  dB.  Actually,  in  spite  of  the 
care  exercised  in  performing  these  measurements,  it  was  not  possible  to  obtain 
a  definite  stroboscopic  reading  of  the  drive  speed  for  the  unperforated  plate 
when  wR/U  =  0.35:  the  motor  was  on  the  point  of  stalling,  and  the  best  that 
could  be  done  was  to  use  the  speed  obtained  for  the  perforated  plate  (for 
which  stalling  did  not  occur) .  This  presumably  accounts  for  the  anomalous 
result  in  Figure  3b  at  wR/U  =0.35.  At  the  higher  frequencies  the  constant 
motor  speed  was  easily  maintained,  and  the  motor  did  not  stall.  If  the 
measurement  at  wR/U  =  0.35  is  disregarded,  the  coefficient  6  of  equation  (10) 
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can  be  adjusted  to  give  the  least  squares  best  fit  to  the  remaining  two  data 
points.  The  optimum  occurs  when  6  =  0.96,  and  the  solid  curve  in  Figure  3b 
represents  the  damping  predicted  by  (10)  for  this  value  of  6.  The  theoretical 
curve  peaks  at  coR/U  =  1.6,  where  the  damping  rises  to  a  maximum  of  about  4  dB. 

4.  FURTHER  WORK 

The  possibility  of  achieving  vibration  damping  of  3  dB  or  more  can  be 
critical  for  the  continued  safe  operation  of  fatigue  sensitive  structures. 

Our  preliminary  experiments  indicate  that  additional  dampings  in  excess  of 
this  may  be  attainable  by  suitably  perforating  a  vibrating  control  surface. 

The  damping  can  be  either  passive,  when  the  surface  is  placed  in  a  tangential 
flow,  or  active,  when  fluid  is  forced  through  the  apertures,  or  perhaps  more 
optimally  a  combination  of  the  two  methods.  We  recommend  performing  two  sets 
of  experiments  and  concurrent  analytical  and  numerical  work.  The  first  set  of 
experiments  should  be  performed  in  the  Harvard  University  flume,  a  (largely 
neglected)  facility  that  is  ideal  for  such  studies;  it  can  be  used  to  validate 
theoretical  models  and  test  new  ideas  suggested  by  the  theory.  A  first 
objective  should  be  to  obtain  a  complete  data  set  of  results,  from  a  carefully 
controlled  series  of  tests  with  perforated  and  unperforated  plates,  to  permit 
validation  of  the  theoretical  models  discussed  in  this  report. 

Second,  wind  tunnel  measurements  of  the  type  discussed  above  for  the 
flume  should  be  performed,  but  in  a  form  more  directly  related  to  the  nozzle 
flap  fatigue  problem  (§1).  One  such  test  is  illustrated  schematically  in 
Figure  4,  in  which  the  vortex  wake  from  an  upstream,  spanwise  oriented 
cylinder  excites  vibrations  in  both  conventional  and  perforated  airfoils. 

This  configuration  is  an  idealized  analog  of  the  excitation  of  an  aircraft 
engine  nozzle  flap  by  vorticity  swept  over  the  flap  from  the  aircraft 
forebody.  Measurements  will  be  made  of  vibration  spectra  for  airfoils  of 
various  sections .  A  comparison  of  corresponding  results  obtained  for  the 
flume  and  wind  tunnel  will  permit  the  dependence  of  damping  on  Mach  number  to 
be  assessed  and  compared  with  theory.  The  influence  of  compressibility  on 
damping  will  likely  be  important  at  higher  Mach  numbers,  and  this  comparison 
will  indicate  when  such  effects  need  to  be  included  in  the  aperture  flow 
modeling . 


145 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


Analytical  work  should  complement  and  facilitate  interpretation  of  the 
experiments,  and  permit  extrapolation  of  measured  results  to  full  scale.  This 
might  include 

•  a  complete  analytical  model  of  the  flume  experiments 

•  the  development  of  a  general  purpose  computer  code  for  the  complex 
j^ayleigh  conductivity  of  an  aperture  of  arbitrary  shape  in  a  mean 
grazing  flow. 

•  numerical  procedures  to  account  for  the  influence  of  finite 
airfoil/flap  thickness,  and  "straking"  of  the  flow-through  cylindrical 
neck  connecting  the  aperture  faces. 

•  an  investigation  of  the  combined  effects  of  grazing  and  bias  flows 
through  the  apertures . 
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Figure  1.  Airfoil  vibrating  in  uniform  mean  stream. 
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Figure  4.  Vibrations  induced  by  a  vortex  wake. 
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ABSTRACT 

A  summary  review  is  given  of  theoretical  predictions  of  the 
damping  of  sound  and  structural  vibrations  by  vorticity 
production  in  the  apertures  of  bias-  and  grazing-flow  perforated 
elastic  plates.  Unsteady  motion  causes  vorticity  to  be  generated 
in  the  perforates;  the  vorticity  and  its  energy  are  swept  away  by 
the  mean  flow  and  result  in  a  net  loss  of  acoustic  and  vibrational 
energy.  Generalized  forms  of  the  thin  plate  bending  wave 
equation  are  given  for  plates  perforated  with  small  circular 
apertures,  and  used  to  predict  the  attenuation  of  sound  and 
resonant  bending  waves.  Acoustic  damping  is  significant  when 
the  fluid  loading  is  small  enough  for  the  plate  to  be  regarded  as 
rigid.  Bending  waves  are  effectively  damped  only  when  the 
fluid  loading  is  large  and  there  is  a  substantial  pressure  drop 
across  the  plate;  when  this  occurs  the  predicted  attenuations  are 
comparable  with  those  usually  achieved  by  the  application  of 
elastomeric  damping  materials. 

INTRODUCTION 

Aerodynamic  sound  is  generated  by  turbulence  and  unsteady 
vorticity,  and  by  their  interactions  with  adjacent  structures. 
Acoustic  energy  can  also  be  dissipated  by  vorticity  production, 
however,  when  sound  impinges  on  a  solid  surface  or  fluid 
region  of  nonuniform  mean  density.  For  example,  sound 
incident  on  the  trailing  edge  of  an  airfoil  in  a  mean  flow  will 
cause  vorticity  to  be  shed  into  the  wake;  at  low  Mach  numbers 
this  leads  to  an  increase  in  the  kinetic  energy  of  the  mean  flow 
and  wake  at  the  expense  of  the  sound  [1,  2].  Similarly, 
perforated  plates  aligned  with  the  mean  flow  through  the  tube 
banks  of  a  heat  exchanger  cavity  can  be  used  to  suppress  cavity 
acoustic  resonances  by  the  production  of  vorticity  in  the 
perforates;  the  kinetic  energy  of  the  vorticity  is  extracted  from 
the  sound,  convected  away  by  the  flow,  and  ultimately 
dissipated  by  viscous  and  thermal  processes  [3,  4].  The  same 
mechanism  is  responsible  for  the  greatly  improved  attenuation  of 


"screech"  tones  in  the  jet  pipe  of  a  jet  engine  when  a  bias  flow 
is  maintained  through  the  perforated  heat  shield  normally  used 
to  protect  the  wall  liner  from  contact  with  hot  combustion 
products  [5,  6]. 

The  effectiveness  of  the  acoustic  attenuation  in  these  cases  is 
critically  dependent  on  the  presence  of  mean  flow  [7  -  10]. 
Analytical  predictions  for  rigid  perforated  screens  indicate  that 
the  attenuation  is  linearly  proportional  to  the  acoustic  amplitude, 
which  is  in  broad  agreement  with  experiment  [6  -  8,  11,  12].  In 
the  absence  of  flow  the  attenuation  is  either  a  nonlinear  function 
of  the  acoustic  amplitude  (and  therefore  weak)  or  is  dominated 
by  the  less  efficient  action  of  viscous  dissipation  in  the 
apertures. 

Theoretical  analyses  have  usually  assumed  the  screen,  plate, 
etc,  to  be  rigid,  although  in  practice  high  acoustic  intensities  are 
accompanied  by  structural  vibrations,  and  a  substantial  part  of 
the  noise  energy  is  contained  in  vibratory  modes  of  the 
structure.  These  modes  would  also  be  expected  to  be  damped 
by  vorticity  production,  and  a  primary  practical  objective  may 
actually  be  the  suppression  of  potentially  harmful  structural 
vibrations  rather  than  sound.  The  simplest  model  problem  that 
exhibits  structural  damping  by  this  means  is  that  in  which  a 
bending  wave  on  a  thin  elastic  plate  is  reflected  at  a  trailing  edge 
in  a  parallel  mean  flow  (see  Figure  1)  [13].  Incident  bending 
wave  energy  is  dissipated  at  the  edge  by  scattering  into  sound, 
and  by  a  transfer  to  the  kinetic  energy  of  the  mean  flow  and 
wake.  Below  the  coincidence  frequency  the  efficiency  with 
which  sound  is  generated  is  usually  very  small,  but  much  larger 
net  losses  are  predicted  owing  to  the  generation  of  wake 
vorticity. 

The  predicted  fractional  bending  wave  power  dissipated  at 
the  edge  is  plotted  in  Figure  2(a)  for  a  steel  plate  in  water  for 
several  values  of  the  mean  flow  Mach  number  M  when  the  edge 
of  the  plate  is  "free".  In  this  figure  the  radian  frequency  w  is 
normalized  by  the  coincidence  frequency  w  .  For  M  =  0  the 

c 


153 


Report  No.  AM-96-006 


Boston  University,  College  of  Engineering 


Figure  1. 


dissipation  can  be  attributed  to  sound  radiated  from  the  edge. 
When  M  >  0  there  is  a  significant  increase  in  power  loss  at  low 
frequencies  caused  by  a  transfer  of  energy  to  the  mean  flow  and 
wake,  with  negligible  changes  in  the  radiated  sound  power  (at 
least  for  M  <  1).  This  is  clear  from  Figure  2(b),  which 
compares  the  total  power  loss  for  a  steel  plate  in  water  (solid 
curve)  with  the  acoustic  power  (dashed  curve)  when  M  =  0.01. 
Inspection  of  Figures  2(a)  and  (b)  confirms  that  the  sound  power 
is  hardly  influenced  by  the  flow.  At  high  frequencies  the 
radiation  of  sound  accounts  for  all  of  the  dissipation;  this  is 
because,  as  o)  ^  co ,  bending  wave  energy  is  contained 
principally  in  the  evanescent  motions  in  the  fluid  on  either  side 
of  the  plate,  that  propagate  with  negligible  plate  motion  at  a 
velocity  which  is  only  slightly  less  than  the  speed  of  sound. 

The  dotted  curve  in  Figure  2(b)  is  the  fraction  of  the 
dissipated  energy  which  appears  as  kinetic  energy  of  the  vortex 
wake.  This  indicates  that,  over  the  range  of  frequencies  within 
which  acoustic  dissipation  is  negligible,  roughly  half  of  the 
dissipated  energy  is  used  to  accelerate  the  mean  flow  and  half 
goes  into  the  wake;  this  is  not  generally  the  case,  but  depends  on 
the  type  of  mechanical  constraints  at  the  edge  [13]. 


THE  BIAS-FLOW  ELASTIC  PLATE 

Consider  an  elastic  plate  perforated  with  small  circular 
apertures.  Let  there  be  N  apertures  per  unit  area  each  with 
radius  R.  A  mean  pressure  differential  is  maintained  across  the 
plate  which  produces  nominally  steady,  low  Mach  number,  high 
Reynolds  number  jets  of  fluid  through  the  apertures  (Figure  3). 
Generalized  forms  of  the  thin  plate  bending  wave  equation  have 
been  derived  in  [14,  15]  respectively  for  bias-flow  and 
grazing-flow  screens  when  the  length  scale  of  the  motion  is 
large  compared  to  both  the  aperture  radius  R  and  the  distance 
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Figure  4. 


between  the  apertures.  For  the  bias-flow  case  it  was  shown  in 
[14]  that  long  wavelength,  small  amplitude,  flexural  displace¬ 
ments  f  (proportional  to  e’'""* )  of  a  plate  of  bending  stiffness  B 
and  mass  m  per  unit  area  immersed  in  fluid  of  mean  density  p, 
satisfy 

{[l-2aa/(l-cr)]BV‘  -  mw  }f 

+  {1  +  2NRKj^[l  -  (l/2Rpw")(BVW)]}[p]  =  0,  (1) 

where  a  =  is  the  fractional  open  area,  a  is  Poisson’s  ratio 
for  the  plate,  V‘'  3  (where  the  z-  or  xy-plane 

coincides  with  the  undisturbed  location  of  the  plate),  and  is 
the  Rayleigh  conductivity  of  the  apertures,  [p]  is  the  net  pressure 
force  on  the  plate  (in  the  direction  opposed  to  that  of  increasing 

f). 

The  conductivity  determines  the  volume  flux  Q  (in  the 
positive  f-direction)  through  an  aperture  as  a  result  of  a  long 
wavelength  pressure  differential  [p]  according  to  the  relation, 

Q  =  KJ<I]  =  KJpJ/ipo),  (2) 

where  ^  is  the  corresponding  perturbation  in  velocity  potential. 
In  this  limit,  in  which  the  length  scale  of  the  motion  of  the  plate 
is  much  larger  than  the  aperture  radius  R,  it  was  argued  in  [14] 
that  K  is  well  approximated  by  its  value  for  a  rigid  plate. 
has  the  dimensions  of  length  and  is  equal  to  2R  in  a  stationary 
ideal  fluid,  but  becomes  complex  valued  and  frequency  depend¬ 
ent  in  the  presence  of  flow,  because  of  the  unsteady  shedding  of 
vorticity  from  the  aperture  edges. 

Introducing  the  Strouhal  number  S  =  coR/U,  where  U  de¬ 
notes  the  vorticity  convection  velocity  in  the  bias  flow  jets 
(approximately  one  half  of  the  mean  jet  velocity),  we  can  write 

=  2R{r(S)-iA(S)},  (3) 

where  T  and  A  are  real  and  positive  (for  S  >  0).  The 


dependence  of  F  and  A  on  S  is  shown  in  Figure  4  (from  [14]). 


Dissipation  of  sound 

Let  a  plane  acoustic  pressure  wave  of  the  form  exp(-i/cz)  be 
incident  normally  on  the  screen  from  z  >  0  (Figure  3),  where  k 
=  oj/c  (c  =  speed  of  sound)  is  the  acoustic  wavenumber.  The 
pressure  can  be  taken  in  the  form 

p  =  exp(-i/cz)  -h  9?exp(ixz)  ,  z  >  0, 

=  Qfexp(-i/cz),  z  <  0,  (4) 

where  9?  and  O'  are  appropriate  reflection  and  transmission 
coefficients,  respectively.  The  displacement  f  of  the  screen  is 
related  to  p  by  f  =  (1/pw  )ap/5z  (z  ^  ±0).  This  together  with 
equation  (1)  supply  two  equations  from  which  9?  and  O'  can  be 
determined.  The  absorption  coefficient  5  =  l-|9?p  ’[Q'M 
is  equal  to  the  fraction  of  the  incident  acoustic  power  dissipated 
at  the  screen.  The  solid  curves  in  Figure  5  show  the  variation 
of  5  with  Strouhal  number  for  sound  incident  on  a  perforated 
aluminium  plate  in  air,  when  the  bias  flow  convection  Mach 
number  M  =  0.05  and  R/h  =  5,  where  h  is  the  thickness  of  the 
plate.  The  maximum  attenuation  of  5  =  0.5  (3  dB)  occurs  at  an 
open  area  ratio  a  =  M  over  a  range  of  Strouhal  numbers 
centered  on  S  «  0.4.  At  very  low  frequencies  the  impedance 
offered  to  the  sound  by  the  plate  becomes  negligible,  the  bias 
flow  jet  velocities  are  unmodulated  by  the  sound,  and  there  is  no 
transfer  of  acoustic  energy  to  hydrodynamic  motions. 

The  corresponding  predicted  attenuations  for  a  rigid  bias- 
flow  screen  are  depicted  by  the  broken  curves  in  Figure  5.  The 
elastic  and  rigid  screen  predictions  differ  when  the  frequency  is 
low  enough  for  the  screen  to  be  regarded  as  acoustically 
transparent.  The  illustrated  behavior  is  typical  of  cases  in  which 
the  fluid  loading  is  small  (for  metallic  screens  in  air,  say).  For 
the  other  extreme  of  a  steel  plate  in  water  the  absorption  of 
acoustic  energy  is  found  to  be  negligible  at  all  frequencies. 
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Figure  6. 


Figure  7. 


Damping  of  bending  waves 

A  flexural  wave  f  ^  with  the  associated  pressure 
loading 

[p]  =  -2p{«^/V(k^-K^)}^e''“ 

satisfies  the  bias-flow  bending  wave  equation  (1)  provided  the 
bending  wavenumber  k  is  a  zero  of  the  dispersion  function 

D(k,w)  =  [l-2aa/(l-(r)]Blc'  -  mu 

-  {2pa)VV(k2-/c^)}{l  +  2NRK^[l-(BIc‘'-mJ)/2RpJ]}  (5) 


In  the  absence  of  the  bias  flow  (K^  2R),  D(k,aj)  has 

precisely  two  equal  and  opposite  real  zeros.  They  satisfy 
I  k  I  >  K  and  correspond  to  undamped  flexural  waves  on  the 
perforated  plate.  To  the  present  order  of  approximation  there  is 
no  damping  due  to  scattering  of  structural  energy  into  sound. 

The  absorption  of  bending  wave  energy  by  the  bias  flow  jets 
causes  the  formerly  real  bending  wavenumber  k  to  acquire  a 
small  positive  imaginary  part  that  accounts  for  the  progressive 
decay  of  the  wave.  The  power  absorbed  per  wavelength  of 
propagation  is  plotted  against  the  Strouhal  number  S  =  ojR/U  in 
Figures  6  and  7  for  a  perforated  steel  plate  in  water  and  several 
different  open  area  ratios  ct  for  an  aperture  radius  R  =  h,  3h 


Figure  8. 


Figure  9. 
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Figure  10.  Figure  11. 

respectively.  The  bias  flow  vorticity  convection  velocity  U  =  3  F  =  2.29S^,  A  ~  1.70S{1  -  1.67S  },  S  ^  1; 

m/s  (M  -  0.002).  Increasing  the  open  area  ratio  increases  the 

attenuation,  whereas  increasing  the  size  of  the  apertures  reduces  F  «  1,  A  «  l/^rS,  S  >  1,  S  =  ojR/U.  (7) 

the  maximum  possible  attenuation. 

To  assess  the  significance  of  these  predictions  we  plot  in  In  applications  to  the  attenuation  of  flexural  and  acoustic  waves, 

Figure  8  the  corresponding  absorption  of  bending  waves  on  a  the  most  significant  effects  occur  at  low  Strouhal  numbers  S  < 

coated  steel  plate  in  water.  For  the  purpose  of  illustration  the  1,  for  which  the  following  rational  approximations  are  adequate 

mass  of  the  coating  is  neglected,  but  the  coating  is  assumed  to  for  estimating  dissipation  rates: 

make  the  effective  stiffness  of  the  plate  complex,  such  that  the 

bending  stiffness  B  of  the  uncoated  plate  is  replaced  by  F  ~  2.29SV(1  4*2.2982),  ^  1. 78/(1  +  1.6782).  (g) 

B(1  -  iri),  where  t)  (>  0)  is  a  loss  factor  whose  value  is 

determined  by  the  properties  of  the  coating,  and  is  typically  of  The  variations  of  F  and  A  given  by  these  approximations  are 

the  same  order  as  shown  in  the  figure  [16].  As  oj  increases  similar  to  those  already  depicted  in  Figure  4  for  the  bias-flow 

towards  a?  both  the  bending  and  acoustic  wavelengths  become  screen.  The  rational  approximation  for  A  is  too  large  by  a 

comparabfe  to  the  thickness  of  the  plate  and  thin  plate  theory  factor  of  3  as  S  oo ,  but  the  behavior  for  8  <  1  should  be  well 

ceases  to  be  strictly  applicable.  A  comparison  with  figures  6  approximated  by  (8). 

and  7  indicates  that  the  attenuations  achieved  by  the  bias  flow  The  flexural  wave  dispersion  equation  is  now 

jets  can  be  comparable  with  those  predicted  for  the  coated  plate.  42  222 

Much  smaller  attenuations  are  predicted  for  bending  waves  [l-2a(T/(l-cr)]Bk  -  mo)  -  {2p(a)-Uk)  /V(k  -k  )} 
propagating  on  a  bias  flow  screen  in  air  (see  [14]  for  details). 

X{1  -H  2NRK^[l-(Bk"W)/2Rp(a;-Uk/]}  =  0.  (9) 

THE  GRAZING-FLOW  SCREEN  Proceeding  as  for  the  bias-flow  screen  we  calculate  the 

The  modified  bending  wave  equation  for  a  perforated  plate  in  typical  attenuations  illustrated  in  Figures  10  and  11  for  a  steel 
a  uniform,  low  Mach  number  tangential  mean  flow  at  speed  U  plate  in  water  at  different  fractional  open  areas  oc  and  aperture 

in  the  x-dhection  (Figure  9)  is  given  in  [15]  in  the  form  radii.  The  results  are  for  bending  waves  propagating  in  die 

direction  of  the  mean  flow  (solid  curves),  and  in  the  direction 
{cj+iUa/ax}  ^{[(l-2Q:(T/(l-(7))BV^-mJ]f  4-  (1 -K2NRK  )[p]}  opposite  to  the  mean  flow  (dashed  curves). 

Comparing  these  predictions  with  Figure  8  for  a  nominally 
-  (NK  /p)(BV^-mw  )[p]  =  0.  (6)  homogeneous  plate  with  one  or  more  layers  of  absorptive 

^  coatings,  it  can  be  seen  that  the  attenuations  are  comparable  with 

The  following  asymptotic  formulae  for  K^/2R  =  F-iA  are  the  best  obtainable  by  coating  the  plate,  at  least  for  8trouhal 

derived  in  [17]:  numbers  in  the  interval  0.1  -  1. 
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CONCLUSION 

The  predictions  presented  in  this  paper  for  bias-  and 
grazing-flow  perforated  elastic  plates  indicate  that  significant 
attenuations  of  sound  and  flexural  vibrations  are  possible 
provided  a  sufficiently  large  fluctuating  pressure  gradient  can  be 
established  across  the  screen  to  cause  unsteady  production  of 
vorticity  in  the  apertures.  When  the  fluid  loading  is  large  (e.g., 
for  metallic  plates  in  water)  an  unrestrained  elastic  screen  is 
effectively  acoustically  "transparent",  and  the  damping  of  sound 
by  vorticity  production  tends  to  be  negligible.  Large  acoustic 
attenuations  can  occur  for  a  lightly  loaded  plate  (in  air,  say). 

The  damping  of  bending  waves  increases  with  the  fluid 
loading,  since  the  surface  pressure  fluctuations  produced  by  a 
wave  of  given  surface  displacement  tends  to  increase  with  fluid 
density.  For  bias-flow  screens  the  damping  is  greatest  for 
aperture  Strouhal  numbers  S  =  wR/U  (based  on  the  vorticity 
convection  velocity  U)  between  about  0.03  and  3,  and  the 
attenuation  experienced  by  bending  waves  on  a  perforated  steel 
plate  in  water  can  then  exceed  or  be  comparable  with  that  for 
waves  on  a  heavily  coated  plate.  Similar  results  are  obtained  for 
grazing-flow  screens,  the  maximum  attenuation  for  waves  on  a 
perforated  steel  plate  in  water  occurs  for  wR/U  in  the  range  0.1  - 
1  (where  U  now  denotes  the  tangential  mean  velocity). 
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